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AN INVESTIGATION OF THE UTITUDE OF UKIAH, CALIF., 

AND OF THE MOTION OF THE POLE. 



By Walter D. Lambert, Mathematician^ U, S. Coast and Geodetic Survey, 



INTRODUCTION. 

The following discussion of the observations of the International 
Latitude Service was occasioned by an article of Prof. Andrew C. 
Lawson,* in which he called attention to the fact that the latitude 
of Ukiah Observatory, Califomia, maintained by the U. S. Coast and 
Geodetic Survey for the International Latitude Service, appears to 
be steadily increasing at the rate of about OfOl a year. The latitude 
of Lick Observatory appears to be increasing even more rapidly, 
although the progressive increase is compHcated by an abrupt 
decrease of 0763 that seems to have occurred in the latter part of 
the year 1903. Obviously if the existence of progressive changes* 
of this order of magnitude can be established beyond question, the 
study of them becomes of great interest to geophysicists, astronomers, 
and geodesists. 

The purpose of this paper is the discussion of the mass of data 
accummated by the observatories of the International Latitude 
Service, in order to ascertain how far the large progressive change of 
latitude is peculiar to Ukiah alone among all the stations of the Lati- 
tude Service and how far the apparent change of latitude may be 
accounted for in other ways than by an actusu shifting of the upper 
strata of the earth^s crust in the vicinity of the station, this being 
Prof. Lawson's explanation. In order to reach any assured con- 
clusion whatever, it is evidently necessary to compare Ukiah with 
other stations; accordingly all the latitudes of all stations of the 
International Service are here analyzed and discussed. The method 
of observation and the program of stars to be observed are the same 
at all stations of the International Service. The observations of the 
service therefore form a homogeneous whole peculiarly adapted to 
discussions of delicate questions of this sort; for one thmg the errors 
of the star places have comparatively little eflfect on the results. 
The observations at Lick Observatory nave not this advantage. No 
serious attempt, therefore, is here made to discuss the evidence from 
the Lick observations except to reproduce the passage on this subject 
from Prof. Lawson^s article, together with some corrections and 
comments, chiefly those kiadly supplied by the director of Lick 
Observatory. 



1 " The mobility of the Coast Ranges of CaUfomiai an exploitation of the elastic rebound theory;" Uni- 
versity of California Publicaticms, Bull. Departaient of Geology, vol. 12. No. 7 Jan. 11, 1921. 

> The word "progressive" is here used in much the same sense that the word "secular" is used in the 
theory of astrononucal perturbations; that is, a progressive chance is one proportional to the time, or nearly 
so. The word "secular," however, carries the connotation of almost inaefluite continuance, as measured 
by ordinary human standards. The changes in latitude with which this report deals have not been proved 
to be of this character; the word ''secular" has thOTefore been avoided and ''progressive" substituted. 
(See p. 37.) 
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In the harmonic analysis of the latitudes the series used are neces- 
sarily short when measured by the number of periods used; the 
harmonic constants of anjr one component derived from such a short 
series are considerably affected by tne presence of other components. 
The approximate formulas of the type used by this Survey in the 
ordinary harmonic analysis of short-period tides for clearing one 
component of the effects of another are not quite adequate in the pres- 
ent mstance and a more accurate process is explained and used. It 
resembles the procedure given by Darwin for clearing one lon^-period 
tidal component from tne effects of the others but is carried one 
step further by the solution of a set of linear equations in general 
terms. The derivation of the final formulas for numerical computa- 
tion is more laborious than for the approximate formulas for dealing 
with short-period tides, but the numerical calculation with the more 
exact formulas, once they are obtained, is quite as simple as with 
the approximate ones. 

In order not to interrupt the main discussion with matters pri^- 
marily mathematical, the development of the formulas for clearing 
one harmonic component from the effects of others and some other 
mathematical developments are given separately in Chapter VI. 
This separation has necessitated a certain amount of repetition. On 
accoimt of the general interest of the subject of latitude variation, 
the attempt has been made to render the details intelligible to a 
fairly wide circle of possible readers. For this reason some matters 
have been treated at greater length than would be necessary in a 
book intended solely for specialists. 

Some results not directly connected with the problem in hand are 
given in CJhapter IV. They were found to be deducible from results 
already obtamed with such a comparatively small amount of extra 
labor that it seemed worth while to obtain and present them. 

Miss Sarah Beall, of the division of geodesy, has given valuable help 
in plotting and reading curves and in malang or checking much of 
the large mass of computation required in this work. Without her 
assistance the appearance of this investigation would have been much 
delayed. I. N. BeaU and M. A. Crews, of the division of geodesy, 
particularly the latter, have also given extensive help with the com- 
putations. 

While this publication was in the printer's hands the author under- 
took a revised presentation of the subject of polar motion with the 
Jurpose, among others, of meeting certain objections raised bj 
roiessor Schlesinger.* As a result it seemed justifiable to omit 
Tschardjui entirely from the determination of the polar motion. The 
rate of displacement of the pole from 1900 to 1917, inclusive, is then 
increased from (y/0050, as foimd in chapter II, section 5 (p. 35), 
to 0f0062 annually and the average direction of the motion of the 
North Pole comes out nearly along the meridian of 90° west. These 
figures should be considered as a revision of those contained in the 
present pubUcation. For details see the author's paper *^The inter- 

f)retation of apparent changes of mean latitude at tne international 
atitude stations," which will appear in the Astronomical Journal, 
Vol. 34, no. 804. The paper contains also a brief discussion of earth 
movements and of the elastic-rebound theory of earthquakes. 

• F. Schleisinger, On Progressive Changes of Latitude, Astronomical Journal, vol. 34, p. 42. 1922. 



Chapter I.— AN EXAMINATION OF PROF. LAWSON'S DATA. 

The data from which Prof. Lawson deduced the progressive increase 
in the latitude of Ukiah are found in an article by Su* F. W. Dyson, 
Astronomer Royal,^ and consist of curves showing the variation of the 
latitude after the annual portion ^ of the variation has been eliminated 
by calculation. The abscissa used in plotting the curve is the date; 
the ordinate, which we shall denote by A</>, is the difference between 
the observed latitude corrected for the annual variation and a certain 
initial latitude, which, for present purposes, may be quite arbitrary. 
The change in latitude still left in the curves is then due to the free 
oscillation of the earth, except in so far as the elimination of the 
annual portion may have been imperfect, or as small effects of other 
motions with periods unknown or uncertain may remain also. The 
curves are given for three stations only — Ukiah, Mizusawa, and Car- 
loforte. The abscissas and ordinates from which the curves were 
plotted are not given in full. The data given for all stations consist 
of the epochs of maximiun and minimum latitude and the differences 
between each maximum and the two minima adjacent; this is not 
sufficient to reconstruct the curves for Tschardjui, Cincinnati, and 
Gaithersburg, so all consideration of these stations is deferred until 
the next chapter. 

The curves shown may be considered as somewhat irre^ar sine 
ciu'ves with variable amplitudes of oscillation. In the Ukiah curve 
an algebraic increase with the time of both the maximum and mini- 
mum ordinates is pretty evident. Prof. Lawson obtained the rate of 
progressive^ increase by drawing a straigjht line as nearly as possible 
through the successive maxima and a similar line through the minima. 
The average slope of the two lines was taken as representing the pro- 

Eressive rate of change of the latitude of Ukiah. For this rate rrof . 
awson finds* 0.29 meter, or 0?0094 second of latitude per year. 
He then continues: 

In arriving at this figure it is to be noted that, since the mean position of the pole 
represented by the zero line is based upon observations of all stations including Ukiah, 
the zero line adopted diverges nbrtheny from the true zero line, and the rate of creep 
inferred from the curve, 0.29 meter per year, is in reality somewhat less than the 
true rate. For, if Ukiah be moving northerly, the increase in latitude due to this 
cause should not enter into a determination of the mean position of the pole, that is 
the position of the zero line in relation to the curve. 

The author doubts whether the statement quoted correctly inter- 
prets Sir F. W. Dyson's procedure. The only statement in regard to 
this is*^ — 

§6. The results for the six observatories given in A. N. 4841 were next examined 
separately. The mean annual terms for the interval 1900.0— 1912.0 are: • * * * 

» Monthly Notices of the Royal Astronomical Society, vol. 78 (1918), p. 462. 

s The annual portion of the variation of latitude includes the annual part both of the motion of the 
pole itself and of all other changes in the latitude, which like the Kimura effect, have an annual periodic 
pcvtion. 

» For the sense in which the word "progressive" is used, see footnote No. 2, p. 1. 

4 See p. 436 of the article referred to on p. 1. 

' L/OC cit ante lo 467 

• These values tabulated by Dyson, but here omitted, are in substantial accord with those found by an 
independent discussion and given on p. 61 of this publication. The exact values are not material to the 
que^ion immediately under discussion. 

3 
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§7. In the discussion of the free period as shown by the different observatorieB the 
annual terms given at the beginning of §6 were extracted and the figures then plotted. 
This is a rather complicated diagram to reproduce, and only the three observatories 
at Carloforte, Mizusawa, and Ukiah are given. 

A reference to Astronomische Nachrichten No. 4841/ which is a long 
article by Przbyllok, shows that the stations are there tabulated 
separately and that the A<^'s for any station are the differences be- 
tween the observed latitudes and some invariable initial latitude, the 
latter having nothing to do with the motion of the pole or with the 
latitudes of the other stations. 

It might be objected, however, that the motion of the pole does 
enter to a certain extent into Przbyllok's A</>'s, for the so-called 
observed latitudes really depend both on instrumental readirigs and 
on the declinations used. The declinations of Przbyllok's demiitiye 
A</>'s are corrected declinations, the corrections deduced from a dis- 
cussion of the differences between the observed latitudes (with the 
original uncorrected declinations) and the adjusted latitudes of the 
International Latitude Service, which depend on all stations. Thus 
to a certain extent Przbyllok's observed latitudes and hence his A<^'s 
for anv station, which are presumably Dyson's also, depend on the 
latitud.es of the other stations. An examination of the corrections 
actually used, as given at the foot of page 282 of Astronomische 
Nachrichten No. 4841, shows, however, that the mean correction to 
the declination is considerably less than 0?001. Furthermore the 
corrections are the same for the earljr part of the series as for the 
later parts and so have no effect in increasing or diminishing any 
apparent progressive change of latitude such as occm« at Ukiah. 

Dyson gives results extending beyond Przbyllok's tabulation, which 
does not go beyond 1912.0. Presumably Dyson's results are based 
on the provisional values in the Astronomische Nachrichten,® but 
whether any correction was applied to them to make them homoge- 
neous with Przbyllok's results Dyson does not state. A lack of 
complete homogeneity would not be particularly important in Dyson's 
work, the purpose oi which was to obtain an accurate value for the 
length of the 14-month period and which was not intended for the 
use to which Prof. Lawson has put it. 

It appears then that neither JProf. Lawson's change of 0.29 meter 
or 0''0094 a year nor the Oi'OOSl presently to be deduced from the 
same data by another method should be augmented for any supposed 
change in the zero line of the A</>'s. 

There may be some objection to the method of using only the 
maxima and minima of the latitude-variation curve in order to aeter- 
mine the rate of progressive change," and it seems certain that the 
method of taking ordmates of the curve at small, equal intervals and 
finding the mean willgive a more reUable result. It seemed desirable 
also to see how far Carloforte and Mizusawa might resemble Ukiah 
in showing a progressive change of latitude. Accordingly, the curves 
of all three stations were read, at intervals of one montn — ^i. e., one- 
twelfth of a year — and the means were taken over a period of 14 

7 Vol. 202 (1916). The article is a careful discussion of the observations of the International Latitude 
Service. 

* For references, see p. 8. 

* For example, the form of a curve representing a periodic fimction mav be such that the average of 
maximum and Tninimntn ordinates differs from the mean of all ordinates; thus in careful tidal work a dis- 
tinction must be made between mean sea level , which is found from the mean of closely spaced , equidistant 
ordinates of the curve dra^Ti by an automatic tide gauge, and half-tide level, which is halfway between 
the mean ordinate of the high waters and the mean ordinate of the low waters. 
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months, or more specifically, 15 consecutive ordinates were used and 
half weight only was assigned to the first and to the fifteenth ordinate. 
The period used (14 months = 426.1 days) is so nearly the free period 
of the latitude variation that the mean over the 14 months should 
be practically unaffected by periodic changes of phase that run their 
course in the free period or its submultiples. For example, if we 
write the principal part of the free oscillation of the latitude in the 
form A cos (k^— a), where A and a are assumed to be constants, i the 
time, and k a constant such that /c^ = 27r when i = 432.5 days, then the 
maximum error due to the incomplete elimination of this term from 
the mean of 15 readings taken as above is easily seen by formula (4) 
(p. 71) to be less than 0.015 A, Now A is 0?20 or less, so the system- 
atic error is at most about 0^003 for any 14-month period. 

The following table shows the means of the values of A</> read from 
the curve for the period the middle of which is shown in the first 
two columns: 

Table 1. — Observed changes in latitude. 



Middle of period. 



Year. 



1900 
1901 
1902 
1904 
1905 

1906 
1907 
1908 
1909 
1911 

1912 
1913 
1914 
1915 



Month. 



7 
9 
11 
1 
3 

5 
7 
9 
11 
1 

3 
5 
7 
9 



Mean value of A<^ for— 



Mizu- 
sawa. 



-0.010 

- .011 

- .042 

- .023 

- .040 

- .012 

- .025 
+ .004 

- .017 

- .047 

- .075 

- .088 

- .034 

- .032 



Carlo- 
forte. 



// 

+0.004 
+ .033 
+ .004 
+ .022 
+ .029 

+ .021 
+ .001 
+ .002 
+ .016 
+ .028 

+ .013 
+ .032 
+ .032 
+ .037 



Ukiah. 



// 



+0.001 
- .040 
+ .001 
+ .047 
+ .054 

+ .068 
+ .066 
+ .047 
+ .058 
+ .070 

+ .063 
+ .065 
+ .118 
+ .152 



The months shown are elapsed months, or more accurately elapsed 
twelfths of a Juhan year; thus in the first Une the '^ 7 '' indicates that 
seven-twelfths of the year 1900 have elapsed, or the date is about 
August 1. 

The values of A<^ in the table were plotted as ordinates against the 
time as abscissa and the results are shown in figures 1, 2, and 3. The 
attempt was then made to fit a straight line to the plotted points by 
the method of least squares; this corresponds to an assumed uniform 
rate of variation of the latitude. 

Each mean value of A<^ gives an observation equation of the form 

A(^ = x + ^2/, (1) 

X being the adjusted value of A<^ at the time 1908 + 9 months and i 
the time reckoned from this epoch; y is the adjusted rate of change 
of the latitude. From the adjustment there was found with the 
14-month period as the unit of time the following values, which are 
given with the probable errors attached : 

For Mizusawa a; = - 0f0338 ± 0^0042 y=- 0?0029 ± 0?0010, 
for Carlof orte a; = + 0.0202 ± 0.0022 y = + 0.0013 ± 0.0005, 
for Ukiah « = + 0.0590 ± 0.0048 ~ y = + 0.0095 ± 0.0012. 



U. S. COASt AND GEODETIC 8UEVBY. 



Fia. 1.— ProgreaalTe chaage of latitude at Ulzusawa bum DjBon'a i 



of latitude at Carlaforte from Dyson's c 



Fia. 3,—Prognealye change otiatllude al UkUli tram Dyson's curves. 
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To reduce the values of y and their probable errors to the year as 
unit they must be multipUed by -J^; we thus find: 

For Mizusawa y = - 070025 ± Of 0009, 
for Carloforte y = + 0.001 1 ± 0.0004, 
for Ukiah y = + 0.0081 ± 0.0010. 

The annual rate of change f oimd for Ukiah is thus slightly smaller 
tlian the rate (0f0094) found by Prof. Lawson and this latter rate, 
as we have seen, should not be mcreased for any supposed change in 
the zero from which A</> is reckoned. Even the smaller rate of change 
for Ukiah is, however, quite large enough to excite interest as being 
rather too large to be due to accidental error. Any attempt to 
explain this rate is, however, postponed until all the six stations of 
the International Latitude Service have been examined for similar 
rates. (See pp. 21-31.) 



Chapter IL— THE STATIONS OP THE INTERNATIONAL 

LATITUDE SERVICE. 

SECTION 1. INTI^ODUCTOBY MATTER. 

The stations of the International Latitude Service and their distri- 
bution in longitude are shown in the following table. West longi- 
tudes are positive. 

Table 2. — Stations of the International Latitude Service. 



station. 



Mizusawa. . . 
Tschardjui » . 

Carloforte 

Gaithersburg 
Cincinnati. . . 
Ukiah 



General locality. 



Japan 

Turkestan 
Sardinia... 
Maryland. 

Ohio 

California. 



Longitude. 



-141 08 

- 63 29 

- 8 19 
+ 77 12 
+ 84 25 
+123 13 



Observations 

ceased at end 

ofyBar — 



1914. 

1914. 
1915. 



1 This is the spelling adopted in the publications of the International Latitude Service. The foUowing 
forms have been found in works of reference published in Great Britain and the United States: Charchm 
Tschaijui, Chaijooee, Chardshuil , and Chaijui. For an English-speaking person the last-given form seems 
to be an adequate transliteration. The longitude is that of the old station. The longitude of the new 
station is -63* 35'. < 

The methods adopted for the discussion of the observations at these 
six stations are essentially the same as those of the preceding chapter 
but with some differences in detail. First of all the observed lati- 
tudes were plotted as ordinates against the dates expressed in years 
and decimals of a year as abscissas. The sources of the observ«ed 
latitudes were the publications of the International Latitude Service 
and the Astrononusche Nachrichten.^ Pages 210 to 219 of the Re- 
sultate, Vol. Ill, and pp. 169-178 of Vol. V, contain the definitive 
results for the years 1900 to 1911, inclusive. The declination system 
of both volumes is that of Volume III, which is maintained for the 
sake of uniformity and convenience.^ TJie same declination system 
is used in the provisional results given in the Astronomische Nachrich- 
ten.* Some of the star pairs have, however, been dropped as no 
longer available on accoimt of their precessional motion, and have 
been replaced by other pairs. In working up the provisional results 
these new pairs were not used, and the results are therefore based 
on the original pairs that still remain available, amounting to about 
five-sixths of the complete program. The results for Tschardjui after 
July 23, 1909, when the new observatory was occupied, were reduced 
to the old observatory by applying to the latitudes observed after 
that date a constant correction of — 05225, as given on page 3 of the 
Resultate (Vol. V) . 

> Resultate des Intemationalen Breitendienstes: Vol. ni, by Albrecht and Wanach, Berlin. 1909; vol. 
V.by Wanach, Berlin, 1916. Astronomische Nachrichtem voL 198, No. 4749; vol. 201, No. 4802; vol. 203, 
No. 4858; vol. 205, No. 4906; Y(A, 208, No. 4969. 

» Vol. V, p. 76. 

* See Astronomische Nachrichten No. 4749, vol. 198. 

8 
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The plotted latitudes for each station were then adjusted graphi- 
cally by drawing a curve to conform as nearly as might be to the 
plotted points without showing sharp vertices or rapia alternations 
in slope. On these curves for the several stations all the subseauent 
work is based; they were therefore drawn with some care and the 
first draft of a curve was always revised by one other person and 
sometimes by two others. The observations are consistent enough 
so that there was very seldom room for any serious difference m 
judgment as to the proper form and location of the curve. In doubt- 
ful cases the doubt was decided in favor of smoothing out the obser- 
vations too little rather than too much. 

SECTION 2. HARMONIC ANALYSES.^ 

Harmonic analyses were made of the curves to detenmne expres- 
sions for the annual and the 14-month components at each station. 
The period of the annual component was assumed to be one Juhan 
year, or 365.25 days. No distinction was made between leap years 
and common years. The period of the 14-month component was 
assumed to be 432.5 mean solar days; this is a mean between Dyson's 
result of 432.2 days® and Wanach's of 432.8 days.' This gives the 
14-month, or free, period of the latitude variation as 1.18412 Julian 
years and the change of phase of the free oscillation in one month 
(one-twelfth of a Jmian year) as 25?33526 and in one vear 304?0231. 
One year equals 0.844509 of the free period and the change of phase 
of the annual component in one-twelfth of a free period (0.098677 
year) is 35?5236, and the change of phase in one free period is 
426?2834. These numbers are here set down for reference as needed. 
When the 14-month component or a 14-month period is spoken of 
hereafter in this report, this should be imderstood as an abbreviation 
for a 432.5-day component, etc., unless it is clear from the context 
that exactly ^^ of a Julian year is meant. 

For the annual component the curve for a station was read at 
intervals of one-twelfth of a year, the first reading coinciding with 
the beginning of a year. For the 14-month component, represent- 
ing the free oscillation, the period was also divided into 12 parts, 
the interval between readings being 0.0987 year, as above; the first 
reading of the series was always taken at the beginning of some 
year. Since we have approximately 5 free periods equal to 6 annual 
periods and 6 free perioos equal to 7 annual periods, with about the 
same degree of approximation in both cases, the obvious lengths for 
a series to be analyzed harmonically are 6 or 7 years. Both lengths 
are here used. A 6-year series for the annual component includes 
72 readings, the time between the first and the last reading being 
5.92 years. For the 14-month component a 6-year series is taken to 
mean one that represents 5 free periods; it has 60 readings with an 
interval of 5.82 years between the first and the last reading. Simi- 
larly, a 7-year series has 84 readings for the annual component, with 
6.92 years between the first and the last reading, and 72 readings for 
the 14-month component, with 7.01 years between the first and the 
last reading. 

» See sec. 2 of Chap. VI, p. 73. 

> Monthly Notices Royal Astronomical Society, vol. 78 (1918), p. 402. 

1 Resultate, Vol. V, p. 209. 
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The mean reading correspondin|g to each twelfth of a period 'was 
found and these means were subjected to the usual process of the 
harmonic analysis for 12 ordinates. The values of a and c in the cus- 
tomary notation, say for the annual component (see p. 76), are 
affected by the presence of the 14-month component, and vice versa. 
To remove the effects, the formulas developed onpages 78-79 were used. 
To show the order of magnitude of the corrections so obtained the 
table of harmonic constants that follows gives the results both from 
the uncorrected and the corrected s's and c's. The epochs or f 's 
for the 14-month component are all reduced to the beginning of the 
year 1900. The quantity f -t-X, where X is the west longitude of the 
station, is added to show how nearly the motion of the pole is uni- 
form and circular. If the motion were uniform and circular and in 
the same direction as the rotation of the earth and if all the results 
were errorless, the quantity f + X would be constant for all longi- 
tudes. The length of the series analyzed is to be inferred from^ne 
dates iQ the second column; thus the dates 1900-05 indicate a 
6-year series with the first reading for both components on 1900.00. 
The dates 1900-06 indicate a 7-year series beginning as before- the 
final reading for the 14-month component of such a series falls on 
1907.01. The harmonic constants in the table below are not quite 
the definitive values resulting from this investigation (see p. 51) but 
they are accurate enough for the immediate purpose. 

Table 3. — Harmonic constants for the variation of latitude. 

[Assumed form, R cos {Kt-i).\ 
ANNUAL COMPONENT. 







Before clearing. 


After clearing. 


zm 




Years of 
series 
inclu- 
sive. 












Station and longitude. 


Ampli- 
tude. 


Epoch. 


Ampli- 
tude. 


Epoch. 


f+X 






R 


i 


R 


f 








n 


e 


tt 


e 


• 


Mizusawa. Janan. X-»— 141*08' 


1900-05 
1900-06 


0.113 
.097 


13.3 
16.1 


0.106 
.107 


15.7 
15.9 


234 6 




234.8 




1906-11 


.138 


4.5 


.132 


10.4 


229.3 




1906-12 


.110 


14.4 


.123 


9.2 


228.1 




1910-15 


.127 


18.9 


.140 


17.4 


236.3 




1911-17 


.116 


20.2 


.111 


13.1 


232.0 




1912-17 


.105 


1.0 


.100 


7.1 


226.0 


Tschardjui, Turkestan. X— -6^'29' 


1900-05 


.104 


334.7 


.101 


340.0 


276.6 




1900-06 


.098 


340.0 


.105 


335.9 


272.4 




1906-11 


.136 


315.6 


.133 


322.3 


258.8 




1906-12 


.122 


331.3 


.131 


323.7 


260.2 




1908-14 


.129 


313.5 


.134 


321.4 


257.9 




1909-14 


.146 


327.0 


.148 


319.8 


256.3 


Carloforte, Sardinia. X= — S* 19' 


1900-05 


.101 


271.3 


.097 


275.9 


267.6 




1900-06 


.094 


283.0 


.101 


278.2 


209.9 




1906-11 


.119 


270.2 


.121 


276.7 


268.4 




1906-12 


.104 


284.6 


.109 


276.3 


268.0 




1910-15 


.107 


273.1 


.122 


268.7 


260.4 




1911-17 


.108 


279.1 


.098 


273.4 


265.1 




1912-17 


.085 


263.9 


.085 


272.1 


263.8 


Qaithersburg, Maryland. X- +77* 12' 


1900-06 


.027 


175.8 


.023 


196.6 


273.8 




1900-06 


.024 


227.0 


.028 


205.4 


282.2 




190a-08 


.024 


266.9 


.013 


255.8 


333.0 




1903-09 


.023 


152.7 


.012 


179.6 


256.7 




1906-11 


.058 


168.7 


.050 


185.8 


263.0 




1906-12 


.049 


203.2 


.052 


182.6 


250.8 




1908-14 


.057 


146.7 


.050 


195.3 


272.5 




1909-14 


.070 


176.3 


.073 


163.1 


240.3 
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Table 3. — Harmonic constants for the variation of latitude — Continued. 

ANNNtJAL COMPONBNT-Contlnued. 





Years of 
series 
inclu- 
sive. 


Before clearing. 


After clearing. 




Station and longitude. 


AmpU- 

tude. 

R 


Epoch. 


Ampli- 
tude. 
R 


Epoch. 


r+x 


Cincinnati, Ohio. X— +84*25' 


1900-05 
1900-06 
1903-08 
190^-11 
1906-12 
1908-13 
1910-15 

1900-05 
1900^ 
1906-11 
1906-12 
1910-15 
1911-17 
1912-17 


It 

.014 
.022 
.041 
.052 
.061 
.065 
.086 

.048 
.024 
.069 
.031 
.067 
.057 
.075 


e 

159.8 
246.9 
283.7 
187.6 
216.9 
182.9 
168.8 

84.6 
88.4 
68.4 
67.1 
96.0 
117.6 
90.6 


n 

.009 
.022 
.029 
.050 
.056 
.053 
.099 

.037 
.035 
.055 
.053 
.082 
.057 
.064 


e 

219.4 
217.2 
261.2 
202.5 
200.8 
178.3 
165.6 

85.7 
92.8 
76.0 
67.8 
96.9 
101.6 
100.5 


e 

303.8 


Ukiah, GaliCornia. X— fl23'13'„ 


301.6 
345.6 
286.9 
285.2 
262.7 
250.0 

208.9 




216.0 
199.2 
191.0 
220.1 
224.8 
223.7 



14-MONTH COMPONENT. 



Mizusawa, Japan. X»— 141*08'. 



Tschardjui, Turkestan. X— -63*29'. 



Cailoforte, Sardinia. X— 8*19'. 



Oaithersburg, l£aryland. X— 1-77*12'. 



Cincinnati, Ohio. X-+84*25'. 



Ukiah, California. X- +123* 13 





ti 


. 


n 


e 


1900-05 


0.142 


179.3 


0.134 


177.5 


1900-06 


.123 


176.4 


.134 


176.5 


1906-11 


.218 


165.8 


.213 


163.2 


1906-12 


.215 


163.6 


.226 


165.3 


1910-15 


.201 


174.8 


.213 


175.6 


1911-17 


.177 


175.8 


.175 


179.4 


1912-17 


.172 


187.8 


.168 


185.2 


1900-05 


.135 


106.8 


.132 


103.0 


1900-06 


.115 


98.9 


.123 


102.3 


1906-11 


.209 


100.7 


.208 


97.3 


1906-12 


.211 


95.0 


.218 


97.8 


1908-14 


.229 


109.1 


.234 


106.2 


1909-14 


.241 


103.7 


.242 


107.0 


1900-05 


.130 


42.9 


.128 


39.4 


1900-06 


.125 


33.4 


.130 


36.3 


1906-11 


.204 


33.8 


.208 


31.2 


1906-12 


.204 


32.6 


.207 


35.0 


1910-15 


.216 


50.5 


.224 


52.2 


1911-17 


.197 


53.8 


.192 


55.7 


1912-17 


.184 


62.3 


.185 


60.2 


1900-05 


.144 


322.1 


.143 


321.3 


1900-06 


.142 


317.8 


.144 


318.6 


1903-08 


.155 


310.5 


.155 


310.6 


1903-00 


.174 


311.7 


.174 


311.4 


1906-11 


.225 


305.0 


.226 


303.7 


1906-12 


.221 


305.2 


.224 


306.4 


1908-14 


.210 


317.0 


.212 


315.4 


1909-14 


.213 


317.6 


.215 


319.4 


1900-05 


.156 


317.3 


.157 


317.1 


1900-06 


.153 


310.5 


.154 


311.3 


1903-08 


.150 


299.2 


.156 


299.3 


1906-11 


.207 


297.5 


.210 


296.4 


1906-12 


.206 


296.2 


.207 


297.8 


1908-13 


.198 


303.1 


.194 


303.6 


1910-15 


.207 


323.2 


.216 


323.9 


1900-05 


.135 


277.5 


.131 


275.7 


1900-06 


.131 


273.4 


.134 


273.0 


1906-11 


.225 


256.3 


.222 


255.5 


1906-12 


.216 


258.6 


.221 


258.4 


1910-15 


.199 


279.7 


.206 


279.7 


1911-17 


.207 


283.2 


.208 


284.5 


1912-17 


.220 


287.5 


.217 


286.2 



36.4 
35.4 
22.1 
24.2 
34.5 
38.3 
44.1 

39.6 
38.8 
33.8 
34.3 
42.7 
43.5 

31.1 
28.0 
22.9 
26.7 
43.9 
47.4 
51.9 

38.5 
35.8 
27.8 
28.6 
20.9 
23.6 
32.6 
36.6 

41.5 
35.7 
23.7 
20.8 
22.2 
28.0 
48.3 

38.9 
36.2 
18.7 
21.6 
42.9 
47.7 
49.4 
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Much might be said in interpretation of the results in Table 3, but 
the discussion may conveniently be postponed (see pp. 51 to 64) 
except for a few remarks bearing directly on the problem of detecting 
a progressive variation of latitude. 

It is at once evident that the so-called ''harmonic constants" are 
only approximately constant. Their precise values have a significance 
only as means for the period of time from which thev were obtained. 
However, if we limit ourselves to comparing merely the harmonic 
constants of the various stations derivea from observations covering 
one and the same period of time for all stations, there are still some 
useful conclusions to be drawn. 

The values of f -f-X for the annual component show such marked 
variations that it is evident that something besides the motion of the 
pole enters into them. This something is generally called the Kimura 
term; it has been variously explained and is doubtless due to a 
combination of causes. It is generally assumed in most discussions 
of the variation of latitude that the Kimura is the same at all stations 
on the same parallel at any given time. This approximate equality 
comes out from the discussions later on (see p. 54) but there is no 
reason except mathematical convenience for assuming that the 
equahty is absolute.* This shows the advantage of dealiug with each 
station by itself, quite apart from any motion of the pole. The 
apparent annual variation of latitude at any given station depends 
not only on the motion of the pole but on the individual Kimura 
term oi that station. From the observations at one station alone 
we can not separate the two, but for the present purpose, namely, 
that of eliminating the annual portion of the variation at that station, 
we do not need to know how much is due to the motion of the pole 
and how much to the individual peculiarities of the Kimura term. 

For the 14-month component, on the other hand, the values of 
R and f +X show a tolerable constiancy from station to station, which 
is favorable to the hypothesis that the variation of latitude iu the 
14-month period is due chiefly to the motion of the pole in an orbit, 
approximately circular, about its mean position. Any individual 
peculiarities of the station as respects the variation oi latitude in 
the 14-month period seem to be quite small. (See p. 61.) 

SECTION 3. MEAN LATITUDES BY CALENDAR YEARS AND BY 14-MONTH 

PERIODS. 

Having obtained the harmonic constants, we could do just what 
was done in Dyson's paper; we could adopt a set of harmonic con- 
stants for a given period and a given station, compute the correspond- 
ing annual portion of the variation of latitude, and by subtracting 
this annual portion from the ordinates of the curve mentioned on 
page 8, which may be called the curve of the complete variation, 
we should get the ordinates of a new curve of exactly the same nature 
as Dyson's curves (p. 3). The new curve would presumably 
contam only the variation of latitude in the 14-month period (this 
statement being subject to the same quaUfications as on p. 3) and 

> The weight of authority seems to i&vot the opinion that the Kimura term is chiefly due to the seasonal 
variation of the refraction, and to changes in the refraction in passing from the observatcny to the outeide 
air rattier more than to refractions due to seasonal barometric gradients. If this be true, it ma^ be con- 
sidered as certain that, owing to differences of climate and local topography, the Kimura term vml not be 
absolut^y alike at all stations. 
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by taking the means of its ordinates over a 14-month period the 
periodic part would disappear and the means would bring to light 
the progressive variation or latitude, if any. Furthermore, the process 
may be reversed; that is, the 14-month portion of the variation may 
be computed from the assumed harmonic constants and subtracted 
from the curve of the complete variation to obtain the ordinates of a 
new curve, which will, presumably, contain only the annual portion 
of the variation. By taking means over a year the annual periodic 
part of the variation should disappear and the yearly means should 
put in evidence the progressive variation of latitude, if any. 

The only difference between the process just outlined and what was 
actually done is that the drawing and reading of a second actual curve 
were dispensed with. The ordinates of the original curve of the com- 
plete variation had already been obtained at suitable intervals in 
the process of making the harmonic analysis. The portion of the 
variation, annual or 14-monthly as the case may be, that is to be 
removed may be computed from the adopted harmonic constants and 
subtracted from the ordinates of the curve of complete variation, thus 
obtaining the ordinates of the new curve and the means of these 
new ordmates may be taken over the proper period, all without 
actually drawing the new curve. Furtnermore, the computation 
may be simphfied, for instead of applying our corrections individually 
to the curve of complete variation and taking the mean of the cor- 
rected result, we may simply apply the mean of the corrections to 
the mean of the imcorrected ordinates, both means covering the proper 
period. The calculation of the mean correction is rendered very 
simple by a formula of the type of (1) on page 71. 

The most obvious weakness of any calculation of this sort is in the 
choice of the harmonic constants for a particular year or period of 
14 months. Since the so-called constants are valid only as mean 
values over a given period of time and it is not easy to determine 
mathematically their values at any given instant, there is necessarily 
an arbitrary element in the assignment of values for any given year 
or 14-montn period. The consequences of this imcertainty are not, 
however, so serious as might be supposed. The calculations will be 
given for this part of the mvestigation in more detail than for most 
of the other parts, and it will not be difficult to satisfy oneself that 
any reasonable choice of the harmonic constants would lead to 
substantially the same result for the progressive change of latitude, 
which is the quantity sought. 

The more important stages of computations for the case where the 
mean latitude is taken by calendar years are given on pages 15 and 16. 
The first column under each station gives the calendar year over 
which the mean is taken and the second the mean itself. In this 
mean 13 values were used, already read off from the curve for the 
purposes of the harmonic analysis. The first and last of the 13 for 
the year 1900, for example, are for the dates 1900.00 and 1901.00. 
These two end values are given only half weight. The mean lati- 
tude evidently applies exactly to the middle of the year 1900. The 
third and fourth columns contain the amplitudes and epochs (R and 
f) of the 14-month component as adopted for the purpose of clearing 
the annual mean. The epochs are reduced to the beginning of 1900 
so as to be comparable with those of the preceding table. The way 
in which these quantities were determined will be explained in the 
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next poraCTaph. Given the values of B and f , the quantity to be 
added to me uncorrected mean to clear it of the outstanding effects 
of tho 14-month component is found to be 



0.174Bcos(332?0-f-55?98r^) 



(1) 



(See p. 73.) In this expression n is the number of years after 1900; 
e. g., for the year 1907, n « 7. The values of the above expression are 
given in the fifth column; the sixtJbi column, which is the sum of the 
second and fifth columns, is the cleared mean which applies to the 
middle of the appropriate calendar year; the trend of tnese should 
bring to light the progressive variation of latitude, if any. 

As has been said, the 14-month component of the variation of 
latitude appears to be due almost entirelv to the motion of the pole, 
but for the present purposes it is enougn to assume that the varia- 
tions with the time in R and f at any station correspond to variations 
in the motion of the pole. In an article by Wanach " there are given 
the instantaneous values of the amplitude and phase of the polar 
motion in the 14-month period. Wanach's quantities are computed 
on the assumption that tne annual component repeats itself exactly, 
an assumption known to be erroneous and one which is therefore quite 
sufl&cient to explain the decidedly erratic variations in the quantities 
tabulated by Wanach. These he calls B, the amplitude of the polar 
motion and P its initial phase reckoned from a point midway between 
maximum and minimum. Although the values of R and P given 
for every tenth of a year are so far from constant, their means over 
a calendar year are more stable and these means *® are given in the 
following table. In taking means the end values were given half 
weight and the year 1917, which is incomplete, is treated as if com- 
plete. 

Table 4. — Mean amplitude and phase of the 14''fnjO'nMy polar motion. 



Calendar year. 



1900 
1901 
1902 
1903 
1904 
1905 
1906 
1907 
1908 



Ampli- 
tude, J2. 


Phase, P. 


// 


e 


0.133 


63 


.166 


55 


.120 


53 


.144 


52 


.139 


52 


.089 


53 


.138 


80 


.189 


94 


.258 


78 



Calendar year. 



Phase, P. 



1909 
1910 
1911 
1912 
1913 
1914 
1915 
1916 
1917 




The values of 22 in Table 4 are directly comparable with those in 
Table 3 and if the means are taken over corresponding periods the 
two sets of values are nearly equal. The quantity IP should be 
nearly equal to 90°— (f +X) for the 14-component and this is seen to 
be very nearly true. An example will explain how the values of R 
and f used on pages 15 to 17 for clearing the annual mean were ob- 
tained from Tallies 3 and 4. 

Let us take the year 1902 at Mizusawa and compare it with the 
harmonic analysis at the same station for the years 1900-05. The 

• Die Chandlersche und die Newoomhsche Periode der Polhewegung. Puh. No. 34 (new series), Zentral- 
bureau der Intemationalen Erdmessung. Berlin, 1919, p. 23. 

10 The ordinary mathematical mean was taken. This process is not quite rigorous in such a case and 
tends to exaggerate a little the vahMB of fi, bat ft>r preaeot purposes it is tuffldwttly accurate^ 
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mean B for the years 1900-05 from Table 4 is 0?132. The analysis 
(Table 3) rives 0ri34. In Table 4 the B for 1902 is 0?120 or 
0.^012 less tnan the mean^ so it is assumed that the B for Mizusawa 
in 1902 should also be 07012 less than 0.^134 or 0n22. In a similar 
way the mean P for tiie period is 55° and for 1902 is 53° or 2° less. 
Since f and P vary in opposite directions from their mean values the 
f for Mizusawa in 1902 should be 2° greater than the f in Table 3 
corresponding to the whole period; i. e., its value is 177?^5 + 2° = 
179?5, or 180 by the ordinary convention for rounding oflF to whole 
degrees. The period 1900-06 and the corresponding harmonic analy- 
sis might be used in a similar way and values for 1902 deduced. The 
results are practically identical with those given. The values, 22= 
Of 122 and f = 180°, are found in the third and fourth colimins of the 
table below. Some years, such as 1911, are included in the periods 
covered by several harmonic analyses; in such cases the mean of 
the several comparisons, which generally differ but little from one 
another, is used. 

Table 5. — Mean latitude by calendar years. 



MIZUSAWA 


., JAPAN 


• 






• 


Year. 


Second3 
of mean 
latitude, 

un- 
cleared. 


Harmonic con- 
stants, 14-mQnth 
component. 


Clear- 
ance. 


Seconds 
of mean 
latitude, 




R 


r 


cleared. 


1900 


n 

3.634 
.627 
.559 
.562 
.590 

3.591 
.613 
.609 
.589 
.546 

3.567 
.559 
.628 
.627 
.612 

3.554 

.571 
.648 


n 

a 134 

.168 
.122 
.146 
.140 

0.090 
.133 
.180 
.248 
.258 

0.222 

■■• .244 

•••.257 

.189 

.174 

0.161 
.124 
.129 




170 
178 
180 
182 
182 

180 
155 
142 
158 
170 

170 
172 
178 
184 
176 

189 
184 
198 


II 

-a 022 

- .004 
+ .016 
+ .024 
+ .007 

-a 010 

- .022 

- .029 
-f .003 
+ .043 

+ao30 

- .004 

- .041 

- .025 

- .001 

+0.026 
+ .019 

- .003 


i 1 
3.612 


1901 


3.623 


1902 


3.575 


1903 


3.586 


1904 


3.597 


1905 


3.581 


1906 


3.591 


1907 


3.580 


1908 


3.592 


1909 


3.589 


1910 


3.597 


1911 


3.555 


1912 


3.587 


1913 


3.602 


1914 


3.611 


1915 


3.580 


1916 


3.590 


1917 


3.645 







TSCHARDJXJI, TURKESTAN. 



1900 
1901 
1902 
1903 
1904 

1905 
1906 
1907 
1908 
1909 

1910 
1911 
1912 
1913 
1914 



10.698 


0.128 


96 


-0.012 


.716 


.161 


104 


-.028 


.643 


.115 


106 


- .008 


.663 


.140 


107 


+ .013 


.636 


.134 


107 


+ .023 


10.602 


0.084 


106 


+0.009 


.536 


.124 


86 


.000 


.617 


.173 


76 


- .022 


.714 


.242 


89 


- .038 


.684 


.256 


100 


+ .002 


ia672 


0.223 


404 


+a034 


.724 


.250 


102 


+ .039 


.737 


.270 


117 


- .006 


.825 


.208 


124 


- .084 


.886 


.196 


116 


- .030 



10.686 
ia688 
10.635 
10.676 
10.659 

10.611 
10.536 
10.595 
10.676 
10.686 

10.706 
10.768 
10.731 
ia791 
ia856 



89016°— 22- 
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Table 5. — Mean latitude by calendar years — Continued. 

CARLOFORTE, SARDINIA. 



• 

Year. 


Seconds 

of mean 

latitude, 

un- 
cleared. 


Harmonic con- 
stants, 14-month 
component. 




R 


r 


1900 


8.909 
9.010 
8.956 
8.935 

8.946 

8.948 
.911 
.904 
.954 
.973 

8.960 

.922 

.932 

' 9.004 

9.027 

9.020 
.027 
.040 


0.130 
.162 
.116 
.142 
.136 

0.086 
.122 
.167 
.236 
.246 

0.218 
.245 
.262 
.201 
.189 

0.176 
.141 
.145 


e 

31 


1901 


39 


1902 


41 


1903 


42 


1904 


42 


1905 


41 


1906 


23 


1907 


11 


1908 


27 


1908 


39 


1910 


41 


1911 


38 


1912 ^ 

1913 


53 
60 


1914 


52 


• 
1915 


66 


1916 


59 


1917 


74 







Clear- 
ance. 



+0. 012 

- .015 

- .020 

- .013 
+ .010 

-fa015 
-I- .019 
+ .009 

- .033 

- .037 

-0.001 
+ .034 
+ .038 

- .004 

- .027 

-0.025 

- .003 
+ .023 



Seconds 

of mean 

latitude, 

<deared. 



8.921 
8. 995 
8.936 
8.922 
8.956 

8.963 
8.930 
8.9J13 
8.921 
8.936 

8.950 
8.956 
8.970 
9.000 
9.000 

8.995 
9.024 
9.063 



GAITHERSBURG, MD. 



1900 

1901 

1902 

1908 

1904 

1906 

1906 

1907 

1908 

1909 

1910 

1911 

1912 4 

1918 

1914.. 



13. 115 
.178 
.238 
.271 
.217 

13.215 
.222 
.208 
.232 

.278 

13.297 
.307 
.288 
.272 
.327 



0.144 
.176 
.130 
.1^ 
.143 

0.093 
.137 
.186 
.249 
.254 

a 217 
.245 
.256 
.184 
.172 



314 
322 
324 
326 
326 

324 
297 
284 
298 
309 

312 
310 
327 
334 
326 



+0.024 
+ .021 

- .005 

- .025 

- .020 

+0.001 
+ .012 
+ .031 
+ .025 

- .023 

-0.038 

- .026 
+ .026 
+ .032 
+ .016 



13.139 
13.199 
13.233 
13.246 
13.197 

13. 216 
13.234 
13.239 
13.257 
13.255 

13.259 
13.281 
13. 314 
13.304 
13.343 



CINCINNATI, OfflO. 



1900 
1901 
1902 
1903 
1904 

1905 
1906 
1907 
1906 
1909 

1910 
1911 
1912 
1913 
1914 

1915 



19.318 
.274 

. Sari 

.436 
.443 

19.292 
.326 
.277 
.341 
.368 

19.408 
.442 
.410 
.392 
.473 

19.548 



a 156 
.189 
.143 
.159 
.153 

0.103 
.132 
.174 
.232 
.235 

0.206 
.233 
.246 
.178 
.188 

a 175 



308 
316 
318 
318 
318 

317 
289 
274 
289 
302 

306 
304 
320 
329 
326 

339 



+0.025 
+ .025 

- .003 

- .025 

- .023 

-a 002 
+ .009 
+ .027 
+ .029 

- .017 

-0.036 

- .027 
+ .021 
+ .031 
+ .017 

-a 018 



19.343 
19.299 
19.389 
19. 411 
19.420 

19.290 
19.335 
19.304 
19.370 
19.351 

19.372 
19.415 
19.431 
19.423 
19.490 

19.530 



STATIONS OF THE INTERNATIONAL LATITUDE SERVICE. 
Table 5. — Mean latitude by calendar yean — Oontiiiued. 

UEIAH, CALIF. 



17 



Year. 


Seconds 
of mean 
latitude, 

un- 
deared. 


Harmonic con- 
stants, 14-month 
component. 


Clear- 
ance. 


Seconds 
of mean 
latitude, 




R 


f 


cleared. 

• 


1900 


It 

12.078 
.024 
.048 
.115 
.142 

12.136 
.118 
.079 
.028 
.103 

12.163 
.205 
.185 
.127 
.182 

12.251 
.258 
.264 


11 

ai33 

.166 
.120 
.146 
.140 

0.090 
.133 
.182 
.250 
.260 

0.221 
.250 
.274 
.212 
.200 

ai87 
.165 
.109 


• 

268 
276 
278 
279 
279 

278 
250 
285 
251 
263 

267 
264 
282 
288 
280 

203 
286 
300 


+aoio 

+ .029 
+ .011 

- .011 
-.024 

-a Oil 

- .006 
+ .013 
+ .043 
+ .012 

-0.027 

- .043 

- .007 
+ .090 
+ .034 

+aoo5 

- .018 

- .028 


tt 
12.068 


1901 


12. 058 


1902 


12.0S9 


1908 


12.104 


1904 


12.118 


1905 


12.126 


1906 


12.112 


1907 


12.008 


1908 


12.071 


1900 


12.115 


1910 


12.136 


1911 


1Z162 


1912 


1Z178 


1913 


12.157 


1914 J 


12.216 


1915 


12.256 


1916 


12.240 


1917 


12.236 







A process the reverse of that set forth m Table 5 consists in taking 
the mean over a period of 14 months (more accurately 432.5 days) 
and clearing this mean of the effects of the annual component. In 
Chapter I me mean was taken over a period of exactly 14 months 
(i4 years) and the difference between this, and the more accurate 
value of 432.5 days, was shown to have a nearly negligible effect. 
The same thing could be done here, but the curve reading used 
would be the same readings as those forming the basis of Table 5. 
We get, therefore, a determination more nearly independent of the 
preceding if we utUize the curve reading taken independently for 
the harmonic analysis of the 14-month (i. e., 432.5-day) component. 
The ^ef disadvantage is that, owing to the practice of banning 
the series to be analyzed at the banning of a year, the 432.5-day 
periods covered do not run consecutively and without overlap. 

Table 6 gives the results of the taking and clearing of means. Its 
arrangement differs only slightly from that of Table 5. The first 
column will be explained later. The second column, which corre- 
sponds to the first column in Table 5, gives the dates in years and 
aecimals of the beginning and ending of the period covered. The 
third column is the mean of 13 reacungs covering this period, the 
initial and final readings having onl^ half weight. The mean then 
corresponds to the midale of the period. The fourth and fifth col- 
umns contain the harmonic constants of the annual component 
adopted for clearing the mean. For lack of a b etter principle oi choice 
the constants in these columns are simply the mean of all available 
determinations in Table 3 that cover the period in question. The sixth 
column contains the clearance of the mean. It is given by the 
formula 

Clearwice= 0.142 R cos (213M -f + 360°/). (2) 
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The R and f of formula (2) are those in columns four and five ; the 
quantity/ is the fraction of a year elapsed at the banning of the period 
to be cleared; thus, for the first line of lifizusawa /= 0.000, for the 
second line /= 0.184, for the third line /= 0.368, etc. The seventh 
column representing the cleared mean is the sum of the third and 
sixth columns. 

Table 6. — Mean latitiuks by periods of4S2.5 days, 
MIZUSAWA, JAPAN. 



Oroup No. 


Period included. 


Seconds 
of mean 
latitude, 

un- 
cleared. 


Harmonic con- 
stants, annual 
component. 


Clear- 
ance. 


Secx>iids 

of in4«n 

latitude, 

cleared. 




J2. 


r. 


1 


1900.000-01.184 
1901.184-02.368 
1902. 36a-03. 552 
1903.552-04.736 
1904.736-05.920 

1905.920-07.104 
1906.000-07.184 
1907. 184r-08. 368 
1908.368-09.552 
1909.552-10.736 

1910.000-11.184 
1910. 73&-11. 920 
1911.000-12.184 
1911.184-12.368 
1911.920-13.104 

1912.000-13.184 
1912. 184-13. 368 
1912. 368-13. 552 
1913. 184-14. 368 
1913.368-14.552 

19ia5S2-14.736 
1914.368-15.552 
1914. 652-15. 736 
1914.736-15.920 
1915. 552-16. 736 
1915.736-16.920 

1916.736-17.920 
1916.920-18.104 


03.633 
.603 
.569 
.585 
.585 

03.604 
.597 
.569 
.580 
.586 

03.612 
.576 
.572 
.557 
.607 

03.612 
.602 
.572 
.615 
.590 

08.602 
.685 
.590 
.602 

.607 
.629 

03,650" 
.665 


0.106 
.106 
.106 
.106 
.106 

0.117 
.121 
.128 
.128 
.132 

0.126 
.126 
.121 
.121 
.121 

0.118 
.118 
.118 
.117 
.117 

0.117 
.117 
.117 
.117 
.117 
.117 

a 106 

.106 


e 

16 
16 
16 
16 
16 

13 
12 
10 
10 
12 

13 
13 
11 
11 
11 

12 
12 
12 
13 
13 

13 
13 
13 
13 
13 
13 

10 
10 


-a 014 

- .002 
+ .013 
+ .012 

- .003 

-a 016 

- .016 
.000 

-1- .016 
+ .014 

-0.017 

- .005 

- .016 

- .001 

- .017 

-a 016 

- .001 
+ .015 

- .001 

- .015 

-i-aoi3 

-1- .015 
+ .013 

- .004 
+ .013 

- .004 

-a 004 

- .016 


03.619 


2 


.601 


3 


.582 


4 


.597 


5 


.582 


6 


03.588 


6 


.581 


7 


.560 


8 


.506 


9 


.600 


9 


03.595 


10 


.571 


10 


.556 


10 


.566 


11 


.590 


11 


03.506 


11...: 


.601 


11 


.587 


12 


.614 


12 


.575 


12 


03.615 


13 


.600 


13 


.603 


13 


.508 


14 


620 


14 


.625 


15 


03.646 


15 


.650 







TSCHARDJUI, TURKESTAN. 



1 
2 
3 
4 
5 

6 
6 
7 
8 
8 

9 

9 

9, 

10 

10. 

10 
11. 
11. 
11. 
12. 

12. 
13. 
13. 



1900.000-01.184 


ia694 


a 103 


338 


-a 008 


1901.184-02.368 


.676 


.108 


338 


+ .008 


1902.368-03.552 


.632 


.103 


338 


- .014 


1903.552-04.736 


.662 


.103 


338 


+ .004 


1904.736-05.920 


.688 


.108 


S38 


- .011 


1905.920-07.104 


ia668 


a 118 


330 


-a 014 


1906. 000-07. 184 


.543 


.123 


327 


- .007 


1907. 184-08. 368 


.608 


.133 


322 


+ .014 


1908. t)0(M)9. 184 


.683 


.136 


322 


^ .006 


1908. 368-09. 552 


.662 


.136 


322 


+ .018 


1900. OOO-ia 184 


ia692 


a 136 


322 


-a 006 


1909. 184-la 368 


.682 


.136 


322 


+ .014 


1909. 552-ia 736 


.702 


.136 


822 


.000 


19ia 184-11. 368 


.721 


.136 


322 


+ .014 


19ia368-lL652 


.720 


.136 


322 


-f .018 


1910-736.11.920 


ia759 


0.136 


322 


-a 018 


1911.368-12.562 


.731 


.136 


322 


+ .018 


1911. 552-12. 736 


.760 


.136 


322 


.000 


1911. 920-13. 104 


.756 


.136 


322 


- .014 


1912.552-13.736 


.741 


.137 


321 


.000 


1912.736-13.920 


ia780 


a 137 


321 


-a 018 


1913.736-14.920 


.875 


.141 


321 


- .018 


1913.920-15.104 


.860 


.141 


321 


- .015 



ia686 
.684 
.618 
.666 
.627 

ia544 
.536 
.622 
.677 
.680 

ia686 



.702 
.735 
.738 

ia741 
.749 
.760 
.7^ 
.741 

ia702 
.857 
.854 



STATIONS OF THE INTEftlTATlONAL LATITUDE SERVICE. 
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Table 6. — Mean latitvdes by periods of 4St.5 days — Continued. 

CARLOFORTE, SARDINIA. 



Oroiip No. 


Period included. 


Seconds 
of mean 
latitude, 

un- 
cleared. 


Harmonic con- 
stants, annual 
component. 


Clear- 
ance. 


Seconds 
of mean 
latitude, 




R. 


r. 


cleared. 


1 


1900 000-01. 184 
1901.184-02.368 
1902.368-03.552 
1903.552-04.736 
1904.736-05.920 

1905.920-07.104 
1906.000-07.184 
1907. 184-08. 368 
1908.368-09.552 
1900.552-10736 

1910 OOO-ll. 184 
1910730-11.920 
1911. 000-12. 184 
1911.184-12.368 
1911.920-13.104 

1912. 000-13. 184 
1912. 184-13. 368 
1912. 388-13. 552 
1913. 184-14. 368 
1913. 368-14. 552 

1913.552-14.736 
1914.368-15.552 
1914.552-15.736 
1914. 736-15. 920 
1915.552-16.736 

1915. 736-16. 920 
1916.736-17.920 
1916.920-18.104 


8.921 
.963 
.923 
.962 
.975 

8.927 
.925 
.907 
.903 
.962 

8.952 
.970 
.954 
.953 
.965 

8.957 
.959 
.980 
.987 
.996 

9.005 
8.984 
&997 
9.000 
.026 

9.031 
.073 
.053 


n 

0099 
. 099 
.099 
. 099 
. 099 

O107 
.110 
.115 
.115 
.117 

0112 
.112 
.107 
.107 
.107 

O103 
.103 
.103 
.102 
.102 

O102 
.102 
.102 
.102 
.102 

O102 
.092 
.092 




277 
277 
277 
277 
277 

277 
277 
276 
276 
274 

274 
274 
274 
274 
274 

272 
272 
272 
271 
271 

271 
271 
271 
271 
271 

271 
273 
273 


II 

-1-0 006 
H- .014 
-f .005 

- .010 

- .013 

-OOOl 
+ .007 
+ .016 
4- .006 

- .012 

+O008 

- .014 
+ .007 
+ .015 

.000 

+O008 
+ .014 
+ .004 
+ .014 
+ .004 

-OOU 
+ .004 

- .011 

- .013 

- .011 

-0 013 

- .012 
.000 


n 

08.927 


2 


.977 


3 


.928 


4 


.952 


5 


.962 


6 


08.926 


6 


.932 


7 


.923 


8 


.909 





.950 


9 


08.960 


10 


.956 


10 


.961 


10 


.968 


11 


.965 


11 


08.965 


11 


.973 


11 


.984 


12 


09.001 


12 


.000 


12 


0&994 




.988 


13 


.986 




.987 


14 


09.015 


14 


09.018 




.061 


15 


.068 







QAITHERSBUFvG, MD. 



1 


1900.000-01.184 
1901.184-02.368 
1902.368-03.552 
1903.000-04.184 
1903.552-04.736 

1904.184-05.368 
1905.368-06.552 
1905.920-07.104 
1906.000-07.184 
1906.552-07.736 

1907.184-08.368 
1907.736-08.920 
1908.000-09.184 
1908.368-09.552 
1908. 920-10 104 

1909.000-10184 
1909.184-10368 
1909.552-10736 
1910 184-11. 368 
1911.368-12,552 

1911.562-12.736 
1911.920-13.104 
1912.552-13.736 
1912.730-13.920 
1913.736-14.920 
1913.920-15.104 


13. 143 
.201 
.254 
.244 
.214 

13.192 
.225 
.229 
.232 
,233 

13.260 
.252 
.245 
.265 
.252 

13.253 
.257 
.263 
.259 
.307 

13.307 
.302 
.329 
.309 
.325 
.331 


O026 
.026 
.019 
.019 
.019 

O019 
.031 
.031 
.033 
.033 

O039 
.039 
.044 
.044 
.044 

O051 
.051 
.051 
.060 
.060 

0.060 
.060 
.061 
.061 
.066 
.066 


201 
201 
209 
209 
209 

209 
201 
201 
202 
202 

200 
200 
194 
194 
194 

181 
181 
181 
182 
182 

182 
182 
180 
180 
179 
179 


+O004 
+ .001 

- .002 
+ .003 

- .003 

+O001 

- .004 
+ .004 
+ .005 

- .004 

+O001 

.000 

+ .006 

- .006 
+ .006 

+O006 

- .001 

- .005 

- .001 

- .008 

-O005 
+ .009 

- .005 
+ .004 
+ .005 
+ .009 


13. 147 


2 


.202 


8 


.252 


4 


.247 


4 


.211 


5 


13.193 


6 


.221 


6 


.233 


6 


.237 


7 


.229 


7 


13.261 


8 


.252 


8 


.251 


8. 


.259 


9 


.258 


9 


13.250 


9 


.256 


9 


.258 


10 


.258 


11 


.299 


11 


13.302 


12 


.311 


12 


.324 


12 


.313 


13 


.330 


13 


.340 
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Table 6. — Mean latituda by periods of4St.5 (fG^»— Oontiiiued. 

CINCINNATI. OHIO. 




Group No. 


Period bKduded. 


Seconds 
of mean 
latitude, 

un- 
cleared. 


Harmonic con- 
stants, annual 
component. 


Clear- 
ance. 


Seconds 

of mean 

latitude, 




R. 


r. 


cleared. 


1 


lOOa 000-01.184 
1901.184-02.368 
1902.368^)8.552 
1908.000-04.184 
1903. 552-OL 736 

1904. 184-05. 368 
1904.736-05.920 
1905.368-06.552 
1905.920-07.104 
1906.000-07.184 

1906.552-07.736 
1907. 184.0A. 368 
1907.736-06.920 
1906.000-09.184 
1908.368-09.552 

1909.184^10.368 
1909.552-10.736 
1910. 000-11. 184 
1910. 368-11. 552 
1910.736-11.920 

1911.184-12.368 
1911.55^12.736 
1911.920-13.104 
1912. 368-13. 552 
1912.736-13.920 

1913. 55^14. 736 
1914.736-15.920 


19.340 
.297 
.325 
.419 
.^ 

19.371 
.316 
.299 
.332 
.329 

19.316 
.323 
.348 
.360 
.301 

19.350 
.352 
.871 
.408 
.406 

19.371 
.430 
.416 
.457 
.424 

19.488 
.513 


aoi6 

.016 
.020 
.020 
.020 

ao2o 

.020 
.038 
.038 
.042 

a042 
.040 
.040 
.049 
.049 

a067 
.067 
.067 
.067 
.067 

a 067 
.067 
.067 
.009 
.069 

a 076 

. Uvtf 


• 

216 
216 
283 
233 
233 

283 
283 
220 
220 
220 

220 
211 
211 
211 
211 

187 
187 
187 
187 
187 

187 
187 
187 
182 
182 

172 
172 


+aoo3 

+ .001 

- .001 
+ .003 

.000 

+aoo2 

- .001 

- .003 
+ .004 
+ .006 

-a 006 

+ .003 

.000 

+ .007 

- .005 

-a 001 

- .007 
+ .009 

- .009 

- .003 

-a 001 

- .007 
+ .010 

- .009 

- .004 

-0.005 
+ .008 


19L343 


2 


.298 


8 


.324 


4 


.422 


4 


.423 


5 


Id. 373 


6 


.315 


6 


.296 


6 


.336 


6 


.335 


7 


19.310 


7 


.325 


8 


.348 


8 


.367 


8 


.296 


9 


19l348 


9 


.345 


10.: 


.380 


10 


.380 


10 


.403 




10l37D 


11 


.423 




.426 


12 


.448 


12 


.420 




19.483 


14 


.621 







UKIAH, CALIF. 



1 
2. 
3 
4. 
5 

6 
6 
7 
8 
9 

9 
10 
10 
10 
11 

11 
11 
11 
12 
12 

12 
13 
13 
13 
14 

14 
15 
16 



1900.000-01.184 
1901. 184-02. 368 
1902.368-03.552 
1903.56^-04.736 
1904.736-05.920 

1905.920-07.104 
1906.000-07.184 
1907.184-06.368 
1908. 368-09. 552 
1909. 552-10. 736 

1910.000-11.184 
1910.735-11.920 
1911.000-12.184 
1911. 184-12. 368 
1911. 920-13. 104 

1912.000-13.184 
1912.184^13. 368 
1912.368-13.552 
1913. 184-14. 368 
1913.368-14.652 

1913.55^14.736 
1914.368-15.652 
1914.652-16.736 
1914.736-15.920 
1916.552-16.736 

1915.736-16.920 
1916.736-17.920 
1916. 920-1& 104 



1Z089 


0.036 


.050 


.036 


.085 


.036 


.115 


.036 


.120 


.036 


12.116 


0.046 


.116 


.048 


.002 


.054 


.104 


.054 


.127 


.063 


12.149 


0.062 


.149 


.062 


.175 


.062 


.176 


.062 


180 


.062 


12.185 


a064 


.186 


.064 


.173 


.064 


.176 


.068 


.182 


.068 


12.195 


0.068 


.257 


.068 


.256 


.068 


.279 


.068 


.239 


.068 


12.229 


a068 


.212 


.060 


.220 


.060 



89 


-a 003 


89 


- .005 


89 


- .001 


89 


+ .004 


89 


+ .004 


81 


-a 002 


79 


- .005 


72 


- .007 


72 


.000 


80 


+ .008 


86 


-a 006 


86 


+ .005 


86 


- .005 


86 


- .009 


89 


- .001 


92 


-a 006 


92 


- .009 


92 


- .003 


100 


- .010 


100 


- .004 


100 


+aoo6 


100 


- .004 


100 


+ .006 


100 


+ .009 


100 


+ .006 


100 


+aoo9 


100 


+ .008 


100 


+ .001 



12.086 
.045 
.064 
.119 
.124 

12.114 
.110 
.085 
.104 
.185 

12.144 
.154 
.170 
.167 
.179 

12.180 
.177 
.170 
.165 
.178 

12.201 
.263 
.262 
.288 
.245 

12.288 
.220 
.221 
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SECTION 4. BATES OP CHANQB OF THE LATITUDES. 

The results in the last column of Table 5 are plotted in figures ia, 
6a, 6a, 7a, 8a, and 9a. An inspection of these suggests that uiere is a 
progressive increase in latitude at the American stations, Gaithers- 
burg, Cincinnati, and Ukiah, and that in the latter part of the period 



FiQ. ia.—pTognsaiva Bbange □[ latitude at Uicusawa— means by calendar years. 



Flo. 46,— ProgrosdTB change ollatltude at Miiusawa— means by W-month pcrioda. 

covered this increase becomes more rapid. Carloforte shows no 
clear indications of either an increas,e or a decrease during most of 
the time, but toward the end of the time there is a fairly w5l-defined 
tendency to increase. Mizusawa shows a slight decrease, if anytliing, 
at the beginning followed by an increase at the end. Tschardjui is 
very irregular with a sharp mcrease toward the end. 
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This changing rate of increase, which in the case of Mizusan-a even 
reverses the sign of the progressive change of latitude, su^ests that 
it may be desirable to modify the mathematical process used in 
Chapter I for determining mathematically these rates of change. 
After some consideration it was decided to assume the following form 
of observation equation: 

A<l> = x+t,y + t^. (3) 



1900 m& jno tsts 

Fio.M.—Piogressirecluiiige of latitude at TacbaidJui—miAiiBbj'cnlenilai Tears. 

In this equation A<^ is the difference, observed latitude mimis some 
fixed arbitrary latitude; the quantity x is the adjusted value of A^ 
for the time 1912.00; y is the progressive rate of change of latitude 
before 1912.00; 2 is the rate after 1912.00; (^ is the time reckoned 
from 1912.00 for dates before that time and la to be taken as zero 
after that time; f, is the time reckoned from 1912.00 for dates after 
that time and is to be taJken as zero before then. These assumptions 
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mean that the endeavor isluade to fit the observed latitudes to two 
straight lines intersecting on the ordinate for 1912.00, the slopes of 
these lines being determined to fit the observations as closely as 
possible. The particular date 1912.00 was fixed on because that is the 
time when the change is made from the definitive results of the Resul- 
tate to the provisional results of the Nachrichten. The observed 



Fia. 5b, — Progressive change ollatltude al TaiAardjuI— meaiu by It-moatb periods. 

latitude depends on the declinations and these depend on the assumed 
proper motions in declination, which are not absolutely known. 

The dropping out of some of the stars after 1912,00, already noted 
at the banning of this chapter^ involves a change in the mean 
error in the proper motions '^ and hence a change in the observed 

lb Uw thMUT orieaat sqUBHes; what Is 
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latitudes. Ttus change should affect all stations alike, barring 
errors of observations and the effect of the weather in causing some 
pairs to be missed at certain stations and other pairs at other statiom. 
An inspection of the plotted observations su^ests also that 1912.0'i 
is not tar from the time when the progressive rate of variation of tht^ 
latitude seems to undergo a change, though if this time w^ere to bt 



Fm. an.— Pragresalve change of latitude at CarlolOrte— means br cakodar years. 



Fio.d*.— ProgresaTechBiigeofiBlltudeatCBrlolmtB— meanaby l*-moDlliperii>j8. 

determined by inspection of the plotted observations alone, it would 
be placed before 1912.00 rather than after it. 

This may also be due to changes in the star program, for although 
all observed stars were used up to 1912.00[ the later program was 
not identical with that used in 1900 and in mtroducing new pairs to 
replace those dropping out on account of precession, the r^uction 
of the new pairs to the original declination system may not have been 
as accurate as was supposed. 
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F10.7O.— Progressive change oIlatitudeatOBitheraburg— means by calendaryears. 

H35 



I IBQS 1910 

Pio. 76.— pnjgreodyB chaoge ollacltiide at GBicheraburg— means by It-mwitli pe 
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On account of the small number of observations in whicli z is 
involved, it is not very well determined. If it is desired to disreganl | 
any |>ossible change m the rate after 1912.00, it is easy to makej 
y = z in the observation equations and the normal equations, thus I 
making all observations contribute to determining a single average 
rate for the ^ole period, as was done in Chapter I. 

Normal equations for each station, based on observation equations 
of the form given by equation (3), were formed and solved for x, y, 
and z. The quantity s was also made equal to y and the resulting 
normal equations solved. The values of y and s, together with the 
probable errors, also the adjusted latitude of each station for 1912.00, 



Fio. Si.— PTOgreaBlve change of lutltude at Cliiclnnatl— means by calendar years. 

which depends directly on the value of x, are given in Table 8. The 
observations at all stations do not cover the same period, some 
stopping as early as the end of 1914. Since there may be some 
advantage in dealing with the same period for all stations, another 
set of solutions was made based on the observations from 1900 to 
1914, inclusive. The results of these solutions are also given in 
Table 8. The discussion of all these results is deferred until some- 
thing has been said about the analogous observations and normal 
equations based on means by 14-month periods. 

The cleared means by 14-month penods in the last column of 
Table 6 correspond to a number of overlapping periods not inde- 
pendent of one another, so that it would not be proper to use each 
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period as the observed quantity in formation of observation equa- 
tions of equal weight. The procedure adopted was to take the means 
\>j groups as indicated in the first column of the table; many of the 
groups contain only a single 14-month period. These groups were 
chosen so that the time interval between the middle of one group 
and the middle of an adjacent group should be as nearly 1.184 years 
(432.5 days) as practicable. Tnere is still some overiapping, but it 
■was not considered worth while to go deeply into the question of 
■weights. The mean was taken for each group, both with respect to 
times and latitudes. The mean latitudes of the groups furnished 
the known quantities for the observation equations, to each of which 



Fta. e£.— Pnigresaiva change oFlalitude at Clnciiiiiati— means by It-month periods. 

was attributed equal weight. In taking the means two very nearly 
coincident perioos, such as are found in group 11 for Mizusawa, 
Carloforte, orUkiah,with dates 1911.920-13.104 and 1912.000-13.184 
are averaged together to form a single mean and this mean is treated 
as of equal weight with the other members of the group in forming 
the final mean lor the entire group. The mean dates oi each group 
and the corresponding mean latitudes are given in Table 7. It wiU 
be seen that Mizusawa, Carloforte, and Ukiah have identical groups. 
The mean latitudes in Table 7 are shown graphically in Figures 46 
to9h. 
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njua WIU 1915 

FlQ. 9b.— Progreaalvs ctwingeodstltudB at UUah— means by calendar years. 



Fig. Vi. — Progressive dunge of latitude at Uklah— meane by ll-moath periods. 
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1*ABLB 7. — Data for obiervatitm equations to determine progressive variation of latitude — 

means by l^-rrwnth periods. 





Date of 

middle 

of group. 


Seconds of mean lati- 
. tude. 

• 


Tschardjui. 


Oaithersburg. 


CincmnatL 


No. of group. 

• 


Mif.u- 
sawa. 


Carlo- 
forte. 


UWah. 


Date of 

middle 

of group. 


Seconds 
of mean 
lati- 
tude. 


Date of 
middle 
of group. 


Seconds 

of mean 

lati- 

tude. 


Date of 
middle 
of group. 


Seconds 
of mean 
lati- 
tude. 


1 


1900.59 
1901.78 
1902.96 
1904.14 
1905.33 

1906.55 
1907.78 
1908.96 
1910.37 
1911.56 

1912.76 
1913.96 
1915. 14 
1916.24 
1917.42 


13.619 
.601 
.582 
.597 
.582 

13.584 
.569 
. 596 
.598 
.561 

13.594 
.601 
.600 
.622 
.618 


08.927 
.977 
.928 
.952 
.962 

08.929 
.923 
.900 
.955 
.962 

08.974 
.998 
.987 

09.016 
.057 


It 

12.086 
.045 
.084 
.119 
.124 

12.112 
.085 
.104 
.140 
.164 

12.176 
.181 
.268 
.242 
.220 


1900.59 
1901.78 
1902.96 
1904.14 
1905.33 

1906.55 
1907.78 
1908.78 
1909.84 
1911.02 

1912.21 
1913.24 
1914.42 


10.686 
.684 
.618 
.666 
.627 

10.540 
.622 
.678 
.695 
.738 

10.750 
.752 
.856 


1900.59 
1901. 78 
1902.96 
1903.87 
1904.78 

1906.26 
1907.46 
1908.63 
1909.82 
1910. 78 

1912.05 
1913.00 
1914.42 


13. 147 
.202 
.252 
.229 
.193 

13.228 
.245 
.254 
.258 
.258 

13.300 
.316 
.335 


1900.59 
1901. 78 
1902.96 
1903.87 
1905.05 

1906.26 
1907.46 
1908.63 
1909.96 
1910.96 

1911.96 
1912.93 
1914. 14 
1915.33 


19.343 


2 


.298 


3 


.324 


4 


.422 


6 


.344 


6 


19.316 


7 


.318 


8 


.337 


9 


.347 


10 


.394 


11 


19.396 


12 


.431 


13 


.483 


14 


.521 


15 



























The observation equations based on Table 7 are, as before, in the 
general form: A«=x + «.y+<^. (3) 

The only difference in meaning of the symbols is due to a difference 
in the epoch at which the rate changes from y to 0, from which i^ and 
<2 are therefore reckoned, and to which x refers. It would have been 
desirable, for the sake of uniformity, to have taken this epoch 1912.00 
jas before. By inadvertence this was not done and different epochs 
Vere used, all slightly earlier than 1912.000. The use of an earlier 
epoch may not be in itself a disadvantage, for it has been noticed 
tnat the cnange of rate may really occur Before 1912.000, and as the 
difference between 1912.000 and the epochs used was not great, the 
original computations were adhered to. The epochs when the rate 
changes, as used in the observation equations and the normal equa- 
tions, were 1911.56 for Mizusawa, Carloforte, and Ukiah; 1911.02 for 
Tschardjui; 1910.78 for Gaithersbiu-g; and 1910.96 for Cincinnati. 
The times are reckoned from these epochs; t^ is zero after the epoch 
stated for the station in question and ^ zero before that time. The 
normal equations were formed as they were when the means were 
taken by calendar years and the corresponding solutions were made, 
namely : (1) y and z unequal, all observations used ; (2) y and z unequal, 
only oDservation for years 1900 to 1914, inclusive, used; (3) y and z 
equal, all observations used; (4) y and z equal, only observations for 

? rears 1900 to 1914, inclusive, used. The results are given in Table 9 
or comparison with the corresponding results in Table 8. The 
value of X does not give immediately the adjusted latitude in 1912.000, 
as that is not the epoch of the change of rate. For convenience in 
comparison, however, the latitude in 1912 is given, by applying a 
correction for the change between the epoch of the adjustment and 
1912.000. 
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Table 8. — Progrufive change of latUvjde, by calendar years^ with results of the adjustment 

and probable errors. 

[The letter A after the name of a station indicates that all observations (see pp. 26 and 81) have been used; 
the letter B that the observations cover only the years 1900 to 1914, inclusive. When y -^ z the value of p 
gives an average annual rate for the period covered and there is no epoch of change of rate. ProbaUe 
errors (indicated by a db sign) are in units of the fourth decimal place of seconds.] 

p AND z UNEQXTAL. 



Station. 


Epoch. 


Latitude 
1912.000 


Annual rate 
of increase 
of latitude 

before epoch, 


Annua Irate 
of increase 
of latitude 

after epoch, 
z. 


Probable 

error of 

a single 

obeervatian, 

e. 


M1ktl<v&wa a 


1912.000 
1912.000 
1912.000 
1912.000 
1912.000 

1912.000 
1912.000 
1912.000 
1912.000 
1912.000 


/ II 

39 06 03.577 
03.575 
10.686 
08.948 
08.945 

39 08 13.284 
19.378 
19.378 
12.147 
12.142 


-0.0Q23i:8 
- .0025dk7 
+ .0042±29 
-1- .0006±11 
+ .0003±12 

-l-0.0090±ll 
+ .0031i:22 
+ .0031±23 
+ .0070±12 
+ .0063±13 


-H>.0060db20 
-1- .0155±40 
-1- .0602±160 
+ .0191 ±26 
+ .0264±60 

-l-0.0218±63 
-1- .0431±84 
+ .0426±130 
-f .0209±29 
-1- .0260±73 


dbllO 


Mimsawa B 


db 94 


Tschardjui A B 


db370 


Carloforte A 


dbl53 


Carlof orte B 


±161 


Gaithersburg AB 


±147 


f;|T<ninTiati A 


±288 


ntncintifttl B 


±300 


TTV^fth A 


±171 


TTkiAh R... 


±168 







p AND 2 EQUAL. 



Station. 



MizusawaA 

Misusawa B 

Tschardjui AB.. 

Carloforte A 

Carloforte B 

Gaithersburg AB 

Cincinnati A 

CincinnatiB 

UkiahA 

UkiahB 



Lat. 1912.000. 



// 



39 08 03. 595 
03.589 
10.736 
08.961 
08.965 



39 08 



13.294 
19.423 
19.408 
12.172 
12.157 



Annual 
rate of 

increase 
of 

latitude, 
y(-z) 



+0.0004±6 
- .0006±7 
+ .0109±26 
+ .0054±9 
+ .0030±11 

-fa0103±9 
+ .0093±19 
+ .0072±20 
+ .0106± 9 
-f .0084±11 



Probable 
error of 
a single 
observa- 
tion, 
e. 



±138 
±114 
±437 
±206 
±186' 

±150 
±350 

±327 
±196 
±179 



Table 9. — Progressive change of latUudCj by 14-monJth periods^ with results of the adjust- 
ment ana probable errors. 



[For explanation see Table 8.] 
y AND z UNEQUAL. 



Station. 


Epoch. 


Latitude. 
1912.000. 


Annual rate 
of increase 
of latitude 

before epoch, 

y- 


Annual rate 
ofincrease 
of latitude 

after epoch, 
z. 


Probable 

error of 

a single 

observationy 

e. 


Mizusawa A 


1911. 56 
1911. 56 
1911.02 
1911. 56 
1911.56 

1910.78 
1910.96 
1910.96 
1911. 56 
1911.56 


• 1 II 

39 08 03.579 
03.580 

10.720 


-0.0027± 7 
- .0027± 8 
+ .0030-1-32 


II 

-f0.0111± 14 
-H .0122± 43 
-h .0502±112 
+ .0174± 23 
+ .0241± 88 

-|-0.0193± 45 
+ .0378± 54 
+ .0389± 79 
+ .0178± 31 
-1- .0168± 72 


± 86 


Mizusawa B 


± 91 


T<w»harH h^l A B , . . . , 


±356 


Carloforie A 


08.951 - .0000±12 
08.953 - . 0001 -t- 16 


±140 


Carlofcxte B 


±187 


Gaithersburg AB 


39 08 13.293 
19.397 
19.398 
12.158 
12.154 


-l-a 0083-1- 14 
+ .0026±20 
-1- .0026±21 
-1- .0076±15 
-f .0072±14 


±150 


Oincinnftti A .. *. , 


±229 


Cincinnati B 


±232 


lTk\^h A 


±185 


UMahl^ 


±153 
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Table 9. — Progressive change of latitude, by l^-month periods, with results of the adjust 

ment and 'probable errors — Continued. 



y AND z EQUAL. 








Station. 

• 


Lat. 1912.000. 


Annual 

rate of 

increase 

of 
latitude, 

ir(-0 


Probable 
error of 
a single 
observa- 
tion, 
e. 


Mizusawa A 


9 t II 

39 08 03.601 
03.585 
ia736 
08.979 
08.962 

39 08 13.395 
19.419 
19.403 
12.175 
12.157 


II 

+0.0013± 7 
- .0011± 7 
+ .0112i:28 
+ .0052±10 
+ .0025±12 

-f0.0103:kl0 
+ .0105i:18 
-f .0081±19 
-f .0106±10 
+ .0080±10 


±142 


Mizusawa B 


±105 


Tschardjui AB 


db427 


Cailoforte A 


±198 


Carlof orte B 


±172 


Qaithfflisbarg AB 


±157 


PInclnnftti A.. . ,. 


±315 


Cincinnati B 


±280 


UkiahA 


±197 


UkiahB 


±151 







The agreement between the correspondmg items of Tables 8 and 9 
seems to be entu*ely satisfactory. As was to be expected, Cincinnati 
and Tschardjui have probable errors considerably larger than those of 
the other stations. The fact that at all stations the increase in 
latitude is more rapid toward the end of the period than toward 
the beginning appears in two ways: First, the values of z are, without 
exception, greater than the values of y for the same station and, 
usually, much greater: second, when z^y, the ''A'' values of y are 
always smaller than the "B'' values, or in other words, the longer 
the period covered the greater the average rate of increase, so that 
tljLC rate must itself be increasing as time goes on. This state of affairs 
roquires that, when the stations are taken together and y and z are 
equal, the same period of time should be involved for all stations: 
this means that either the ''B" values should be taken when all 
stations are included or that it the "A'' values are taken, only the 
three stations, Mizusawa, Carloforte, and Ukiah, should be used, as 
these three alone nm from 1900 to 1917, inclusive, Cincinnati stopping 
with 1915. When y and z are unequal the ^' A" and the *^B'' values 
of y are practically identical. This is as it should be, for the addition 
of a few more observations after 1914 should have little effect on the 
fitting of a straight line to the observations between 1900 and 1911 
or 1912; the only mathematical connection between the equations 
in which y occurs and those in which z occurs is the common value 
of X, The values of z are too poorly determined, as appears from their 
probable errors, to afford any very definite information. 

SECTION 6. EXPLANATION OF THE CHANGE OF LATITUDE BY A MOTION 

OF THE POLE. 

It appears at once from Tables 8 and 9 that the size of the increase 
which Prof. Lawson found in the latitude of Ukiah, 0f0094 a year, 
or the rate of Of 0081 found by a better method from Lawson's data, 
is in no way exceptional. 

If we adopt Lawson's explanation, namely, a creep of the surface 
strata, and apply it to all stations, it is evident that the accurate 

89016°— 22 3 
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determination of astronomical quantities will be attended with very 
great difficulties ; it is hardly likely, moreover, that all stations except 
Slizusawa should show a tendency to creep toward the pole, and 
that even Mizusawa should promise to develop such a tendency. 
It is more natural to seek an astronomical explanation, particularly 
in the astronomical quantities that enter into the computed declina- 
tions on which the ODserved latitudes are based. 

The reductions of a star from mean to apparent place include 
corrections for precession, nutation, aberration, and proper motion. 
The annual precession of a star in dechnation is of the form n cos a, 
n being practically constant and equal to about 20 seconds and a 
being the right ascension. Since the star program of the International 
Latitude Service is pretty uniformly distributed in right ascension 
and since the right ascension of the stars observed depends upon the 
season of the year, the effect of an error in n would appear as a term 
of annual period with an increasing amplitude and would be combined 
with other effects having the same period. The annual part of the 
variation of latitude has, however, been eliminated. The effect of an 
erroneous constant of aberration would be to rive rise to an annual 
term with a smaller term of a six-month period, which would go out 
along with annual terms when means were taken by calendar years. 
An error in the constant of nutation gives rise to something like an 
annual term, though one varying slowly in amplitude and phase with 
the revolution of the moon's node. The effect of such an error should 
nearly disappear when the mean over a year is taken, and the mean 
of the small residual effect should nearly vanish in a period of 18 
years.^^ The terms of short period in tne nutation are practically 
without effect even in series much shorter than those here dealt with. 

The effects of errors in the precession, aberration, and nutation 
are then practically periodic and are eliminated from the y^s and z's 
of Tables 8 and 9. An error in the proper motion in declination would 
enter directly into the observed latitudes for its full amount. If 
the average annual proper motion used in computing the declinations 
were too great, the observed latitudes would show an annual increase 
equal to the error. Since the star program is the same for all stations, 
all stations would be affected aUke By such an error in the proper 
motion, and indeed by any error in the star positions which may have 
been overlooked in the previous discussion, except in so far as the 
weather may cause the stars actually observed to be different at the 
different stations. 

The values of y, however, are far from being alike, but show signs 
of being dependent on the longitude. When y and z are unequal, y 
starts with a negative value at Mizusawa and becomes nearly zero 
at Carloforte, the value for Tschardjui with its big probable error 
counting for little either one way or the other; on the American 
continent y is clearly positive. The same algebraic increase from 
Japan through Europe to America appears when we consider the 
values of y for the case when y and z are equal, though the value of y 

u The literature of the variation of latitude abounds in articles dealing with the determination of the 
constant of aberration from observations for the variation of latitude, a subject closely connected with that 
of the Kimura term. The nutation constant has been deduced from the observations of the Intematiooal 
Latitude Service by Prcbyllok in a work entitled: " Die Nutaticmskonstante abgeleitet aus denBeobadi- 
tungen des Intemationalen Breitendienstes'' (publication No. 36, new series, <n the Zentrat-bureaa dcr 
Inteniati(Hxalen Erdmessung). Berlin, 1920. 
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starts from a different value at Mizusawa. From a purely geometrical 
point of view the simplest way to make y a function of the longitude 
IS to attribute it to a motion of the pole. 

To formulate this idea mathematically let us consider the North 

Pole and take as origin of coordinates some mean position of the 

pole and as axes of x and y two lines in a plane tangent to the earth's 

surface at the origin and lying in the meridian planes of Greenwich 

axid of 90® west longitude, respectively." Let u and v be the annual 

rates of displacement of the pole resolved along the axes of x and y, 

respectively; let w be the annual increase of latitude common to all 

stations alike and due presumably, in great part, to errors in the 

assumed proper motions, but including any other effects that may 

produce such a change common to all stations; further let X denote 

the west longitude oi a station. Then the increase in latitude at a 

station is equal to Wy the portion common to all stations, plus the 

component of polar motion toward the particular station. Stating 

this m an equation we get 

u cos X+v sin X+tt7=y. (4) 

This is the type of observation equation which may be used to deter- 
mine values of u, v, and t^ so as to represent, as nearly as possible, 
the observed values of y. The quantity z might be substituted 
for y; it is evident that z corresponds to a different program of stars 
so that the w obtained from the z's should differ from that obtained 
from the y's. The values of u and v, however, should be the same 
in both cases, on the supposition that the progressive motion of the 
pole continues unchanged. Unfortimately the probable errors of 
the ;s's are so large as to make a calculation of this sort of rather 
doubtful value. In one case, however, the s's have been used as a 
partial verification. (See p. 36.) 

As is evident from the probable errors in Tables 8 and 9, the values 
of y are not all ecjually reUable; these probable errors are some guide 
as to the appropriate weights for the observation equations of form (4) , 
'but they are not necessarily an infaUible guide. In forming the 
normal equations from the observation eauations (4), two systems 
of weighting were adopted; one in which the weights were taken, in 
accordance with the theory of least squares, to be inversely propor- 
tional to the squares of the probable errors ^* of the y's m Taoles 
8 and 9, and the other in which the weights were arbitrarily taken 
as unity except for the stations Tschardjui and Cinciimati, to which 
were assigned weights of one-half. Since there are two systems of 
weighting, two periods according to which means are taken, the 
annual and the 14-month, and two values of y according as y and z 
are equal or unequal, there are eight possible solutions for u, v, and w. 
All these solutions were carried out and the results with probable 
errors attached are set down in Table 10. 



i> These are the axes used in the publicatiaiu of the International Latitude Service. 
14 That is, approximately. Convenient round numbers that satisfied this cimdition fairly well were 
used fen* the weights. 
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Table 10.— ifbfion o/ North Pole. 
[Probable errors (indicated by a i: sign) are in units of the fourth decimal place of seconds.l 



Solution No. 



1 
2 

3 
4 

5 
6 

7 
8 



Annual component of south- 
ward motion of North 
Pole along- 



Meridian of 

Greenwich. 

tt 



It 

+0.0008d: 9 
+ .0014d:15 

-f .0002±11 
+ .0016±16 

+ .0007±16 
+ .0015d:20 

+ .0002±15 
+ .0018±21 



Meridian of 
90* west. 

V 



tt 

+aO0G9:fe 9 
+ .0050^13 

+ .00&5d:ll 
+ .0053±13 

+ .0043i:14 
+ .0034±17 

+ .0047±13 
+ .0036±18 



Correction to 
average an- 
nual proper 
motion of 

stars in decli- 
nation. 

w 



•t-a0019d: 7 
+ .0043i:10 

+ .0013± 8 
+ .0040±10 

-h .0026^:11 
+ .0052±13 

+ .0022±10 
+ .0050±14 



Probable 

error of sta- 

ti<ni of unit 

wei^t. 



±12 

±23 

±21 

±20 

±23 

±» 

±22 
±30 



SoXvXixm No, /.- 
years. 

Solution No, 2.— 

Solution No, S,- 
periods. , 

Solution Jvb. 4.- 
X>eriods. 

Solution No, 6,- 
calendar years. 

Solution No, 6.- 
calendar years. 

Solution No, 7.- 
14-month periods. 

Solution No, 8,- 
14-month periods. 



CONDITIONS. 

-Weights according to probable errors in Table 8, y and z unequal, means by calendar 

Weights according to probable errors in Table 8, y and z equal, means by calendar years. 
-We^hts according to probable errors in Table 9, y and z unequal, means by 14-inonth 

-Weights according to probable errors in Table 9, y and z equal, means by 14-month 

-Tschardjui and Cincinnati half weight, others unit wei^t, y and z unequal, means by 

-Tschardjui and Cincinnati half weight, other unit weight, y and z equal, means hj 

-Tschardjui and Cincinnati half wei^t, others nit weight, y and z unequal, means by 

-Tschardjui and Cincinnati half weight, others unit weight, y and z equal, means by 



The value of e is given in Table 10 in order to aflford some idea of 
how well the y's of the individual stations are represented. In all 
solutions either Gaithersburg or Ukiah, or both, had imit weight and 
the probable error of their y's may be compared with e. It is seen 
that the e^s in Table 10 are, m general, larger than the probable errors 
of the y's in Tables 8 and 9. This is to be expected, since the probable 
errors of the y's come simply from the failure of the latitudes to con- 
form exactly to a linear law of increase, each station being treated 
bv itself and not bemg tied up with other stations by the hypothesis 
01 a polar motion. 

The stations that show large residuals on any scheme of weighting 
are Tschardjui and Cincinnati. The former is so irregular at alTtimes 
that it is difficult to make any plausible suggestion that will bring it 
into line with the other stations, and our uncertainties are increased 
by the enforced change in the position of the station and the com- 
paratively weak connection between the old and the new stations. 
Although Cincinnati shows an increase in latitude between 1900 
and 1912, it is noticeably smaller than its neighbors, Gaithersburg 
and Ukiah. An examination of figures 7a and 76 shows that the 
latitudes of Cincinnati for the years 1902, 1903, and 1904 seem abnor- 
mally high, and these years being near the beginning of the period 
treated tne high values would tend to bring down the rate of increase 
in latitude to its apparently low figure in Tables 8 and 9. A similar 
phenomenon, though not so mareed, occurs at Gaithersburg and 
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XJkiah, as appears from figures 8a and 9a, but it is not easy to see 
any counterpart for it in the old-world stations such as womd exist 
if it were due to the motion of the pole. It may then be ventured 
that the phenomenon is due to variations in refraction, stronger at 
Cincinnati than at Ukiah or Gaithersbui^ because the climate at the 
former place is more markedly continental than at the stations nearer 
tlie sea coast. The suggestion of refraction is only a guess and does 
not justify any arbitrary emendation of the observed results at Cin- 
cinnati. It seems, however, appropriate to call attention at this 
point to the importance of refraction when dealing with such minute 
quantities as we are here concerned with and to the existence of large 
and persistent anomalies in the latitude which it is difficult to attribute 
to anything except refraction, however the latter may have been 
brougnt about, rrzbyllok, a careful student of the subject, saj^ :^^ 

The existence and behavior of the zenithal refraction is of great significance in the 
determination of fundamental astronomical constants. The observations of the 
International Latitude Service have shown that the mean latitude of an evening 
jnay be in error by nearly a second of arc. From the figures in column 6 of Table 2 
it may be seen tliat latitude determined bv a month's observations may be affected 
by an error of nearly Y^; even the mean of a year's work may give values erroneous 
by about Of 1. 

Since we are especially interested in the latitude of Ukiah, it may 
be reinarked that the values of u, v, and w, in Table 10, always rep- 
resent the value of y satisfactorily, the difference between the '^ ob- 
served" in Tables 8 and 9 and the ''adjusted" value being often 
less than the probable error in Tables 8 and 9 and never greatly 
exceeding it. The adjusted rate of increase at Ukiah deduced from 
u, Vj and w is always less than the observed rate, except in one case. 
This is due to the mfluence of its neighbor, Cincinnati, which has a 
small rate of increase. 

It is not easy to choose among all the solutions in Table 10, but in 
order to present a single definite statement as the outcome of this 
part of the investigation, we may arbitrarily take the mean values 
of u and v from all eight solutions. These means are: 

u=+orooio 

v= +0.^0047 

The resultant of these two component motions is in magnitude 

V(0f0010FT(0f0047)^ = 0ro048, or 

in round numbers, Of 0050 = 1/200". 

The direction of the resultant motion is southward along the meridian 

X = tan~M jjTr^^wTTT^ 1 = 78° west of Greenwich. The North Pole then 

appears to be shifting toward the North American Continent in the 
general direction of the east coast of Hudson Bay, or of Washington 
City, at the rate of 0?0050 (or in linear units, 6 inches or 15 centi- 
meters) a year. A mean value of w would have no general interest 
or significance, since it depends on the particular stars in the program 
of observation and the proper motions assumed in calculating their 
declinations. 

The values of w correspond to the increase in the nonperiodic por- 
tion of the Kimura or z term, as obtained by the International Lati- 

u Translated from Astronomische Nachrichten No. 4840-1, toI. 202 (1916), pp. 286-287. 
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tude Service. This increase may be seen by comparing the following 
mean values for z for various calendar years as derived from the vari- 
ous publications of the service.^® 



Year. 


Mean 2. 


Year. 


Meanz. 


1900 


-0.006 
- .010 
+ .028 


1913 


It 


1905 


1915 


+ .oa 


1910 


1917 


+ .114 









It is, of course, evident that the rates of change at Gaithersburg are 
the main support of a conclusion of this sort. Ukiah is suspected 
of being on imstable groimd and unless the deduced motion of the 
pole has the same general character regardless of whether UHah is 
mcluded or not, the conclusion rests on a questionable basis. An 
inspection of the values of y in Tables 8 and 9 suggests that, in order 
to accoimt for the change in latitude at Gaithersburg, even ^^heu 
Ukiah is omitted, the values of u^ v, and w must in general be similar 
to those in Table 10. Actual calculation confirms this view. 

In order to use only the same periods of time at all stations the 
results in Table 10 are based on tne ''B" values in Tables 8 and 9, 
when y and z are assumed to be equal. When y and z are assumed to 
be unequal it is almost indiflFerent whether the "A^' values or the ** B" 
values are used. The results in Table 10 are, therefore, average 
results for the years 1900-14, inclusive, the period corresponding to 
the *'B" values. Since Cincinnati extends only through 1915, the 
only considerable extension of the period covered will depend on the 
results at the three stations, Mizusawa, Carloforte, and Ukiah. These 
three are lUst sufficient to determine u^ v, and w. 

Two solutions were made with these three stations; one using the 
means of the ''A'' values of y from Tables 8 and 9, where y and z are 
assumed equal; the other using the mean of the "A" values of z 
from Tables 8 and 9. 

The first result is then a mean one for the years 1900-17, inclusive. 

^^^ u= +0.^0007, 

t?=+0f0066, 
-M;=+0f0055. 

The second solution from the 2's evidently applies to the years 
1915-17 and may be considered as an extension of the results of 
Table 10; it is ^=+0^0032, 

t;=+oroo62, 

w;= +0.^0160. 

The resemblance of the two sets of values of u and v to each other 
and to the corresponding quantities in Table 10, suggests that the 
polar motion continued at about the same rate and in about the 
same direction dining the years 1915-17 as during the years 1900-14. 
The values of w are different, as might be expected from the difference 
in the programs of stars observed. The two sets of values just 
given have, of course, the weakness of depending on the suspected 

u The I does not have the same meaning here as in Tables 8 and 9. 
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station of Ukiah. This is unavoidable if the results are to be ex- 
i^ended through the year 1917, and the chief value of the values just 
given lies in then- resemblance to corresponding quantities in Table 10 
t>hat are not essentiallv dependent upon Ukiah. 

For a revision of the polar motion deduced in this section see 
paragraph at end of introduction (p. 2) added while the pubhcation 
-was m press. 

SECTION 6. THE SIGNIFICANCE OF THE POLAR MOTION. 

The rate and direction of the polar motion, as just determined, 
agree roughly with Wanach's " determination. He found from the 
oDservations of 1900-1911 a motion nrobably not exceeding 0.*0030 a 
year and in the general direction of Newfoundland." The rate here 
determined is more than half as large again as that of Wanach and 
the direction diflFers by 22°, but in view of the additional observations 
included in the present discussion and the great difference in the 
method of treatment, the difference in the results for so small a 
quantity is not particularly surprising.** 

Wanach remarks that the polar motion could be accounted for by a 
combination of oscillations of long period as satisfactorily as by a 
secular (progressive) change. In astronomy secular changes are, 
in general, merely small parts of periodic changes whose periods are 
very long; thus the distinction between secular and periodic motion 
is chiefly one of convenience. Probably, however, Wanach means 
to say that the polar motion could be accounted for by a combination 
of periodic motions whose periods might be lon^ in comnarison with 
a year or a 14-month perioa but short in comparison with the periods 
oi the secular perturbations of astronomy, periods of perhaps from 
10 to 60 years. Even with the larger motion here found and the 
longer period covered, we can not be sure that the suggested combina- 
tion of motions of moderate period is not the true explanation. The 
Question would then arise as to the origin of these periodic changes; 
tne answer might be that the forces causing these changes were of 
meteorological origin and the periods were the rather indefinite 
peiriods, such as those of 6, 11, or 30 years, that climatologists are 
studying. With our present knowledge, the representation of the 
polar motion, as proportional to the time, and m one direction, is 
simply an empirical formula; extrapolation from such a formula is 
an uncertain process. The geophysical significance of such a motion 
will be touched on shortly. 

It would be quite possible, however, for a polar motion of this 
mf^nitude to have been going on since the beginning of historic time 
without producing any perceptible change in climate. In 10,000 
years the change m latitude would be at most 50 seconds, or about 1 
statute mile, and this only in the regions like eastern North America 
and Mongolia, which are near the meridian of maximum change, the 
meridian being considered as a complete circle. In western Europe, 

H Resultate des Intemationalen Breitendienstes, Vol. V, pp. 219-220. 

i> This is Wanach's tann. of statement. The longitude of central Newfoundland may be taken as 56" west. 
The details of Wanach's method for obtaining the direction and velocity of the polar motion are not given. 

>• Wanach connects the stations with one another from the start in order to obtain the coordinates of the 
pde and then studies the variations of these coordinates. In this investigation, each station is treated 
separately to obtain its rate of change of latitude, and the stations are combined omy at the last in order to 
obtain the motion of the pole. Wanach's treatment gives full weidit to Tschardjui and Cindnnati, instead 
of half weight or less, aduterenoe that results m deo^tsiDgthe deouoed motion of the pole. 
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which lies at nearly right angles to this meridian, the change in latitude 
would be nearly zero. It is, however, precisely in western Europe 
that we have nearly all our long series of accurate astronomical 
observations. It is not improbable, therefore, that these long series 
would have sufficed to bring to light such a motion of the pole as the 
one here supposed had they been taken somewhere else^ but being 
where they are, they have been insufficient. 

This is a convenient place to give the formulas for the change in 
latitude of a point ana the change in the diflFerence of longitude 
between two points, due to the motion of the pole. Suppose that in a 

fiven interval of time the North Pole has moved toward the Equator 
y a small amount 6 and along the meridian whose longitude is a. 
liet (f> and 0' denote the latitudes of a place at the beginningand end 
of the motion, respectively, and let X be the longitude. Xhen the 
change of latitude is given oy 

</>'-0 = ^ cos (X-a). (5) 

Let Xj and X^ denote the longitudes of two places at the beginning of 
the motion, </>3 and <t>i the corresponding latitudes and X'^ and \\ the 
longitudes at the end of the motion. Then the change in difference 
of longitude is given by 

(X'2 — X'l) — (X2 — Xi) == ^ [tan (t>2 sin (Xj — a) — tan <t>i sin (Xi — a)]. (6) 

Equations (5) and (6) come from the differential variations of spheri- 
cal triangles and applv only when ^ is a small quantity and when the 
distances from the pole of the places concerned are very large as com- 
pared with 6. In the right-hand sides of (5) and (6) the primed 
(quantities may be substituted for the unprimed. The left-hand 
side is expressed in the same unit as d. 

Prof. Wegener ,^^ of Marburg, has suggested that a comparison 
of older and more recent observations on differences of longitude 
shows the difference of longitude^® between Europe and North 
America to be increasing. Formula (6) shows that a motion of the 

?ole along the meridian of 78° west woiild have an effect of this sort, 
he magnitude of the effect for a motion of 0f0050 a year is, how- 
ever, much smaller than the difference that Wegener finds, being 
only about one-thirtieth as much. 

In order to see whether the motion of the pole from 1900 to 1917 
is a continuing phenomenon or not, it would be desirable to have a 
careful investigation made of existing long series of latitude observa- 
tions, particularly those made in -fmerica and India. Too much 
should not be expected from such series, as the periodic variations 
are large enough to mask the secular variation unless the former 
have been carefully studied and eliminated. The results at a single 
observatorj, unless obtained from circumpolar stars observed at 
both culmmations, are so dependent on the declinations and on the 

1* Petermann's Mitteiluneen, vol. 58 (1912), p. 307. This is a paper running through several numben 
of Petermann, which was afterwards elaborated in book form under the title Die Entstehung der KoDti> 
nente und Ozeane, of which a second edition has recently appeared, but which the author of this investi- 
gation has not seen. 

*"The difference of longitude Washington-Paris determined by wireless in 1913-14 is 0*09 greater than 
the one dependent upon the trans- Atlantic cable determination of 1866 and 0*05 less than the one in 1892. 
(See Publications of the U. S. Naval Observatory, second series, Vol. IX, pp. £ 98, E 99, and U. S. Coast 
and Cieodetic Survey Report for 1897, App. 2, The Adjustment of the Longitude Net of the United States 
in Connection with that of Europe, 1866-1896.) A study of the various trans-Atlantic longitude results as 
given in the above reports seems to show that little information of value can be obtained from tbem which 
will throw light on the motion of the pole or the creep of the continental blocks. 
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alues adopted for the fundamental constants of astronomy, that 
leir value is doubtful unless discussed with extreme care and there 
airways the possibility of reasoning in a circle — ^from the latitude 
3 tlie declination and from the declination back to the latitude. 
. strong plea for observations on circumpolar stars at both culmina- 
Lons is made by Prof. Chelli,^^ of Turin Observatory. 

It is unfortunate that the number of observatories of the Inter- 
national Latitude Service should now be reduced to three. This is 
ust enough to determine the quantities u, v, and w, but there is no 
lieck from any supernumerary station. At least four stations are 
»ssential to any satisfactory determination and the original number 
>f six, or an even larger number, would be highly desirable. 

The hypothesis of tne displacement of the pole within the earth 
las been put forward from time to time to explain the observed 
'acts of geology and paleontology. Often the hypothesis has 
i^simved forms that appear decidedly extravagant to astronomers 
a,nd mathematicians.^^ It is, therefore, appropriate to call atten- 
tion to the physical consequences of even the comparatively moderate 
polar displacement of Of0050 that has just been deduced. The 
subject oi a possible displacemnt of the pole has been treated mathe- 
matically by several writers, and most satisfactorily perhaps as re- 
gards secular changes of the pole, by Darwin.^ 

The eflfect of sucn a progressive motion of the pole on its periodic 
motion is small, amounting at most to a change in the amplitude and 
phase of the periodic motion necessary to produce a change of about 
OfOOl in the position of the pole.^* With observations of the present 
degree of acciu'acy, a quantity of this size can not be separated from 
the other elements that enter into the periodic motion. 

Let us consider the effect on the position of the pole of a transfer 
of matter from one portion of the surface to the other, with especial 
reference to the idea that such a transfer, due to the erosion and 
deposition of matter by the rivers of the world, might explain the 
displacement of the pole. Suppose a mass m to be transferred from 
a point on the earth s surface, whose latitude is </>i and whose longi- 
tude is Xj, to another point on the surface whose geographic coor- 
dinates are <t>2 and W let c be the earth's mean radius and G and A 
the maximum and minimum moments of inertia of the earth about 
axes through its center, and let B be the angular displacement of the 
earth's axis of figure, that is, of the axis about which the moment 
of inertia is a maximum, in short, the displacement of the pole, due 
to the transfer of mass m; then 6 is given by ^* 

^^ 2iC-A) ^^^^ 2 01 + sin^ 2 </)2 — 2 sin 2 </>! sin 2 </)2 cos (X,— XJ. (7) 

But from the known eUipticity of the earth, (7-^ = 0.00110 Mc^^ 
approximately, M being the mass of the earth, so that we get, numeri- 
cally, 

455 m y-5 

^== —y V sm^ 2 01 + sin^ 2 02 — 2 sin 2 0i sin 2 02 cos (X2 — \) . (8) 



11 Agtrononrii8che Nacfarichten, No. 4631, vol. 103 (1913), p. 405. 

u See an article by Uie late Prof. Joseph Barrell of Yale in Science, vol. 40 (1914), p. 333. 

S3 The influence of geological changes on the earth's axis of rotation, Philosophical Transactions, 



vol. 167 (1877), p. 271, or Scientific Papers, Vol. m, p. 1. 



Parti, 



s«Seep. 

» Thisformulais readily derived fromsection 11 of Darwin's article referred toin footnote 23, above. 
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The above value of d is in radians and must be multiplied bj 
206265 to reduce it to seconds of arc. The maximiim value of ^ 
radical evidently occurs when ^j, Xj, and <fr3, X^ are antipodal points i^ 
latitude ± 45*^. This maximimi value of the radical is then eaual to? 

Some rough estimate may be made of the amount oi mattd 
transported oy the river systems of the world.^ From Grikie| 
Geology (p. 589) it will be seen that it is almost certainly a^ 
overestimate to assimie that the land, as a whole, siiffers an averagi 
annual loss from erosion of 1/1000 foot or 0.3 millimeter; aUowins 
149 million square kilometers to the land area of the earth we M 
for the mass eroded 4.47x10^® cubic meters, or 1.21x10*^ metii 
tons,^ if a density of 2.7 be assumed. Let us suppose all this mas 
to be, originally, m the immediate vicinity of some point in latitu(k 
H-45° ana to be transported to the antipoaal point in latitude —45' 
This is certainly a very gross exaggeration of the effect of the matter 
transported by rivers, for no river carries matter anything like half- 
way around the earth and there is certainly considerable mutual com- 
pensation of the effects of different rivers on the position of the pole. 
Even maldng, however, the extreme supposition just stated, m 
putting Jf= 6.03x10^^ metric tons, we still find for 6 a value (i 
0^0038, a value less than the one deduced from observations. 
The suppositions made are so extreme that it is probable that the 
effect of erosion and deposition is not equal to the tenth part of tk 
value given. 

To the amount of matter eroded by the rivers, there might be 
added the amount of matter removed m mining operations without 
greatly increasing the computed displacement of the pole. From 
Fairly recent statistics the world^s output of coal amoimts to about 
1,200,000,000 metric tons per year, and of iron and steel combined 
about a tenth of that. These two materials make up the bulk d 
the mineral output, so, if we put the total output oi the minerals 
at 2 billion (2 xlO*) metric tons, we shall almost certainly make an 
overestimate, and these 2 billion tons are but a small part of tb 
1.21 XlO" or 121 billion tons estimated to be removed by the rivers 
of the globe. 

Although the amount of matter removed hj the rivers is insuf- 
ficient, even under the most favorable conditions to account for 
the observed <iisplacement of the pole, yet this displacement may be 
accounted for by a widespread redistribution of matter having a 
thickness of the same order of magnitude as the average layer that 
was assumed to have been removed by the rivers. Let us consider 
the approximate expression for the radius vector, r, of the eartii 
in terms of the latitude, </>, 

r^a—afsia^ <t>, (9) 

a being the ccjuatorial radius and / the flattening or ellipticity. 

By differentiation, -, r - r. . j . 

•^ ' ar = — of sm 2<^ d(t>. 

Putting (Z<^ = 0f0050 and using known values of a and /we find for 

dr in rmllimeters , r^ ^^ - r.. 

dr= —0.51 sm 2<^. 



M For data of this sort, see, among other works of reference: GeUde, Text Book of Geology, 4th Ed. 
Vol. I, pp. 487-496 and 586-597. Rudiki, Physik der Erde, pp. 470-485. 
ST Tnisls a very high estimate, equivalent to the effect of about 290 Miasissippis or 1,500 Dannbes. 
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f ^we reckon longitude (X) from the meridian along which the pole 
J assumed to be approaching the equator and suppose that at 
very point in the earth there is added algebraically a layer of 
ariable thickness Ar, in millimeters given by 

Ar = - 0.51 sin 2(t> cos X, (10) 

hen the ellipticity of the earth will remain unaltered, but the position 
>f the pole of figure will shift with respect to the continents by just 
he Of 0050 used in building up the equation. It is as if we took a 
jeographic globe flattened at the poles and squeezed it so that the 
uxLount of flattening remained imchanged, but the flat places were 
2.0 longer under the Arctic Ocean and the Antarctic Continent. 

This simple geometrical calculation impUes other transfer of 
matter than on the smface merely; every equipotential surface, 
within the earth is affected in much the same way as the outer surface, 
though to a diminishing extent as we approach the center. To follow 
out the analogy of the squeezed globe, conceived as solid, it is evident 
that there must be movements of matter in the interior in ord^ to 
permit the globe to take on its new figiu'e with respect to the con- 
tinents painted on its surface. The deformation is evidently a 
plastic one and would leave no record of itself in the rising or falling 
of the continents with respect to sea level. 

The effect on the position of the pole of a layer of matter of thick- 
ness given by equation (10), and of the density of ordinary smface 
rock (<r) may be computed and compared with the foregoing result. 
We assume that the layer is simply added algebraically to the 
radius and simply rests on the surface as on a ngid body without 
causing movement of matter in the interior. For generality, we 
assiune a maximum thickness of h instead of 0.51 millimeter as in 
(10). We then have, ^^ using as much of the previous notation as is 
necessary, 

4 
By using (7— A = 0.001 1 Mc^ and Jf = « ir c^ p, where p denotes the mean 

density of the earth, we get for 6 in radians 

^ = 364--, 
c p 

or A = 0.00275 c ^ -• 

This gives the maximum thickness of the stratum. With -=2.1 

and d = 0f0050 we find ^ = 0.89 millimeter. 

The full amount of the matter added or removed does not appear 
as a rising or sinking with respect to sea level. The addition or 
removal oi matter raises or lowers the sea level. By a fundamental 
theorem in the theory of attraction the apparent rising or sinking 
with respect to sea level will be the thickness of matter added or 

> See sec. 11 of Darwin's article referred to in the footnote on p. 39. 
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removed, respectively (assuming the law of thickness to be of ilj 
form (10), with the general h in place of 0.51), multiplied by 1 — t- 

Op 

that is, by about -j^. 

If the earth be treated as elastic, rather than as plastic or rigid 
it appears that the effect of an additional layer of given tbicknessi 
less than it would be if the earth were rigid.'* The elastic yieldicj 
of the earth causes a rearrangement of matter, and to estimate tlii 
effect of a surface lavor we must, for the earth, apply a reducim 
factor of about | to the effect of such a layer when placed on a ri^ 
body. The foregoing estimates have supposed an actual transfer d 
matter from one part of the surface to another, or internal change] 
equivalent to such a transfer. If the changes are merely cbanges (i 
density in a relatively thin outer shell of the earth, then the eflfect of aa 
elevation or depression caused by an expansion or contraction d 
portions of this shell is less than the effect of an equal elevation 
or depression caused by transfer of matter.^® 

The purpose of the loregoing calculations and remarks regarding 
the effect of the transfer of matter on the position of the pole is merelj 
to give some idea of the order of magnitude of the quantities involved 
and to call attention to some of the elements that enter into the 
problem. Large and presumably widespread elevations and sub- 
sidences of the earth's surface have occurred in the past and changes 
are still occurring, but we know so little about their causes and laws 
that it would probably be difficult to connect known geologic facts 
with a possible motion of the pole. 

SECTION 7. ADDITIONAL OBSERVATIONS AT UKXAH. 

While the manuscript of this pubUcation was beingprepared for the 
printer recent results for Ukiah were received from Dr. H. G. Van de 
Sande Bakhuyzen, secretary of the International Geodetic Associa- 
tion. The transmitting letter is dated September 26, 1921. These 
results are given in Table 10a. 

Table 10a. — Observed latitudes at Ukidh latitude station. 




1918. 



Apr. 10 
Apr. 27 
May 27. 

June 28 
July 30. 
Aug. 31 
Oct. 13. 
Nov. 12 
Dec. 20. 



,39" 08' 12f 24 
12.25 
12.25 

39" 08' 12. 38 
12.34 
12.35 
12.20 
12.23 
12.26 



Date. 


Latitude. 


Date. 


Latitude. 


Date. 


1 

1919. 




1920. 




1921. 


Jan. 20 


39° 08' 12*25 


Jan. 15 


39° 08' 12f32 


Jan. 17 


Feb. 22.... 


12.21 


Feb. 14.... 


12.30 


Feb. 10.... 


Mar. 10 


12.26 


Mar. 14.... 


12.42 


Alar. i<5. ... 


Apr. 4 


12.12 


Apr. 2 


12.28 


Apr. 2 


Apr.28...- 
May28..... 


12.02 


Apr. 29 

May 23 


12.18 




12.01 


11.97 




June 27 


39° 08' 12. 04 


June 17 


39»08'n.94 




July 24..... 


12.11 


July 25 


U.99 




Sept.l 


12.20 


Aug. 31 


n.88 




Oct. 14 


12.22 


Oct. 14 


11.90 




Nov. 17.... 


12.24 


Nov. 17.... 


12.05 




Dec. 20 


12.31 


Dec. 21 


12.22 





Latitude. 



39'08'12fi7 
12.2; 
12.21 
12.23 



2* See article by Jeffreys entitled "Causes contributing to the annual variation of latitude" in Monthly 
Notices of the Royal Astronomical Society, vol. 76 (1916), p. 504. 

vo See sec. 23 of Darwin's paper already cited. It may oe said that Laplace's law of density used by 
Darwin is now believed to overestimate the yielding of the matter of the earth's interior to compression. 
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These values were plotted as explained on page 8 and readings 
Bre taken to obtain nannonic constants from a 6-year series cover- 
g the years 1915-20, inclusive. The curve for the missing period 
itween 1918.04, the last result for Ukiah published in the Astro- 
:>inische Nachrichten, and April 10, 1918, was drawn from values 
>inputed with the harmonic constants as given by the series 1912-17, 
he results of the analysis corrected for the effect of one component 
a another, but not corrected for ^' slope'' are: 



Annual component. 


li-month component. 


Amplitude 
E. 


Epoch 


Amplitude 


Epoch 


II 

0.121 


e 

93.8 


II 
0.201 


e 

276.9 



The means were taken by years and by 14-month periods and 
cleared as by the process described on page 73. The harmonic 
constants used in the clearance were those just obtained. The 
results are: 

Mean latitudes of Uhiah. 



BY CALENDAR YEARS. 



Calendar year. 



1918. 
1919- 
1920. 



Seconds 

of mean 

latitude, 

uncleared. 



II 



12.243 
12. 174 
12.104 



Clearance. 



n 



-0.021 
+0.011 
+0.034 



Seconds 
of mean 
latitude, 
cleared. 



n 



12.222 
12.185 
12.138 



BY 14-MONTH PERIODS. 



Period included. 



1917.664-1918.848 
1918.848-1920.032 
1920.082-1921.207 



Seconds 

of mean 

latitude, 

uncleared. 



// 



12.226 
12. 187 
12.120 



Clearance. 



n 



+0.017 
+0.007 
-0.011 



Seconds 
of mean 
latitude, 
cleared. 



n 



12.243 
12.194 
12.100 



The results show that the latitude of Ukiah after a long period of 
increase is decreasing again, but nothing conclusive can Be inferred 
from this comparatively short period of observations at a single 
station. Whether the result is due to a creep of the ground, to me 
declinations used, to a movement of the pole, or to abnormal refrac- 
tion can not even be plausibly guessed without the observations at 
other stations for comparison. 



Chapter IIL— THE LATITUDE OF LICK OBSERVATORY. 

It does not seem practicable at the present time for tlie U. S 
Coast and Qeodetic Survey to attempt any extended discussion o| 
a possible progressive change in the latitude of Lick Observatort 
In order, however, to put before the reader the facts now availably 
there is quoted below from Prof. Lawson's article * the passage n 
which he deals with the subject. The table of latitudes has oed 
corrected in accordance with a communication from the director oi 
Lick Observatory, Dr. W. W. Campbell, and some remarks, chiefly 
extracts from letters written by him, have been added.^ . 

The extract from Prof. Lawson's article reads : 

Meridian Circle observarions have been maintained at Lick Observatory siiice 
1893 by Astronomer R. H. Tucker, who has kindly supplied me with the following 
memorandum of the values for the latitude of the station. The results are based o 
observations both on circumpolar and zenith stars for an average of sixty-eight night? 
per year during the period from 1893 to 1918, except for the interval 1908 to 1912, 
when Mr. Tucker was in South America. The values tabulated are the averages by 
quarter years, as near as possible, and the number of nights in each quarter durin? 
which observations were made averj^es seventeen. The results have been corrected 
for the revised refractions adopted at Lick Observatory and for the variation ol 
latitude. 

Average qrwrterly values fof the latitude of Lick Observatory based on Meridian CM 
observations by Astronomer R. IL Tucker from 1893 to 1918. 



Epoch. 


Latitude. 


Epoch. 


Latitude. 


Epoch. 


Latitude. 




1893.82 


37" 20' 25f48 


1899.12 


37" 20' 25*86 


1906.11 


37* 20' 25f 58 




1894.12 


.69 


.32 


.82 


.45 


.63 




.33 


.47 


.63 


.58 


1907.75 


.56 




.51 


.44 


1900.38 


.73 


1908.33 


.18 




.83 


.43 


.89 


.64 


1912. 18 


.79 


1895.11 


.61 


1901.46 


5.94 


.51 


.60 


.33 


.57 


.57 


6.03 


.75 


.38 


.63 


.52 


.92 


5.71 


.92 


.78 




.83 


.56 


1902.13 


.90 


1913. 16 


.68 




1896.11 


.82 


.37 


.76 


.53 


.32 




.40 


.72 


.62 


.87 


1914. 14 


.68 




.67 


.43 


.88 


.68 


.37 


.67 




.91 


.52 


1903.12 


.65 


.63 


.62 




1897.13 


.41 


.35 


.76 


.90 


.69 




• .31 


.47 


.60 


.88 


1917.18 


.34 




.72 


.62 


.92 


.25 


.39 


.58 




.92 


.67 


1904.11 


.48 


.62 


.54 




1898.15 


.94 


.30 


.76 


.88 


.51 




.36 


.65 


1905.59 


.26 


1918. 11 


.55 




.83 


.77 


.77 


.42 









I have taken these values and plotted them on coordinates in figure 1, B.' It is 
apparent frcan an inspection of the figure that, notwithstanding the variations in the 
values, there is in general a fairly steady increase in the latitude from 1893 to 1903. 
This increase is expressed in the mean line A-B, It amounts to (y\A in ten years, 
or at the rate of 0''''.04 or 1.24 meters per year.* 

Between 1903.60 and 1903.92 there is an exceptionally large drop of 0^^.63 in the 
value for the latitude. In this interval occurrea the earthquake of August 2, 1903, 
which was rather severe at the lick Observatory. It seems not improbable that 
the rather large drop in the value of latitude may be due to a shift of the ground at 
the time of this earthquake. The values for latitude since 1903.60 clearly fall into 
a grouping distinct from the grouping of the values for the period preceding that 
date. The mean expression for the values, as indicated by the line C-D, is, however, 
not so satisfactory as the line A-B^ owing largely to the interval of no obsorvatioDs 

1 "The mobility of the Coast Ranges of California, an exploitation of the elastic rebound theory:" Uni* 
versity of California Publications, Bulletin of the Department of Geology, vol. 12, No. 7 (Jan. 11, 1921). 

s These extracts are used with Dr. Campbell's permission. 

* Fieure 10 takes the place of figure 1, B. It contatoks Uie corrections communicated by the ditecter 
of lAek Observatory. 

4 See p. 49 of this publication for revised figures obtained by least-square adjustment. 
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from 1908 to 1912. This line rises from 1903 to 1915, but at a less rapid rate than 
the rise of the line A-B, The rate of increase of the value for latitude for this period 
is 0^^.022 or 0.70 meter per year.* Beyond 1915 the data are insufficient for the 
determination of a mean position. It is interesting to note that there is no indication 
of sudden change of position at the time of the slip on the San Andreas fault in 1906. 
The rate of latitude of Lick Observatory on either side of 1903.6 is much greater 
than that at Ukiah, and there may be doubt as to the true value of the rate, owing 
to the fact that the lines A-B and C-D each represent the mean position of a series 
of x)oints which depart notably and irregularly from the adopted mean. There can, 
however, be no doubt as to the general signincance of the observations, particularly 
between 1893 and 1903. They accord with those of Ukiah in pointing to a northerly 
creep of the region. The increase of latitude both at Ukiah and Moimt Hamilton, 
taken together with the legitimate deductions that may be made, in the light of the 
reboimd theory, from the results of the geodetic surveys in the intervening territory, 
seems to establish the fact of a northerly creep of the Middle Coast Ranges. 

With regard to the abrupt change m latitude that seems to have 
occurred between 1903.60 and 1903.92, Dr. Campbell on April 29, 
1921, writes as follows: 

In supplying the data to Prof. Lawson, Mr. Tucker endeavored to base them, as far 
as practicable, upon quarterly mean values* that is, about four latitudes per anniini, 
He was not aware that the date 1903.6 mignt be a critical one in Lawson's problem, 
and the latitudes for 1903 were grouped unfortunately. Mr. Tucker is of the opinion 
that the four data for the year 1903 should be regrouped for only three mean epochs. 
You will see the reason for this from the following tabulation of Mr. Tucker's results 
for that year as grouped by months: 

January, 8 nights, mean latitude 25. 76 

February, 8 nights, mean latitude 25. 61 

March, 2 nights, mean latitude 25. 57 

Mean for 3 months 1903.10 25.67 

April, 6 nights, mean latitude 25. 68 

May, 11 nights, mean latitude 25. 86 

June, 4 nights, mean latitude 26. 16 

Mean for 3 months 1903.37 25.86 

Mr. Tucker was absent on the Atlantic coast from July 16 to September 16. 

October, 2 nights, mean latitude 25. 61 

November, 2 nights, mean latitude 25. 26 

December, 8 nights, mean latitude 25. 25 

Mean for 3 months 1903.92 25.31 

Again in the same letter he writes: 

The plottings for the four data of the year 1903 should be changed to conform to 
the three substitute data. 

Mr. Tucker *8 series of latitude values are unfortimately not on a homogeneous 
system. They are merely by-products of his extensive observations made for the 
purpose of determining the acciurate positions of stars widely distributed over the 
northern sky. The bases of fundamental stars are quite variea, having been selected 
to meet the requirements of the successive programs upon which he was engaged, 
and without any thought whatsoever as to the latitude values which could be extracted 
as by-products. The four quoted results for 1897.72 to 1898.36, inclusive, depend 
upon circumpolar stars observed at both upper and lower culminations. All the other 
results are based upon a combination of certain so-called standard systems and cir- 
cumpolars observed at both culminations. These systems are as follows:* 

Auwers Berliner Jahrbuch+circumpolars 1893.82-1896.40 <^n =25.57 

Boss U. S. Northern Boundary Commission circumpolars 1896.67-1897.31 

inc <^ 4=25.46 

Circumpolars exclusively 1897.72-1898.36 inc 4> 4=25.72 

Auwers 3034-circumpolars 1898.83-1900.38 inc *t> 5=25.75 

Newcomb+circumpolars 1900.89-1904.30 inc <^i8=25.74 

Auwers Berliner Jahrbuch+circumpolars 1905.59-1918.11 inc 02t =25.53 

* See. p. 49 of this publicftticm for revised finires obtained by least-square adjustment. 
B The subscripts otw f9 in4icat« the number of quarter-yearly values in tl\e list on p. 44 contribatloc 
to the latitude, 
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Mr. Tucker is of the opinion that the variations in his tabular latitudes are due in 
part to differences in the fundamental s^tems and in the remaining part to accidental 
errors. My opinion, on the contrary, if I may venture it, is that Lawson has foimd 
evidence of a change in the Mount Hamilton latitudes at a time approximating that of 
our very severe 1903 August 2 earthquake. *« You will note that while the latitudes 
for the two Auwers Berliner Jahrbuch periods are substantially equal in the mean, 
yet the results increase from 1894 to 1896 and again from 1905-8 to 1912-15. The 
results on the Newcomb basis are on the average large, but we can not fail to note 
that they drop suddenly following 1903.6 and that these later vsdues are in agreement 
with the 1905-8 Auwers results. In my opinion, none of the observed variations 
would call for serious consideration were it not for the sudden drop in the values in 
the latter half of 1903 and the conformity of these latter to those of the succeed- 
ing years. The subject at least merits most careful consideration. 

Dr. Campbell takes up the subject in another letter, dated May 12, 
1921, saying: 

Perhaps I ought to repeat, in order to avoid any misunderstanding of interpreta- 
tion, that my belief in a variation of the latitude of Mount Hamilton refers only to 
the 1903 epoch, perhaps as an accompaniment of our severe local earthquake of 
August 2, 1903. The observed latitudes preceding and more immediatelv following 
that date all depend upon one and the same fundamental basis, Newcomb s. I think 
you may find value inTrof. Tucker's results for individual nights of 1903, and they 
are as in the following table.^ 

Mount Hamilton (Lick Observatory) latiiudes in 190S by individual nights. 



Date. 


Observed lati- 
tude, 0. 


correc- 
tion. 


Mean 
latitude 
of group, 


Date. 


Observed lati- 
tude, 0. 


0-0. 

correc- 
tion. 

■ 


Mean 
latitude 
of group, 

00- 


Jan. 5 
8 
10 
12 
13 
. 14 
15 
16 


+ 37 20 24.96 
25.97 
25.48 
26.43 
25.98 
25.43 
26.25 
26.09 


- 0.08 

• 

-0.17 


25.76 
25.61 


Apr. 22 
23 
24 
27 
29 
30 


O 1 II 

+ 37 20 25.64 
25.93 
25.28 
25.97 
25.99 
26.45 


- 0.20 

- 0.19 

- 0.16 


25.68 

25.85 
26.16 


1903.32 

May 1 

4 

5 

8 

11 

12 

14 

15 

25 

27 

28 


25.88 

28.17 
25.77 
25.92 
26.26 
26.44 
26.06 
25.47 
26.44 
25.98 
25.92 
25.96 


1903.03 

Feb. 10 
16 
17 
20 
24 
25 
20 
27 


25.84 

25.83 
25.83 
25.59 
25.26 
26.35 
25.99 
25.57 
25.85 


1903.14 

Mar. 
11 


25.78 

25.61 
25.93 


1903.37 

June 2 

3 

10 

11 


26.04 

26.30 
26.50 
26.29 
26.19 


1903.19 


25.77 


- 0.20 


25.57 


1908.10 






25.67 


1903.43 


26.32 






1908.37 






25.86 



K a In expressing the opinion that Prof. Lawson had found ''evidence of a change" at or about the epoch 
of the sevore earthquake, 1903.6. 1 have not for a moment thought that any actual change could equal the 
value0.''63, which is the vertioalinterval between Lawson's two inclined lines at that epoch. This interval 
wasinterpreted by me to be equal to a mudi smaller actual change plus a term of such size as could easily 
arifle from accidental errors. It is recommended agson that the reader note carefully the observed latitude 
vahies for the individual nights in 1903, comparing those obtained after 1903.6 with those obtained before 
that epoch.— W. W. C, April 4, 1922. 

• For convenience in printing, this table has been slightly leananged fixxm the form given in the letter. 

89016*»— 22 4 
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Mount Hamilton (Lick Observatory) latitudes in 1903 hy individval nights — Continued. 



Date. 


Observed lati- 
tude, ift. 


correc- 
tion. 


Mean 
latitude 
of group, 


Oct. 6 
6 


+ 37 20 25.54 
25.48 


+ 0.10 
+ 0.13 

+ 0.07 


25.61 
25.26 

25.25 


1903.76 

Nov. 10 
25 


25.61 

25.33 
24.93 


1903.88 

Dec. 1 
7 
14 
21 
22 
24 
28 
29 


25.13 

24.60 
25.18 
25.59 
25.08 
25.46 
25.62 
25.30 
24.80 


1903.97 


25.18 


1903.92 






25.31 



The data by nights are not corrected for the variation of latitude as determined 
simultaneously at the international latitude stations. Mr. Tucker has computed the 
corrections for the several months on this account, and the corrected latitudes as set 
down for the several months in the last column are those quoted in my recent letter 
to you. You may want to apply slightly different corrections for the latitude varia- 
tions. It is of interest to note that tne latitudes for a considerable proportion of the 
nights in the quarter preceding the earthquake are above 26 seconds, and for several 
nirfits following the earthquake a number of the results are below 25 seconds. 

Prof. Tucker searched for a refraction or other similar effect, which would be re- 
sponsible for smaller observed latitudes in the fall months than in the first two quarters 
of the year, on the basis of his observations running through about twenty-five years, 
and as a result he has no reason to suspect that such an effect exists. 

There is notliing improbable in a sudden change of latitude as the 
result of an earthquake. Such changes are known to have occurred 
in California during the earthquake oi 1906 .J The dislocation of sm*- 
face strata due to that earthquake is, however, rather smaller than 
the change of latitude of Lick Observatory apparently accepted by 
Prof. Lawson as not improbable. The greatest of the dislocations of 
strata deduced for the earthquake of 1906 is about 6 meters or 0.*20, 
and most of them are considerably less, while the change of latitude of 
Lick Observatory late in 1903 seems to have been from 0.^40 to Of 60. 
Moreover, the triangulation of the Coast and Geodetic Survey gives 
little evidence of such a displacement of Lick Observatory. Between 
1887 and the time of the resurvey, shortly after the earthquake, the 
latitude of Lick Observatory appears from the surveys to be im- 
changed, while the longitude has changed but Of 005 or 0.12 meter. 
Even this change is considered doubtfm by Hayford and Baldwin. 

The main object of this report is not, however, the discussion of 
sudden changes of latitude due to the dislocation of surface strata, 
nor of the geological questions involved in a gradual creep of the 
surface strata; the object is rather the examination of the astro- 

"> Coast and Geodetic Survey Report for 1907. Appendix 3, Earth Movements in CalifomlA Earthqoaka 
of 1906. bv Hayford and Baldwin; also Report for 1910, Appendix 5, Triangulation in CaUfcnnia, by Duvall 
and Baldwin. Other references are given in Prof. Lawson's paper. 

« See references in the preceding footnote; p. 76 of the Report for 1907, and p. 195 of the Report for 1910. 
See also addendum on p. 108. 
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nomical evidence for or against a progressive change of latitude. If 
the increase in the latitudes of Gaithersbiu^, Cincmnati, and Ukiah 
deduced in the preceding chapter is due mainly to a secular shifting 
of the pole, the latitude of Lick Observatory should show a similar 
increase practically identical with that of Ukiah. In the extract 
quoted at the beginning of the chapter Prof. Lawson presents the 
evidence for such an increase even larger than the one at Ukiah. 

Instead of drawing the lines AB and CD by eye we riiay make a 
least-squares adjustment, using observation equations of the form pre- 
viously adopted, namely, 

in which A^ is the difference between the observed latitude and some 
convenient but arbitrary initial latitude, x is the adjusted latitude 
at some arbitrary epoch, t the time reckoned from the epoch, and y 
the adjusted rate of change of latitude. In making these adjust- 
ments the observations of 1903 were regrouped as su^ested by Dr. 
Campbell (see extract on p. 46), and the latitude for 1908.33 was 
omitted in accordance with a further suggestion from hun. The 
values of y deduced from the adjustment are — 

For the period 1893.82 to 1903.37, y= +0!'0361 ±0f0052 per year. 
Andfor the period 1903.92 to 1918.11, y= +0f0040±0.*'0045 per year. 

These lines AB and CD are shown with these rates of change or 
slopes in figure 10, which takes the place of Prof. Lawson's fis[ure 1, B, 

These rates of increase are not very different from those found by 
Prof. Lawson in the passage quoted and are of the same order of mag- 
nitude as tile rate for Ukiah. Unfortunatelv, however, calculations 
of this sort are nearly meaningless unless they cover a long period, 
and this is so because, as Dr. Campbell states on page 46, the obser- 
vations are not based on a homogeneous system of declinations. If 
the period covered were long, the progressive change of latitude would 
then presumably be too great to oe masked by cnanges in the decli- 
nation system, but this is hardly true of the Lick observations. It 
may be noted, however, that the 1 1 quarter-yearly means from 1893.82 
to 1896.40 and the 21 means from 1905.59 to 1918.11 are both on 
the Berliner Jahrbuch system and both give practically the same 
latitude, the mean seconds of the earUer value being 25."'57 and of 
tie later ones 25f53 or 26f55 when the suspected value for 1908.33 
is omitted. 

If we accept a sudden decrease of the latitude due to the earth- 
quake in the autumn of 1903, then the fact that the latitudes have 
become equal again must mean that the sudden decrease has been 
compensated by a gradual increase, which would, on Prof. Lawson 's 
theory, be due to a creep of the surface strata, or might, on the ten- 
tative hypothesis in the preceding chapter, be due to a shifting of 
the pole. In round numbers, this rate of gradual increase would be 
about 0f02 a year, provided we assume a sudden decrease of 0.*4 
in 1903; the rate, of course, would be subject to all imcertainties due 
to errors in the assumed sudden decrease, observation errors, and 
errors due to incomplete identity of the declination systems which 
have been taken as absolutely identical. Perhaps all that may be 
safely asserted is that it would be difficult to disprove the existence 
at Lick Observatory of a progressive rate of increase of latitude of 
the same order of magnitude as the rate for Ukiah. 
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It is questionable whether it would be worth while to try furttier 
refinements by reducing the Lick Observatory results as nearly as 
possible to a common system. Dr. Campbell writes on May 16, 1921 : 

I do not think it is practicable to determine a set of systematic corrections which 
wiU reduce the observed latitudes based upon the various fundamental systems to a 
homogeneous system. Proi Tucker's work as published seems to ^ow the following 
differences in the declination systems: 

Publications, Lick Observatory, VI, p. 13; Berliner Jahrbuch— Boss, U. 8. Northem Boundary C<mi- 

raission,— +or02. 
Publications, Lick Observatory, VI. p. 148; Berliner Jahrbudi—Auwers 303— —0734 at— 13* dec. 
Lick Observatory Bulletin, X p. «; Berliner Jahrbuch— Newoomb— — ori6 at equator 

_s -i-of 08 at+25 dec. 
Lick Observatory Bulletin, X, p. 42; Berlin^ Jahrbucb— Boss Preliminary General Catalogue 

•= +0f 11 at equator. 

Mr. Tucker's result for the quantity (Auwers new Berliner Jahrbuch system minus 
old Berliner Jahrbudi system) equals —Of 02 at the Equator for 110 staro observed. 
The two declination systems are sufficiently identical for present purposes. 

The observed latitudes obtained by giving equal weights to the circumpolars and 
fundamental stars can be reduced to approximate homc^eneity — that is, to the ecarlt/ 
Berliner Jahrbrich system — ^by applying one-half ihe differences as quoted above to 
them. 

The Lick results for individual nights through May, 1904, have 
been published and will be found in volumes IV, VI, and X of the 
publications of the Lick Observatory. 

The difficulties arising from the declinations and proper motions 
of the stars suggest that the two most practicable ways of studying 
the secular variation of latitude, where the changes involved are 
necessarily very small, are (1) a plan substantially the same as that 
of the International Latitude Service, in which all observatories have 
the same latitude within a few seconds and use the same program 
of stars and the same methods of observing and (2) the systematic, 
long-continued observation of circumpolar stars at both culminations. 
It is to be hoped that the work of the International Latitude Service 
will be extended, rather than restricted, and that the observatories 
will be maintained without interruption for long periods. If it 
should appear that the creep of surface strata is a common occurrence 
it will be necessary to select stations where, from geological con- 
siderations, little or no such creep is to be expected. At worst, if 
this were impossible, it might be necessary to greatly multiply the 
number of stations in the hope that, in the long run, the northward 
creep of some stations would nearly neutralize the southward creep 
of others. The plan of observing circumpolar stars at both culmina- 
tions does not require international cooperation, as the other plan 
does. Any observatory with a good meridian circle can do the work. 
It would be necessary to observe on every clear day in order to get 
a good hold on the periodic, or approximately periodic, portion of 
the variation of latitude. The annual term in particular is almost 
certain to be peculiar to the station in question. Part of the annual 
term is no doubt due to the seasonal changes in refraction. Effects 
of this sort might be expected to be large, since one culmination 
would almost always be observed bv dayUght when the other culmina- 
tion occurred at night. Very probably the annual term for one cir- 
cumpolar would be different from the term of another circumpolar 
culminating at different altitudes, so each star would have to be 
analyzed separately. Since available circumpolar stars are not 
numerous, it might be desirable to take circum-meridian observations 
at both culminations in order to increase the number of pointings. 



Chapter IV.— ADDITIONAL RESULTS OF THE DISCUSSION. 

The results that form the principal object of this publication have 
already been gjven in preceding chapters. It was practicable to 
obtain also, with comparatively little extra labor, certain other 
results regarding the variation oi latitude, and these are set forth in 
this chapter. Their chief interest is that, in general, they confirm 
results oTbtained by other methods of discussion or from observations 
covering a shorter time. A brief summary of this chapter and of 
all preceding chapters will be found in Chapter V, page 67. 

SECTION 1. THE HARMONIC CONSTANTS. 

Table 11. — Harmonic constants^ final valties. 
[Epcx^ of 14-month compcment reduced to 1900.00.] 





Years of 
series 
inclu- 
sive. 


Slope in 
seconds 
per year. 

TO 


Harmonic constants. 


Station and longitude. 


Annual compo- 
nent. 


14-month compo- 
nent. 




Ampli- 
tude. 
R 


Epoch. 

r 


Ampli- 
tude. 
B 


Epoch. 

r 


Mi«usawa,X»»-141"08' 


1900-05 
1900-06 
1906-11 
1906-12 
1910-15 
1911-17 
1912-17 

1900-05 
1900-06 
1906-11 
1906-12 
1908-14 
1909-14 

1900-05 
1900-06 
1906-11 
1906-12 
1910-15 
1911-17 
1912-17 

1900-05 
1900-06 
1903-08 
1903-09 
1906-11 
1906-12 
1908-14 
1909-14 

1900-05 
1900-06 
1903-08 
1906-11 
1906-12 
1908-13 
1910-15 

1900-05 
1900-06 
1906-11 
1906-12 
1910-15 
1911-17 
1912-17 


tf 

-0.0067 

- .0067 
.0000 
.0000 

+ .0083 
+ .0091 
-h .0093 

.0000 
.0000 
+ .0233 
-f .0229 
+ . 0186 
-h .0267 

- .0040 

- .0040 
-h .0100 
-h .0100 
-f .0150 
+ .0157 
-h .0187 

.0000 
.0000 
+ . 0117 
+ .0107 
+ .0117 
+ .0120 
+ .0140 
+ .0150 

.0000 
.0000 
-h .0067 
-h .0125 
+ .0154 
+ .0183 
+ .0292 

+ .0142 
-f .0100 
-h .0133 
+ .0121 
-1- .0183 
-h .0171 
+ .0192 


n 

0.105 
.106 
.132 

!ll3 
.101 

.101 
.105 
.130 
.128 
.132 
.144 

.098 
.102 
.118 
.106 
.117 
.093 
.079 

.023 
.028 
.010 
.012 
.058 
.051 
.057 
.072 

.009 
.022 
.027 
.057 
.054 
.052 
.099 

.042 
.038 
.059 
.057 
.087 
.062 
.060 


o 

14.7 
14.9 
10.4 

9.2 
18.3 
14.4 

8.7 

340 
335.9 
325.3 
326.7 
323.8 
322.8 

275.7 
27a 
277.3 
276.9 
269.3 
274.4 
273.4 

196.6 
205.4 
256.2 
163.5 
182.4 
178.5 
191.4 
159.2 

219.4 
217.2 
261.7 
199.3 
196.4 
171.9 
160 2 

84.5 
91.2 
75.9 
68.4 
05.4 
99.3 
98.3 


0135 
.135 
.213 
.226 
.214 
.178 
.170 

.132 
.123 
.217 
.226 
.231 
.232 

.127 
.129 
.209 

.208 
.218 
.189 
.188 

.143 
.144 
.158 
.177 
.222 
.220 
.217 
.220 

.167 
.154 
.159 
.205 
.201 
.199 
.217 

.126 
.130 
.217 
.216 
.211 
.207 
.210 


o 
178.6 


Tschardjul, X= -63° 29' 


177.5 
163.2 
165.3 
176.4 
179.4 
184.2 

103.0 


Carloforte,X=--8°19' 


102.3 

97.4 

97.9 

104.7 

106.8 

39.0 


Gaithersburg,X»+77*'12' 


35.9 
32.1 
36.0 
52.0 
57.2 
62.1 

321.4 


ninnlnnat*, X=-f-84°25'.. ... ...x. 


318.6 
309.5 
310 6 
304.2 
307.0 
316.2 
318.7 

317.1 


Uklah,X=.+123°13' 


311.3 
298.8 
296.8 
298.3 
304.9 
320.9 

276.6 




273.5 
264.9 
258.0 
278.4 
282.8 
286.6 
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The harmonic constants in Table 3 are not the final values^ though 
they are evidently amply accm^ate for the purposes for which they 
are used. The fact that there is in general a progressive increase in 
the latitude of a station affects the narmonic constants deduced in 
the ordinary way, and these reouire a correction, as explained in 
Chanter VI, section 3, page 80. The constants of Table 11 are those 
of Table 3 corrected in this way and are to be considered as the 
definitive values of these constants resulting from this investigation. 
The quantity m, in the notation of Chapter VI, section 3, which 
was used in deducing the harmonic constants of Table 11 from those 
of Table 3, is also given. The quantity m is given for the year as 
unit of time and is the slope of tne line of mean latitudes or annu^ 
rate of increase of mean latitude. It was obtained separately for 
each series by a graphic adjustment of the data plotted m figures 4a 
to 9a and 46 to 96. A straight line was drawn by eye that fitted the 
plotted points of the six or seven years of the series as nearly as possible 
and the slope of this line computed; this gave the value of m for the 
annual component. This value was mmtipUed by 1.184, the ratio 
of the perioas, to obtain the value of m for the 14-month component. 
The values of m should resemble in a general way the values of y or 
in Tables 8 and 9. On accoimt of the character of the observations 
the resemblance is far from close in many cases. 

SECTION 2. THE ANNUAL COMPONENT OF THE MOTION OP THE POLE 

AND THE KIMURA TERM. 

It is of interest to determine how far the annual variations are due 
to a motion of the pol6 and how far to causes peculiar to each station. 
For this purpose we may subject the constants to the adjustment 
described in section 4 of Chapter VI (p. 82), and then study the resid- 
uals of the separate stations. Since the harmonic constants are 
valid only as average values for a given period, only those series 
should be used that cover identical periods of time for all stations 
concerned. 

The adjusted values of the constants of the north polar motion, as 
deduced irom all available series, are given in Table 12, For the 
present the discussion applies only to the portion of the table above 
the double horizontal line; the portion below wiU be treated later. 
This table gives the constants for two methods of representing the 
motion. For the first method there are given the quantities a, 5. 
a, and j3, in the notation of Section 4 of Chapter Vl. The annual 
component along the ineridian of Greenwich is represented by a term 
a cos (/c^— a), where fis the time reckoned from the beginning of the 
year. The quantity k depends on the unit of time used, but for the 
annual component ic^ = 360®, when t is put equal to the number of 
time units in a year. The annual component along the meridian of 
90° west of Greenwich is similarly represented by a term of the form 
5 cos {d — p) . The quantities a and P are the values of d at the time 
of the maximum excursion of the pole away from its mean position 
toward the equator along the meridians of Greenwich and 90*^ west, 
respectively. For example, the value of a being on the average 
somewhat greater than 240°, the maximum southward excursion of 
the pole along the meridian of Greenwich occurs after 240/30 or 8 
months have elapsed; that is, in September. The maximum south- 
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ward excursion along the ninetieth meridian west occurs after i3/30 
months have elapsed; that is, in June. The annual motion of the 
North Pole about its mean position is thus seen to be from west to 
east. • 

As is well known, a harmonic motion of the kind just specified means 
that the pole describes an ellipse. The elements oi this elliptic 
motion are given in the right-hand portion of the table. The columns 
for the lengths of the semiaxes and for the eccentricity are self- 
explanatory. The direction, 6. of the major axis gives the longitude 
of the meridian along which the equatorward excursion has its maxi- 
mum numerical value, A. In longitudes ^±90® the equatorward 
excursion has its minimum numerical value, B. The epoch, /a, gives 
the time of year at which the maximum excursion takes place, /z/30 
being the number of months elapsed from the beginning of the year 
to tins maximum. When the elfipse is nearly circular the values of 
and fi are not very accurately determined, and small changes in a, 6, 
a, and P will make them vary widely. The quantity 6 + /x, however, 
remains nearly constant and nearly equal to a, approaching equality 
•with a as the ellipse approaches a circle. An example of this is the 
series for 1912-17. This is an exceptional series, for the directions of 
the major and minor axes are nearly interchanged as compared with 
their directions for the other series. The value of d + fx is, however, 
almost the same as for the other series. 

Table 12. — Harmonic constants of the anniuil motUm of the pole. 





Harmonic constants for merid- 
ian of— 


ElemeMs of the elliptic motion. 


Yean olseries indnsive. 


Greenwich. 


90' west. 


Semiaxes. 


Eccen- 
tricity. 

E 


For major axis. 




Ampli- 
tude. 

a 


Epoch. 
a 


AmpU- 
tnde. 

6 


Epoch. 


Major. 

A 


Minor. 
B 


Direc- 
tion. 

9 


Epodi. 


1900-05 


It 

a 075 
.079 
.099 
.090 
.119 

.087 
.076 
.087 
.087 


e 

254.3 
256.2 
248.1 
247.0 
240.0 

24a6 
238.3 
250.8 
238.8 


a 047 
.052 
.076 
.072 
.099 

.082 
.063 
.062 
.074 


e 

168.3 
172.1 
164.9 
163.2 
153.2 

155.1 
145.8 
166.1 
152.2 


n 

a076 
.079 
.100 
.091 
.119 

.089 
.083 
.087 
.068 


It 

a047 
.052 
.076 
.070 
.098 

.080 
.076 
.061 
.074 


0.78 
.76 
.66 
.63 
.56 

.43 
.42 
.71 
.54 


e 

+ 4.1 
+ 6.8 
+ 12.1 
+ 12.7 
+ 8.4 

+ 25.1 
+103.5 
+ 7.5 
+ 10.4 


• 
251.7 


1900-060 


251.8 


1906-11 a 


238.9 


1906-I2a 


237.1 


1910-15 b 


233.0 


1911-17 c 


217.7 


1012-17 c 


133.5 


1900-11 d 


245.5 


1900-17 « 


230.0 






1900-11/ 


.088 
.086 


243.7 
243.6 


.076 
.076 


165.3 
152.3 


.092 
.066 


.072 
.076 


.63 

.47 


+ 27.4 
- 4.8 


221.7 


1900-17 P 


247.8 







a stations used. 

6 4 stations used. 

e 3 stations used. 

d 6 stations used. Mean harmonic constants for series 1900-05, and 1906-12. 

« 3 stations used. Mean harmonic constants for s^es 1900-05, 1906-12, and 1912-17. 

/ 6 stations used. Mean harmonic constants for series 1900-05, and 1906-12, with the Kimura term 
included in the adjustment. 

9 3 stations used. Mean harmonic constants for series 1900-05, 1906-12, and 1912-17, with the Kimura 
term included in the adjustment. 

The elliptic elements in Table 12 are far from representing the 
annual components of the several stations, and the outstanding 
residuals have enough interest to be given in detail. This is done in 
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Table 13, which contains the diflferences between the observed and the 
adjusted values of 2J cos f and R sin f , the residuals being taken in tlie 
sense *^ observed '' minus "adjusted." Only the residuals for the series 
in Table 12 above the double line ar^ given m Table 13. The adj us ted 
value of R cos f is given by R cos f =* m cosX+p sinX, the notation 
being that of page 87, and the values of m and p being from the 
equations on page 90. The '* observed '* value of R cos f and R sin f 
are readily f ouna from the harmonic constants of Table 11. Sinularly 
the adjusted value of i2 sin f is given by: R sin f = n cos X-f- j sin X. 
For Drevity the residual in jK cos f, with sign as noted above, will 
be here denoted by v^ and the residual in fi sin f by v,. 

Tablb 13. — Residuals of the individual stations for the elliptic elements of the annual 

motion given in previous table. 









[Units in the fourth decimal of seconds.] 










Years of series 


Mizusawa. 


Tschardjui. 


Carloforte. 


Oaithersburg. 


Cincinnati. 


Ukiah. 


inclusive. 


Vl 


vt 


Vl 


v% 


Vl 


V% 


Vl 


Vt 


Vl 


Vt 


Vl 


Vt 


1900-05 


+563 
+556 
+550 
+510 
+332 
+293 
+260 
+338 


-239 
-278 
-352 
-316 

- 73 

- 97 

- 60 
-144 


+626 
+579 

+574 
+618 


+ 65 
- 22 
-152 
-ISO 


+232 
+253 
+408 
+375 
+444 
+388 
+344 
+450 


-245 
-244 
-231 
-207 

- 91 
-128 

- 79 
-190 


+279 
+290 
+224 
+236 


+ 3 
-20 
-14 
- 6 


+414 
+353 
+233 
+201 
- 2 


- 79 
-131 
-296 
-278 

- 6 


+316 
+321 
+557 
+591 
+328 
+287 
+255 
+331 


— 60 


1900-06 


— 06 


1906-11 


—100 


1906-12 


— 96 


1910-15 


— 62 


1911-17 










— 05 


1912-17 














— 69 


1900-17 














—141 



















The pairs of residuals are seen to have a general resemblance to 
one another at all stations. The residuals for a given station would 
naturally be expected to vary according to the stations used ; the table, 
however, does not immediately reveal an effect of this sort, for the 
latter is not separated from the changes in residuals due to the Varia- 
tion in the harmonic constants from one period to another. The effect 
could be brought out, however, for the first four series in the table by- 
using less than the 6 stations and determining the elements of the 
elliptic inotion with the resulting residuals for tne new elements. 

The Eamura term, or the z term in the discussion of the motion of 
the pole, is a term common to all stations.^ This term consists of 
several parts, one of which has been found to have an annual period; 
in so far as the values of v^ and Vj for a ^ven series remain constant 
in passing from station to station, their constant values may be 
used in forming the expression for this annual portion of the Eamura 
term. It is seen from the table that the assumption of a Eamura 
term constant from station to station, even when only a single period 
is considered, is only roughly approximate. The enect represented 
by the Eamura term is now supposed to be chiefly due to refraction,' 
and there is no reason for supposing such an effect to be constant at 
all stations, except in so far as the climate, the local topography, 
and the structure of the observatory * are alike at all stations. 

1 The variation of thelatitude of a station from its initialvalue is written z cos X+y sin X+z; here x and y 
are the rectangular coordinates of the p<de (in seconds) and z is the Kimura term. 

s Besultate, Vol. V, p. 189. 

* The refraction to which the Kimura term is attributed is probably for the most part room refiractton; 
that is, refraction between the air within the observatory and the air outside. The assumption that the 
displacement of the apparent zenith implied in the Kimura effect is due entirely to the nonhorizontaUty 
ona large scale of the surfaces of equal air density, leads to quite improbable basometric gradients. 
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Table 13 is then presented for its possible use in studying, for the 
individual stations, the annual portions of the latitude variation not 
due to the motion of the pole rather than for the purpose of extracting 
from it a Eamura term common to all stations. If we wish such a 
term, we may take as representative the data on the last line, which 
represents the mean of 18 years at the three stations which have 
observations over the entire period. The mean v^ is +373 and the 
mean V3 is — 158. This indicates that the mean annual portion of 
the Eiinura term for the 18 years and the three stations may be 
written in the form 

2 = -f Of0373 cos Kt- 0^0158 sin d, \ ... 

or « = 0?0405co8 (ice-337?0). J ^^^ 

Here id has the same meaniog as on page 52 ; that is, if the time t 
be reckoned in days from the beginning of the year^ then^ for the 
deeree as imit^ k = 360/365. 26. 

For comparison with this value we have Wanach's detenmnation * 
for the years 1900-11, inclusive, determined from all six stations by 
quite another method than the one here used. For the annual 
part Wanach gets a result equivalent to 

z = 0f047 cos (irf- 336?8) , 1 ,^. 

or 2= +0f0432 cos jce-0f0185 sin O. J ^^^ 

A result more properly comparable with (2) may be had by taking 
the mean values of Vi and V2 from the first and third lines of Table 13, 
thus covering exactly the same stations and the same period as 
Wanach. The mean Vi is +415 and the mean t;, is — 117 so that the 
Eamura term in this case is given by 

2= +0f0415 cos ic^-0f0117 sin d, \ ,^, 

or s= +0f0431 cos (jce-344?3). j W 

Expression (3) is almost identical with that used by Dyson ^ which 

reads ^^ +0f041 cos irf-OfOll sin /cf, 1 ,,, 

or 2=0f0424co8(K«-344?5). jW 

The source of Dyson's expression is not given; it presumably covers 
the years 1900-11. The agreement between (2) and (3) is entirely 
satisfactory when the difference between the methods of treatment 
is considered, and the agreement between (2) or (3) and (1) may be 
considered good in view of the difference in the period covered and 
in the stations used. 

We might determine from the adjustment in addition to the 
constants of the annual motion of the pole the constants of the annual 
portion of the Kimura term assumed to be the same at all stations. 
The formulas for this case are given on page — . This has been done 
for two of the more comprehensive series and the results are given 
below the double line in Table 12. The Kimura terms deduced 
from the adjustment resemble the expressions (1) and (3). Corres- 
ponding to (3), we have for the series 1900-11 

2= +0f0448 cos ic^ -Of 01 53 sin /cf, 1 ,^. 

or 3=0.^0473 cos (ice-341?l), J W 

* Resultate, Vol. V, p. 188. 

•MoatlUy Notices, Boyal Astronomical Society, Vol. 78 (1918), p. 457. 
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and corresponding to (1) we have for the series 1900-17, 



or 



2= 4-0f0381 cos ic<-0f0162 sin d, 
2 = Of 04 14 cos (irf-337?0). 



1(6) 



The values of the harmonic constants below the double line differ 
from the corresponding ones above it chiefly in giving a larger oscilla- 
tion alons the ninetieth meridian (the y axis) and an annual path 
of the poie more nearly circular. It is probable that these values I 
are to be preferred to the values above tne horizontal line since the 
existence of the Kimura term is unquestionable and its approximate 
constancy from station to station is lairly well established. It seems 
more logical, then, to include the Kimura in the adjustment from the 
beginning; this process avoids giving imdue w:eight in the series 
1900-11 to the region around the two neighboring stations, Gaithers- 
burg and Cincinnati. 

At any rate the constants below the horizontal line for the series^ 
1900-11 agree better with other determinations for the same period.. 
For comparison we have two different deterxninations by Wanach*' 
and one by Dyson.^ The results are as follows: 





Harmonic constants. 










Designation of constant. 

• 


Wanadi. 


Dyson. 




(1) 


(2) 


a 


0*0841 
050759 

25592 
16592 


050847 
0507eO 

24590 
16490 


05069 
050TO 

24290 
159?0 


b : 


a 


a 




. 



The aCTeement between the results for 1900-11 below the line with 
the results of Dyson and Wanach is as good as the agreement of the 
latter among themselves; the chief discrepancy is in the value of /?• 

The results for the two series below the line are presumed to repre- 
sent average conditions better than any of the other results in Table 
12 • they are, therefore, used in deducing the motion of the imdisturbed 

Eole of inertia in section 5 of Chapter VI (p. 104). One conclusion to 
e drawn from this discussion of the annual motion is that the 
elements of the annual path of the pole depend considSraUy on the 
method of treatment used, and that it is very desirable^. have a 
sufficient number of stations to get a good hold on the EiniirllLtenn. 
The 3 stations now maintained are just sufficient to determinate^® 
constants of the annual polar motion and the annual portion of ^® 
Emura term, the latter being assumed to be the same at all 3 statioi 
There is no station to furnish a check on this assiunption, a conditio 
of affairs which it is hoped will soon be remedied; the original 
stations were none too many. 

The semiannual portion of the variation of latitude has been \ 
included with the annual as being a mere harmonic of it. The same i 
harmonic analysis that gives the annual terms gives the semiannual 

• Resultate, Vol. V, pp. 211 and 217. 

T Monthly Notices of the Royal Astronomical Society, vol. 78, 1918, p. 456. 
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terms also. Material for forming these terms is given in Table 14. 
Its form is somewhat different from that of other tables of harmonic 
constants. Only the s's and c's are given, not the amplitudes and 
epochs (R^s and f's). This is because the amplitudes and epochs for 
the same station differ widely for different series, as indeed the s's 
and c's do also, but in such a case it is not correct to take the arith- 
metic mean of the amplitudes and epochs for different series to get a 
mean amplitude and a mean epoch for the entire period, while such 
a taking of means is entirely proper for the s's and c's. When a 
mean s and a mean c have been found, the corresponding B and f 
may be found in the usual way. (See p. 74.) 

The quantities in the table are not the imimediate results of the 
harmonic analysis, but have been corrected for ** slope" (see pp. 61 
and 80) in the same manner as the quantities in Taole 11, ana the 
effect of the 14-month component has oeen eliminated. The 7-month 
component, being nearer in period to the semiannual component 
than the 14-monm, would have a much larger effect if it existed, but 
as its existence is quite doubtful (see p. 61) no attempt has been made 
to correct for it. 

Tablb 14. — Harmonic constants for the semimmual portion of the variations of latitude. 





Years of 










Years of 






station and longitude. 


series in- 
clusive. 


9t 




Ct 


Station and longitude. 


series in- 
clusive. 


8t 


Ct 






// 




// 






ff 


// 


Minisawa, Japan, 
\ — 141* 08^ 


1900-06 


-0.004 


( 


).000 


Oaithersburg, Blaryland, 


1900-06 


-0.006 


0.000 


1900-06 


- .003 




.000 


1900-06 


- .008 


- .004 




1906-11 


+ .004 


~ 


.012 




1903-08 


•h .002 


+ .001 




1906-12 


+ .004 


— . 


.013 




1903-09 


+ .004 


-1- .001 




1910-15 


+ .012 


_ 


.019 




1906-11 


+ .003 


- .012 




1911-17 


+ .017 


— 


.018 




1906-12 


-H .007 


- .008 




1912-17 


-t- .014 


^ 


.020 




1908-14 
1900-14 


-f .011 
-f .013 


- .011 

- .OU 


Tschardjui, Turkestan, 


1900-06 


+ .008 


+ 


.005 










X— 63« 29'. 


1900-06 


+ .007 


+ 


.006 


Cincinnati, Ohio, 


1903-06 


-f .001 


— .008 




1906-11 


- .009 




.002 


X-+84" 25'. 


1900-06 


- .004 


- .002 




1906-12 


-t- .001 


+ 


.004 




1903-08 


-f .001 


- .001 




190&-14 


- .006 




.005 




1906-11 


- .002 


- .007 




1909-14 


- .007 


•^ 


.003 




1906-12 
1908-13 


-H .001 
+ .008 


- .002 

- .006 


Oarloforte, Sardinia, 


1900-06 


- .003 


— 


.004 




1910-15 


-1- .013 


- .010 


X- -&• 19'. 


1900-06 


- .003 


— 


.006 












1906-11 


+ .001 


_ 


.006 


Ukiah, California, 


1900-05 


•f .006 


— .007 




1906-12 


+ .004 


— . 


.005 


X=+123» 13'. 


1900-06 


+ .004 


— .008 


• 


1910-16 


+ .005 


_ 


.006 




1906-11 


.000 


- .010 




1911-17 


-1- .012 


^ 


.005 




1906-12 


+ .001 


- .008 




1912-17 


+ .018 




.003 




1910-15 
1911-17 
1912-17 


-1- .012 
+ .008 
-f .010 


~ .002 

- .008 

- .006 



The values of s^ and c, in Table 14 are evidently small and their 
significance is therefore liKely to be obscured by accidental errors of 
one sort or another. There appears to be in general a numerical 
mcrease in these quantities between the first and the last series for 
each station. The values of «, and c, for the later series at different 
stations somewhat resemble one another and do not definitely re- 
verse in sign on opposite sides of the pole, as they would if the semi- 
annual term were due chiefly to the motion of the pole. 

As the results in this chapter are incidental to the main purpose of 
this report, it has not seemed worth while to imdertake any very 
elaborate discussion of the semiannual component. Elements of 



> 
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the elliptic motion might be computed from Table 14^ for the sama 
series of years as are shown in Table 12, but the only series ac- 
tually considered are two representative ones; the first is the 
mean result, for all stations, for the years 1900-11, inclusive; the 
second is the mean result, for three stations, for the years 1900-17, 
inclusive. From the values of a, ft, a, and j8 as found from the for- 
mulas on pages 82 and 86, we may write fpr the semiannual terms in 
the motion of the North Pole along the axes previously specified (p. 33) 
from the series 1900-11 



motion along x axis = 0f0016 cos (2/rf 
motion along y axis = Of 0038 cos (2d 

from the series 1900-17 



311? 6), 
174? 0), 



motion along x axis = Of 0046 cos (2/rf— 351? 2), 
motion along y axis = Of 0010 cos (2irf— 45? 4). 



} (7) 
} (8) 



In these equations, k and t have the meanings already specified for 
the annual component (p. 52). A result that should be comparable 
with (7) is that foimd by Wanach," which covers also the years 
1900-11 and is equivalent to 



motion along x axis = Of 0032 cos (2k^-23°), 
motion along y axis=Of 0031 cos (2k^— 175°). 



} (9) 



Equations (7) and (9) agree satisfactorily, difference in methods 
considered, for the motion &lon^ the jf axis, but rather poorly for 
the motion along the x axis. Ine residuals of (7^ and (8) for the 
separate stations m the sense ** observed *' minus '^ aajusted " are given 
in Table 16. 

Table 15. — Residuals for the semiannual element of polar motion. 



station. 


Residual by for- 
mula (7) — 


Residual by for- 
mula (8)— 




Incs. 


In ft. 


Incs. 


In«s. 




-73 
-22 
-68 
-22 
-16 
-46 
-41 


-11 
+ 6 
+ 3 
-15 
- 5 
+20 
+ 


-65 


+44 


Taohardiui 




Carioforie 


. -88 


+59 


O aithersbuTE. 




Cincinnati ^ . , ^ , ^ . . , ^ , , „ 






Ukiah 


-64 
-72 


+43 


Umn 


+49 







The portion of the semiannual variation of latitude peculiar to 
each station as given by Table 15 is seen to exceed, in jgeneral, the 
portion due to the senuannual motion of the pole as given by for- 
mulas (7) or (8). 

Just as in Table 13, the similarity of the residuals in Table 15 for 
the different stations suggests that a common value of all the resid- 
uals should be represented in the Kimura term; this part of the 
term would be 

i;iCos2ic^4-i;3sin 2irf, (10) 

• Resultate, Vol. V, p. 217. 
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the mean values of v^ and Vj being found as in Table 13. Such a term 
might, like the annual portion oi the Kimura term, be chiefly due to 
some peculiarity of the seasonal climatic cycle common to all stations, 
or it might be due to periodic errors in me declinations used. The 
alteration with the time of the values of S3 and C3, to which attention has 
been called and which is reflected in the difference between the two 
sets of residuals in Table 15, suggest periodic errors in the declina- 
tion as a more probable cause^ for the star program necessarily 
changes with the time, and periodic errors may not have been so 
completely eliminated from the later star places as from the earlier. 
The magnitude of such period errors, as judged by the systematic- cor- 
rections necessary to reduce from one star catalogue to another/ is 
in general considerably lai^er than would be ^ven by an expression 
like (8) with v^ and v, o^ ^he order of magmtude of the values in 
Table 15. No mvestigation, however, of such periodic errors in the 
declinations used by the International Latitude Service has here 
been attempted. 

It does not seem practicable to formulate in a few words the 
results of this investigation in so far as they concern the semiannual 
term, this term being small and likely to be obscured by accidental 
errors. The numerical results are presented in Tables 14 and 15 and 
in formulas (7), (8), and (10), which are given for what they may be 
worth. 

Formulas (7) and (8) do not agree at all well with each other, and 
the discrepancy would lead one to doubt the reality of a semiannual 
term in tne polar motion given by either one of them. Wanach ® 
discusses the two halves of the series 1900-11 separately (by a 
method quite different from the one here used) and finds, as regards 
the semiannual term, a passable agreement between the two halves. 
He thence infers that the semiannual term is real and approximately 
of the value found for the series 1900-11. To be sure, the result for 
the 24-year period, 1890.5-1914.5, does not agree well with the 
other result, but he attributes the discordance to the nature of the 
material available before 1900. 

The change with time in the semiannual terms for the individual 
stations seems pretty well established by an inspection of Table 14, 
also the fact that, particularly for the latter part of the period 1900-17, 
there is a semiannual portion in the Kimura term. The mean value 
of the semiannual portion may be represented by expression (10), 
with the mean values of v^ and V2 from Table 15, or as follows: 
Semi-annual portion of Kimura term = + Of 0072 cos 2 k^ + Of 0049 sin 2 d 

= 0f0087 cos (2 /c<-145°.8). 

SECTION 3. THE 14-MONTH COMPONENT OF THE MOTION OF THE POLE. 

From the harmonic constants of the 14-month component of the 
separate stations we may deduce bv a least-squares adjustment the 
harmonic constants for so much 01 the variation as is due to the 
motion of the pole. This has been done and the result is set forth 
in Table 16. This table has the same form as Table 12, except that 
the probable error of 6 (indicated by the ± sign) is shown in connec- 
tion with 6, This matter will be considered later. 

' See for example Boss's Preliminary General Catalogue of 6188 Stars, etc., Washington (Carnegie Insti- 
tution), 1910, Appendix in. 
• Besultate, Vol. V, pp. 218-219. 
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Table 16. — Hcarmonic constants for the li-monthly motion of the Pole, 
[For notation see pp. 82-80. Epochs reduced to 1900. 00.] 





Harmonic constants for 
meridian of— 




Elements of elliptic motion. 






Greenwich. 


9QP W. 


Semi- 


-axes. 




For maior axis. 


Years of series 














Eccen- 








inclusive. 




















• 


Ampli- 
tude. 

a 


Epoch. 
a 


Ampli- 
tude. 

b 


Epoch. 


Major. 

A 


Minor. 
B 


tricity. 
E 


Direction. 
e 


Epoch. 




M 


e+$i 




If 


o 


II 


o 


II 


II 




e o 


o 


o 


1900-05O 


a 131 


32.7 


a 141 


310.6 


a 146 


ai2s 


a 52 


59. 4 ± 4. 


337.4 


3ft.8 


l900-06a 


.131 


aai 


.139 


307.7 


.145 


.126 


.50 


57. 3 ± 4. 8 


336.7 


34.0 


1905-lla 


.20ft 


23.3 


.217 


293.4 


.217 


.206 


.31 


89.6 ±18.5 


293.7 


23.3 


1906-12a 


.210 


25.9 


.218 


295.6 


.219 


.210 


.28 


94.2 ±22.2 


291.6 


25.8 


1910-15 ft 


.210 


40.6 


.219 


312.2 


.220 


.209 


.32 


74.6 ±16.3 


326.8 


41.4 


1911-17 c 


.179 


47.1 


.208 


310.0 


.212 


.175 


.57 


110.1 ±2.7 


293.2 


43.3 


1912-17 •! 


.179 


52.8 


.205 


312.0 


.213 


.109 


.61 


116. 6 ± 1. 5 


29a4 


47.0 


1900-lld 


.173 


29.8 


.176 


300.2 


.177 


.173 


.21 


80.5 


309.5 


30.0 


1900-17« 


.170 


35.5 


.179 


303.0 


.180 


.168 


.36 


lias 


283.9 


34.2 



a 6 stations used. 

ft 4 stations used. 

c 3 stations used. 

d 6 stations used . Mean harmonic constants of separate stations found from series for 190(M)5 and 1906-11 
uid ^ments of elliptic motion determined from mean constants. 

« 3 stati(ms used. Mean harmonic constants of separate stations found from series for 1900-05, 1006-11 
and 1912-17 and elements of elliptic motion determined from mean constants. 

For comparison with these results there is available the results 
deduced by Dyson and Wanach.® Wanach's result is deduced on 
the assumption of circular motion. The radius and epoch for a motion 
of this sort are given under a and a, respectively, and these should 
be approximately equal to the quantities i{a + b) and i (a+i3— 270°) 
of tne elliptic motion for a corresponding period. 

Harmonic constants. 





Designation of constant. 


Wanach, 1900-11. 


Dyson. 




(1) 


(2) 


1900-05 

o:i27 
ori3i 

3196 
30993 


1906-11 


d 


on 738 


orio7 


Of 206 


6 


Of 213 


or 


2691 


2690 


2594 


a 


294?2 












Dyson's results were translated into the notation of this report 
by the use of his own value for the free period, 432.2 days. Part of 
tne slight difference between the epochs for corresponding periods 
may be due to the difference in the free periods assumed, Wanach's 
value being 432.8 days and for this report 432.5 days. The results 
obtained in this report agree satisfactorily with the results of Dyson 
and Wanach for corresponding periods of time. The ellipses deduced 
from Dyson's values oi a. 6, a, and j8, will be foimd to be of the same 
general character and to nave their axes in the same general direction 
as those for like periods in Table 16. 

It is of interest in this connection to mention a verification of the 
a and a for the series 1912-17 coming from outside the International 

* See references in footnotes 6 and 7 on p. 56. 
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Latitude Service. Observations at Greenwich with the Cookson 
floating zenith telescope give a curve of latitude which, when analyzed, 
yields the following expression for the 14-month term for 1912-17, 

Of 191 cos (irf-49?9). 

The ampUtude differs only Of 01 2 from the corresponding a in Table 
16 and the epoch differs 2? 9, a satisfactory agreement for both. The 
Greenwich results may be introduced into tne determination of the 
14-monthly polar motion, but, as might be expected, their introduc- 
tion does not change substantially the figures m Table 16. 

Since the latitude of a place depends Doth on the direction of the 
earth's axis and on the direction of the plumb line at the place, it is 
conceivable that the plumb line of each station may have an individual 
oscillation of its own in a 14-month period, thus causing variation in 
latitude independent of the polar motion. The observations here 
discussed, however, lend Uttle or no support to this idea. An examina- 
tion was made of the residuals from the adjustment that furnished 
the data for Table 16; that is, of the residuals from the 14-month 
Components of polar motion analagous to those that furnish the 
£imura term for the annual variation. The residuals were found 
to be much smaller than for the annual component and to vary 
greatly from one series to another, giving evidence of being mainly 
due to accidental errors. The resmts are, therefore, not given in 
detail. The magnitude of the residuals was such that the probable 
error of « = iJ sin f or of c = i2 cos J", as obtained from the adjustment 
of the harmonic constants of the 14-month component, may be 
estimated at ± 0f006 or ± Of 007. 

There is one reason for the existence of a 14-month oscillation 
peculiar to each station that may present itself to one considering 
the subject. The motion of the pole changes the centrifugal force 
over the earth's surface, and the ocean endeavors to conform to the 
resulting changes in the equipotential surface, thus producing a 
14-monflily latitude-variation tide which has been investigated and 
detected. If the earth wwe covered with water, this tide would 
have the same character at all stations in the same latitude; but the 
earth, not being so covered, the plumb lines of all such stations are 
not equally affected by the gravitational action of the tidal load. 
A little calculation will show, however, that the difference of latitude 
due to an effect of this sort can not exceed OfOOOl or 0f0002 and 
may therefore be neglected.^® 

The hannonic an^ysis of the 14-month component gives reaxiily 
enough the data for a possible term in the latitude variation with 
a 7-month period. The. theory of the variation of latitude does not 
indicate the existence of such a term in the motion of the pole, and 
since, as regards the local periodic deflections of the plumb line, we 
have n^lected the 14-month term we should naturally neglect with 
greater reason the 7-montli term, which would be merely a narmonic 

i«The latitude of a place is subject to an indirect influence, due to the effect of the displacement of the jxde 
on the direction of the plumb line. The displacement of the pole causes equal changes in the eeocentric 
latitudes of all points in the same meridian, but does not cause equid chan^ in the geoerai>hic latitudes, 
which are the quantities directly measured. From one point of view this slight inequality in the change 
of geographic latitude along the same meridian may be considered as due to the effect on the plumb line of 
the swellings and depressions of the earth tides produced by the displacement of the pole and the attendant 
changes in the centrifugal force. A correction for this fact would be necessary onlv if the motion of the pole 
were to be deduced from two sets of latitude stations situated on widely separated, parallels of latitude; Uie 
amount of the ccMrection necessary to reduce stations to latitude 45* as a standard is of the order of mag- 
nitude (rf OTOOOS and therefore unimportant with observations of the present degree of accuracy. 
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of the 14-month term. The 7-month terms deduced from the analy- 
sis give every indication, in spite of all refinements of computation, 
of being due merely to accidental errors. The results of the compu- 
tation of the 7-month terms are therefore not given here. 

The mathematical theory of the variation of latitude (see p. 93) 
assumes that small changes occur in the position of the eartli^ axis 
of figure; that is, the axis about which the moment of inertia is 
greatest." Changes of this sort of a roughly periodic character 
would occm* because of the unsymmetrical annual distribution of 
snowfall, barometric pressure, etc. These periodic changes cause 
what may be called a forced periodic oscillation of the pole of figure, 
the period being a year. The pole of rotation does not in general 
coincide with the pole of figure, but is subject to a forced osculation 
in an orbit of its own, different from the orbit of the pole of figm^e, 
but having the same period. The pole of rotation will have in 
addition wiat is called a free oscillation in a period determined, if 
the earth were a rigid body, merely by its moments of inertia and 
the angular velocity of its rotation, or determined for om* actual 
earth by the foregomg Quantities and, in addition, the mean elastic 
moduli of the earth. This free period is the 14-month period, 
or more definitely 432.5 days. 

Furthermore, theory shows that if the forced oscillation of the 

Eole of figure were strictly periodic, the exact period or periods 
eing imimportant provided none of them coincides too nearly with 
the free period, and if the two principal equatorial moments of 
inertia were equal, then, whatever the character of the forced oscLQa- 
tion, simple or complex, the free oscillation would always consist of a 
uniform circular motion of the pole of rotation. If, on the other 
hand, the principal equatorial moments of inertia were unequal, the 
forced oscillation bemg still strictly periodic, the free oscillation 
would take place in an ellipse with its major axis lying in the meridian 
of that equatorial axis about which the principal moment of inertia is 
the larger; that is, in the meridian of the shorter equatorial axis of 
thegeoid considered as an ellipsoid of three unequal axes.^ 

With a view to detecting a possible inequality in the equatorial 
moments of inertia the assumption of uniform circular motion was 
avoided in making the adjustment to determine the elements of the 
polar motion given in Table 16. By such a process, of com^e, the 
motion, even if exactly circular and uniform, would come out elliptic 
owing to the presence of accidental errors.^^ In such a case, how- 
ever, there would be no evidence of even approximate constancy in 
the direction of the major axis of the ellipse, represented by ^ iu 
Table 16. The variability of ^ in this table, and the size of some of 
the probable errors," show the presence of considerable accidental 

■ ^ ■ ■ ■ I M ■ l^^i^.^ I I I ^—1 ^l.*.^— — i^^»^— ^1— ^^1^— — ■ Mill I P»«.IM» ■ ■! ■■■■ ■■■■^■1 !■■■ ^ m , I — ■ Ml P ■ I ■■ .1 ■ ■».» 

11 This axis may also be called the axis of inertia, and a pole of this axis may be termed indifferently 
either the pole of figure or the pole of inertia. ' 

u Foe the rdation between moments of inertia, inegualitles in the radii of the ge<^d, and inequalities 
in gravity on the surface of the geoid, see Helmert, Hbhere Geodasie. Vol. II, chap. 2. 

1" Departures ^om perfect penodicity in the free and forced oscillations are necessarily treated as acci- 
dental errors. 

i« The larger probable errors should not be taken as anything more than rude ai)proximations; their 
size is due to the small eccentricity of the ellipse, which necessarily renders uncertain the location of its 
major axis, rather than to the greater inaccuracy of the underlying observations. For these larger wrors 
the assumption made in deriving the formulas, namely, that changes in $ may be treated as dinerentials, 
no longer gives a good approximation for a change in d of the size of the larger probable errors. In comieo- 
tion with all these probable errors, large and small, it should be remembered that thev include neither 
the errors of the assumption that tne variation of latitude is representable by harmonic lerms of the kind 
here used nor the errors of the conventions by which the fictitious observation equations were introduced. 
(See p. 87.) 
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error, nevertheless there seems to be a tendency for the major axes 
to cluster aromid a direction not far from 90° west of Greenwich. 

For comparison with this result there is Helmert's determination 
of the same thing by a very different method/^ namely, from gravity 
observations. (S these observations Helmert made various adjust- 
ments differing in the data used and the weights assigned. The 
determinations of the longitude of the equatorial axis with the 
larger principal moment oi inertia (the shorter axis of the geoid) 
varied correspondingly from 94° W. to 122° W. His adopted result 
was 107° W. with a probable error of ±4°. The results from the 
latitude variation are seen to be affected with a greater uncertainty 
tJian those from gravity observations, but are not inconsistent with 
the latter, confirming tne general character of Helmert's result as far 
as it concerns the position of the shorter axis of the geoid. A rough 
confiima-tion of this sort is far from useless, for uxe character of 
Helmert's result is rather surprising, since any great iuequality in 
the equatorial moments of inertia or equatorial axes seems incon- 
sistent with our generally accepted notions about the condition of 
hydrostatic equiliorium prevailing within the earth, or, in other words, 
with the theory of isostasy. Furthermore, Helmert's whole method 
might be questioned on accoimt of systematic errors in the gravity 
observations arising from the topography and the surface geology in 
the vicinity of the gravity stations and because the gravity observa- 
tions are at present confined to one-fourth only of the earth's surface, 
that is, to the land. 

The difference between the two principal moments of inertia 
and the corresponding difference in the equatorial setniaxes are 

Juite as important to know as the longitude of the shorter axis. 
[elmert's various adjustments yield results for the difference be- 
tween the semiaxes varying from 154 meters to 510 meters. His 
adopted result is 230 meters. The corresponding difference between 
A and B, the principal moments of inertia about axes in the plane 
of the equator is given by 

J5-^ = 0.000024 Jfa^ 

where Jlf is the mass of the earth and a the mean equatorial radius. 
Calling C the moment of inertia about the polar axis, we find 
(5 — .d)/((7~i(J.+J5) = l/46, a quantity which enters into the deter- 
mination of the form of the ellipse of polar motion. Schweydar^^ has 
compute the form of the elhpse of the polar motion corresponding to 
Helmert's adopted result; he finds an elUpticity of 0.016^ that is an 
eccentricity of 0.18, which is smaller than any eccentricity m Table 16, 
though not far from the smallest. The result from the variation 
of latitude as regards the difference between the equatorial semi* 
axes of the geoid seems again to be a confirmation of Helmut's 
result in a general way, though one suggesting a rather greater differ- 
ence than me one that Helmert finally adopted. 

In regard to the difference between the equatorial semiaxes 
two facts should be noted, which point in opposite directions. 
First, the difference between the eauatorial semiaxes varies as 
the square of the eccentricity of the ellipse of polar motion, so that 



^ Neae Formeln fdr den Verlauf der Schwerkraft im Meeresniveau beim FesUande. Sitzangsberichte 
9r Koniglicn Preussischeii Akademie der Wissenschaften, 1915, p. 676. 
I'Astronomische Nachrichten, No. 4855, Vol. 203 (1916), p. 101. 

89016^—22 5 
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an eccentricity like 0.36 in Table 16, which is just twice Schweydar^s 
value deduced from Hehnert, implies a difference between the 
semiaxes of the earth of four times that found by Helmert, or 920 
meters. Second, there seems to be a tendency as series are made 
longer for the eccentricities to decrease, indicating that the larger 
eccentricities are at least partly due to accidental errors. Of course, 
if all eccentricity were due to error of this sort, we should expect 
the directions of the major axes to be scattered all over a semicircle, 
but this is not the case. This tendency toward eccentricities dimin- 
ishing as the length of series increases may be noted in the series for 
1900-11 which gives the smallest eccentricity in the table, smaller 
than either of its component parts, *^ the series for 1900-05 and 
1906-11. The series 1900-17 shows the same tendency, yielding an 
eccentricity less than the eccentricities of two of its three component 
parts,^* the series for 1900-05, 1906-11, and 1912-17. 

By balancing these two opposite tendencies we are led to see in 
the figures of Table 16 a rough confirmation of Helmert's general 
result, which is all that the data now available will yield. Perhaps, 
on the whole, the latitude observations indicate a difference in the 
equatorial semiaxes rather greater than Helmert's adopted result, say 
300 meters or even more, rather than 200. These figures and the figure 
90° west longitude for the direction of the axis of the ellipse of free 
polar motion are both subject to a correction for the unsymmetrical 
distribution of land and water and its effect on the yielding of the 
ocean waters to the centrifugal forces arising from the displacement 
of the pole. (See p. 101.) No definitive results are therefore stated 
until this correction has been further investigated. It appears to 
be rather small. If the number of latitude stations could be main- 
tained at six or more over a considerable period, it should be possible 
to get a result for the figure of the earth having as great an apparent 
degree of accuracy as Sie gravitv observations give, and free from 
the suspicion of systematic error due to the fact that gravity observa- 
tions, being at present confined to land, cover only one-fourth of 
the earth's surface and are much affected by local topography and 
geology. 

SECTION 4. COMPONENTS OF LONGER PERIOD. 

The annual chmatic cycle appears to be the fundamental cause of 
the forced annual oscillation of the pole, though there are difliculties 
in obtaining an accurate evaluation of all the effects involved.^^ 
Other climatic cycles are known or suspected and it is natural to look 
for their effects on the motion of the pole. An extended study of 
the subject would be guite beyond tne scope of this report. All 
that is here attempted is to utiuze a period of observation available 
by extending slightly some work done by Wanach.^^ Attention was 

17 This result is not unreasonable, but arises partly from the difference in direction of the malor axes. 
A simple extreme case will illustrate. Suppose two equal concentric ellipses of small eccentricities "with 
major axes at right angles to one another, and let a curve be drawn with radius vector equal to the mean of 
the two radii vectores . Neglecting the fourth powers of the eccentricity, the new curve will be a circle, i.e., 
an ellipse of eccentricity zero. 

" Tl^ eccentricities, from the first two component parts, as set down in the same table, come from six 
stations, not from three, and are therefore not strictly comparable with the series 1900-17. \ 

I'Schweydar: Zur Erklarungder Bewegungder Rotationspole der Erde; Sitzungsberichte def^reussi- 
schen Akademie der Wissenschaften, 1919, p. 357. Jeffreys: Causes contribirtory to the annual Vl^tion 
of latitude: Monthly Notices of the Royal Astronomical Society, Vol. 76 (1916), p. 499, 

» Resultate, Vol. V, p. 217, 




J 



ADDITIONAL BESULTS OF THE DISCUSSION. 65 

drawn by Angenheister ^^ to periods of 3 and 6 years which he found 
in the heights of the barometer at a chain of 9 stations encircUng the 
globe. He gave also some estimates to show that these inequalities 
m barometnc pressure might produce sensible effects in the motion 
of the pole. There is a period of about 3 years in the behavior of 
sunspots with which the 3-year barometric mequality is presumably 
connected. 

It is not difficult to utilize the monthly readings of the latitude, 
already used in obtaining the annual component, for the purpose of 
obtaining terms in the latitude whose period is some mmtipie of a 
year. To study a possible 6-year term the readings of three consec- 
utive months were grouped together and the mean treated as a 
single reading; this gives 24 readings in a 6-year period, and at 
Mizusawa, Carloforte, and Ukiah, the 18 years of observation gave 
just three such periods. The second harmonic of the 6-year term 
gives a 3-year and the third harmonic a 2-year term, which was 
also derived, although no 2-year period was suggested by Angen- 
heister. The augmenting factors for the ampUtudes and the reduc- 
tion of the epochs to the beginning of the year 1900, both being made 
necessary by the use of 3-month groups as single observations, are 
easily found.. (See p. 79). 

From the harmonic constants of the 3 stations the corresponding 
motion of the pole was then deduced. Let Xj, Xg, and Xj, denote the 
resolved portion of the polar motion on the x axis (Meridian of Green- 
wich) due to the 6-vear, 3-year and 2-year terms, respectively; the 
subscripts denote the harmonic order of the term with the 6-year 
term as fimdamental; let y^, y2, and 1/3 denote the resolving portions 
of polar motion of the same periods on the y axis; for these terms the 
following values were deduced with which are shown Wanach's 
values for comparison. 



Period. 


This investigation. 


Wanach. 




/xi= Of 0031 cos (K^/e-SO") 


-0f0223cos(«</6-106"). ' 
-0r0083 COS (irf/6-268"). 

-Of0038 cos (i<</3-125'»). 
»0:0170cos(«</3-9"). 

-070084 cos (irf/2-49"). 
-0r0053 cos \Ktf2-r). 


6 years 


\yi-0f0121 cos (*</6-272*) 

r3<r-0f 0107 cos («</3-225') 




3years 


it/Qos 0^0146 cos ?«i/3— 41*") 


2 years 


ili-OfOU? cos (irf/2-204«) 




\y3— 070136 COS (irf/2— 284*) 



The /c in these equations is, of course, the same as in the expressions 
for the annual portion of the Kimura term on page 55. Wanach's 
values are from observations covering the years 1900-11, inclusive, 
with all 6 stations. The other values depend on 3 stations but cover the 
longer period 1 900-1 7, inclusive. There is some resemblance in the two 
sets of values for the 3-year and the 6-year terms that has survived the 
diflFerence in method of treatment, in the stations used, andin the period 
covered. Stace, furthermore, there appears to be a priori reason for 
expecting such terms, we may consider their existence as fairly prob- 
able, but how accurately the polar motion may be given by either 

^ "^ber die dreijahrige Luf tdruckschwankung und ihren Zusammenhang mit PolschwaDkungen; Nach- 
richten der K5mglichen Oesellschaft der Wissenschaften £U Gdttingen, Math.-ph]^. Klasse, 1914. p. 1. 
The paper is described as a preliminary cammunication, but nothing further on the subject from Angen- 
heister nas come to the author's attention. 
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set of expressions is another question. The 2-year terms seem to be 
merely tne result of manipulating figures. The 3-year terms in the 
expressions for the polar motion reproduce better the terms of like 
period at the several statioi^s from which the former were derived 
than do the terms in the polar motion with periods of 2 or 6 years. 
The mean residual (without regard to sign) in the s's and c's of the 3 
stations is 0f0058 for the 6-year terms, 0f0027 for the 3-year terms, 
and 050088 for the 2-year terms. 

The smnmation for the 6-year component afforded a partial check 
on the niunerical work. The sixth harmonic is, of course, the annual 
component. The annual term so obtained, when all corrections have 
been appHed, should be equal to the mean result of the three con- 
secutive 6-year series already obtained, namely, 1900-05, 1906-11, 
and 1912-17. In taking the mean of the three it is more accurate to 
average the c's and s's than the ampHtudes and epochs. The resiilts 
of this check were entirely satisfactory; the epochs agreed within a 
fraction of a degree and the amplitudes withm a few units of the 
fourth decimal oi a second, that being the last decimal place retained 
in the computations. 

The fifth harmonic of the 6-year component comes out very large; 
it is explained by the fact that the period of this harmomc, 1.20 
years, is so nearly that of the 14-montn term, or 1.184 years, so that 
the two periods dfo not separate much in the course of 18 years. He 
fifth harmonic can, in fact, be made to furnish an approxmxate value 
of the 14-month term, thus affording a rough check. 



Chapter V.— SUMMARY AND CONCLUSIONS. 

The rate of apparent increase in the latitude of Ukiah found by 
Prof. Lawson, 0f0094 per year, is reduced to 050081 per year by 
the appUcation of a more accurate method to the same material as 
Prof. Lawson used. Neither rate is to be further increased on 
account of any supposed change in the zero frdm which changes in 
latitude were reckoned. 

To make sure that the apparent increase is not due to the star 
places used, the other stations of the International Latitude Service 
were examined for possible changes of latitude similar to those at 
Ukiah. The definitive latitudes of the International Latitude 
Service from 1900.00 on were used so far as available and then the 
provisional latitudes through the year 1917. The details of the 
metliods used are given in Chapter II, pages 8 to 30. The results 
for all stations are given in Tables 8 and 9. These tables indicate that 
the increase in latitude at Ukiah is in no wav exceptional, that all 
stations except Mizusawa show, on the whole, apparent increases 
in latitude throughout the period covered, and tnat toward the 
latter part of that period Mizusawa also shows an increase. The 
average rate of increase of Gaithersburg usually exceeds the averagje 
rate of Ukiah over like periods of time. The rate of increase of lati- 
tude at Ukiah varies from 0f0063 to Of0260 a year accordihg to the 
method used and the period covered. The rate at Gaithersburg 
varies from Of 0083 to 0f0218 a year. The larger values just men- 
tioned for the two stations are too uncertain to be used as a basis 
for any trustworthy conclusion. The more rehable rates are not 
far from the rate of 0f0081 already found. (For particulars, see 
Tables 8 and 9.) 

A general apparent increase in latitude at all stations might mean 
an approach of these stations toward the pole, but a much more 
probable explanation is an increasing error in the star places due to 
erroneous values of the proper motions, an error that would affect 
all stations alike, since the star program is aUke at all stations. An 
inspection of Tables 8 and 9 shows a comparatively large increase at 
American stations and a small increase or an actual decrease at Asiatic 
and European stations, which are on the opposite side of the pole 
from America. This suggests that superposed on the purely apparent 
change of latitude due to star-places there may be a change of latitude 
due to a displacement of the North Pole toward the American con- 
tinent. A least-squares adjustment on this hypothesis of the rates 
of change given in Tables 8 and 9 gives different out fairly consistent 
results according to the data and weights used. These results are 
set forth in Table 10, page 34. An arbitrary mean of these results 
was adopted as one defimtive result of this mvestigation, namely a 
progressive shifting of the North Pole at the rate of 0.^0050 (equal 
to 15 centimeters or 6 inches) per year in magnitude and directed 
toward the Equator along the meridian of 78 west of Greenwich. 
This is a mean result for the years 1900-14, inclusive; the observa- 
tions may be interpreted as meaning that a polar motion of the same 
general character continues through 1917. A revision of the calcu- 
lations was undertaken while this pubhcation was m press and was 
completed too late for the results to be thoro^hly incorporated in it. 
The revision led to the complete rejection of Tscnardjui in the deter- 
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mination of the polar displacement. The effect of the rejection is 
to increase the annual rate of displacement of the North Pole for the 
period 1900-1917 to about 0f0062. The direction of the displace- 
ment is a few degrees nearer to the meridian of 90° west than the 78® 
just mentioned, the direction apparently varying somewhat during 
the period tested. For reference to the article containing the 
revised calculation see the last paragraph of the introduction, page 2. 

In regions situated on mendians nearlv at right angles to the 
meridian of 78° west the shifting of the pole would have practically 
no effect on the latitude. One such region is evidently western 
Europe, which is just the region where the longest series of accurate 
observations havft been made. These observations would then be of 
little use in decidii^ whether or not a polar displacement of the 
general nature of the one here deduced has been going on for some 
time in the past and might therefore be expected to continue for 
some time in the future. The shifting of the pole during the years 
in question may be due to a combination of periodic terms having 
periods of a few years or a few decades, terms whose existence has 
not yet been established but which might correspond to meteoro- 
logical cycles of one sort or another. Verv recent observations at 
UEiah might indicate that this shifting of the pole is not continuing 
and might even be reversing itself, but the evidence of Ukiah alone 
on this point is not conclusive. 

This shifting of the pole toward the American continent combined 
with a change in the average proper inotion of the stars used of 
about the same magnitude (see values of t^ in Table 10) appears to 
be a fairly satisfactory explanation of the apparent increase in lati- 
tude at Ukiah, and at all other stations of the international Latitude 
Service during the period in question, without invoking the effects 
of a creep of surface strata. No attempt is made to pronounce on 
Prof. Lawson's hypothesis in general; the question in this investi- 
gation is simply the interpretation of the astronomical evid^ace with 
special reference to Ukiah. The hvpo thesis of a surface creep at 
iJkiah is not absolutely excluded, but if a creep exists, its effect 
appears to be subordinate to the effects of polar motion and of errors 
in the assumed proper motions of the stars. A brief discussion is 
given of the geophysical questions involved in a possible progressive 
shifting of the jjole (pp. 39 to 42). 

The passage in Prof. Lawson's article dealing with the latitude of 
Lick Observatory is reproduced on page 44. A rate of change of 
latitude of .Lick Observatory before and after September, 1903, is 
deduced by least-squares adjustments from the data used by Prof. 
Lawson. Ilates of mcrease as found from these adjustments resem- 
ble those deduced by Prof. Lawson by a graphic process, but the 
value of the entire result is rendered doubtnil by a lack of homoge- 
neity of the declinations used. There is some evidence in the lati- 
tudes themselves of a sudden decrease in the latitude of Lick Obser- 
vatory in September, 1 903, about the time of a severe local earthquake. 
Some decrease of this general nature is accepted as fairly plausible 
by the director, Dr. Campbell, but is doubted by Mr. Tucker, the 
astronomer who made the latitude observations. Triangulation exe- 
cuted both before and after 1903 gives no clear evidence of a shift 
in the position of Lick Observatory. If this sudden decrease be 
accepted, there is then evidence, based solely on stfirs belonging to 
a single declination system, that the latitude of Lick Observatory may 
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have increased progressively at about the rate due to the motion of 
the pole given above. The lack of homogeneity of the declinations 
for any extended period renders any very definite conclusion on the 
subject rather difficult. 

The harmonic constants derived in the course of the work to repre- 
sent the variation of latitude at the various stations are corrected for 
aU known sources of error and the definitive results are given in 
Tables 11 and 14 (pp. 51 and 57). 

The portion of the variation of latitude having an annual period is 
composed of the motion of the pole with an annual period and to 
various other effects not due to the motion of the pole and to a cer- 
tain extent different from station to station. These other effects are 
lumped together as the annual part of the Kimura term. The annual 
portion of the motion of the pole, as deduced from an adjustment of 
the harmonic constants of the several stations, varies according to 
the way in which the Eamura term is treated. The results for the 
annual motion considered to be on the whole the most satisfactory 
are given in the two lower lines of Table 12, on page 53. They agree 
satisf actorilv with results deduced by other methods by Dyson and 
Wanach. Expressions for the annual portion of the Kimiu*a term itself 
are given bv equations (1), (3), (5), and (6), (pp. 55 and 56) ^ agreeing 
satisfactorily with Wanach's and Dyson's results, equations (2) and (4) . 

Harmonic constants for the semiannual term are given in Table 14 
(p. 67). Attention is called to an increase in the amplitude of this 
term toward the end of the period treated. It is suggested that some 
of the semiannual effect may be part of the Kimura term and may 
be due to systematic errors m the declinations having a period of 12 
hours of right ascension. The annual motion of the pole is shown in 
figures 14 and 15. 

The motion of the pole in its ''free" period of 432.5 days is deduced 
from the harmonic constants for this term at the several stations. 
There is no evidence from the observations of any effect of this period 
in the Kimura term and no a priori reason for an effect of this sort 
large enough to be considered. Harmonic constants for the motion 
of the pole in the freeperiod are given in Table 16 (p. 60). The 
results of Dyson and Wanach deduced by other methods are given 
for comparison; they agree satisfactorily as far as they go with the 
results m Table 16. 

The chief point of interest in Table 16 is the direction of the major 
axis of the eUipse of polar motion. The results vary considerably 
and are affected by large probable errors, but a direction not far from 
the meridian of 90^ west of Greenwich is pretty clearly indicated. 
This longitude is the longitude of the axis of the larger principal 
equatorial moment of inertia, if there is any perceptible difference in 
the principal equatorial moments. The lai^er principal moment 
correspon(t to the shorter axis of the geoid considered as an eUipsoid 
of three unequal axes. The eccentricity of the ellipse of polar 
motion affords a means for deriving the difference between the two 
equatorial moments of inertia or between the corresponding equa- 
terial radlii of the geoid. The eccentricity of Table 16 vary consider- 
ably but seem to indicate, on the whole, a difference between the two 
equatorial radii of the order of magnitude of 200 or 300 meters.^ 

- ■ .III I ■! Illllll. I Ill 

1 These results are subject to a correction not yet precisely evaluated but relatively rather small for the 
mobility of the ocean waters. (See pp. 99 and 101.) For this reason the results me stated only tenta- 
tively and in round numbess. 
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These results (radius of the equatorial ellipse in longitude 90® W., 
200 or 300 meters shorter than the radius m longitude 0° or 180°) 

Xe roughly with the results of Helmert from gravity observations 
rtest radius at 107° W., difference in radii of 230 meters). The 
determination of these quantities from the latitude observations is 
necessarily rough, and an accurate result would require a long series 
of observations, but even the rough result is of value as fending 
plausibility to Helmert's conclusions because the latter appear to con- 
flict with our ideas about conditions in the interior of the earth and 
to be perhaps open to suspicion on accoxmt of systematic topographic 
or geologic effects in the gravity observations. 

The observations were exammed with a view of detecting possible 
terms in the motion of the pole having periods of 3 or 6 years, there 
being some a priori reasons for expecting terms of this sort. The 
results are obtained from three stations over an 18-year interval. 
They are compared with results obtained by a different method bj 
Wanach from a 12-year interval but based on all 6 stations. The 
two results of the two discussions are compared on page 65. Some 
little resemblance between the two has survived the ainerence in the 
stations used, in the method, and in the period covered. This resem- 
blance, combined with a priori considerations, appears to be just about 
sufficient to render fairly probable the existence of terms oi this sort, 
but the values of their amplitudes and epochs are far from certain. 

The motion of the undisturbed pole of figure as deduced from the 
annual motion of thepole of rotation is given on pages 105 and 106, 
figures 14 and 15. Tne amplitude of the motion of the undisturbed 

f)ole is much smaller than that of the pole of rotation, and its motion 
ess accurately determined. 

It appears to the writer, as a result of this investigation, that it 
would DC desirable to repeat some of the calculations wnen the defini- 
tive latitudes become available for the years after 1912.000. 

It also appears unfortunate that the number of latitude stations 
should have Ibeen reduced to three and that one of these shoidd be in 
a re^on where the ground is supposed to be unstable. Three stations 
are just enough to determine the harmonic constants of the annual 
motion of the pole with the annual part of the Kimura term included, 
and additional stations to serve as a check are highly desirable. 

Three stations are just enough also to determine a progressive 
shift of the pole combined with a constant correction to tne average 
proper motion in declination (the w of Chap. II, p. 33). Again, addi- 
tional stations to serve as a check are imperative if the results are to 
be considered as really; certain. 

The comphcations introduced by the declinations in delicate 
researches on the motion of the pole are such that it is desirable to 
be independent of declinations if possible, even at the expense of the 
inaccuracy of daylight observations. It, therefore, appears desir- 
able to determine latitude from observations on circumpolar stars 
at both culminations. The observations should be made system- 
aticallv and as frequently as possible, and the results for each star 
at each observatory should be discussed separately. In this way it 
is not necessary to have a chain of latitude observatories all on the 
same parallel; any well-equipped and well-staffed observatory out- 
side of title Tropical Zone would be able to contribute results of value 
in the study of the variation of latitude. 



Chapter VL— MISCELLANEOUS MATHEMATICAL DEVELOP- 

MENTS. 

This chapter contains various mathematical developments regard- 
ing the details of the calcidations previously referred to. The mathe- 
matical notation of each section of the chapter is independent of 
that of the other sections. 

SECTION 1. SUMS OF CERTAIN TRIGONOMETMO SBBIBS— OLBABANCE 

FORMUIA8. 

The following formulas for the simmiation of a series of sines or 
cosines of angles in arithmetical precession are continuallv needed 
in the developments of the harmonic analysis. The proofs will be 
found m ahnost any work on analytical trigonometry. 

sina + sin(a+/3)+sin(tt + 2/S)-|-sin(a + 3iS). • • • -hsin[a+(7i-l)j8] 

(. ^~lio\ • ^i^ P /1\ 

cosa + cos(a-l-/S) +eos(a+2j8) +cos(a + 3/3) -fees [a+ (n- l)fi] 

= cos( aH — 2~"i8 jsm -^ cosec^- (2) 

When 90 equals 27r or a multiple of 2t the sum of the terms of the 
left-hand members of (1) and (2) is evidently zero. 

Let us multiply equation (1) by i, then in (1) replace a by a+ jS, 
multiply the equation after the replacement by i and add to the first 
result. Every term in e&ch series nas a term in the other exactly like 
it except the nrst term of the first series and the last one of the second 
series; by combining like terms we find 

|sina + sin (a + jS)+sin (a+2jS)+sin (a + SjS)- • • • -hsin [ac+(n — l)jS] 
-hj sin (a + njS) 

==isin(«+!Lzi^)sin^oosec§ 



+ 



, . / , n+l^\ . fi/3 j8 

J sm ( ac+— o— i^ ) sm —■ cosec ~f 



or by obvious transformations of the right-hand side 

i-sin a + sin (a+/3)-|-sin (a + 2/3) • • • • +sin [a+(n—l)/3] 
+ isin (a+njS) 

= sin (a+^fijsm — cot |- (3) 

In a similar way from (2) 

i cos a + cos (a + /3) 4-cos (a-f 2^) .... +cos [a+ (w— 1)0] 
+ i cos (ac + 9ijS) 

= cosf a+^i3jsinycot^- (4) 
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Formulas (3) and (4) apply to the case frequent in practice where 
a set of ordinates is to be summed with only naif weignt assigned to 
the ordinates at the beginning and end. 

On page 5 formula (4) is applied to an expression of the form 

A cos (irf — a) . Fifteen ordinates are used^ or 7i = 14. It is assumed 

that irf = 27r when < =432.5 days, this being taken as the length of the 

free period of the latitude variation. Then dividing this period into 14 

432.5 432 5 

parts, corresponding to ^ = 0, ^1 > 2 X ■. / > etc, days, we find 

\ A cos (0 — a) -\-A cos Tt^ — aW^ cos (2 Xj? — a ) 

H-^cosri3XY^-aj+i^cos (^14Xy7-a) = 0. 

Again we use the same value of k, but instead of extending the 
ordinates over 432.5 days we extend them over exactly 14 mondis or 
rather 14/12 of a Julian year, that is, 426.125 days. Then putting 

7==^||^X^=0.140751ir = 25?3353 
we find 

i^ cos (0 — a) +-4 cos (7 — a) + J. cos (2 7 — a) • • • • +^ cos (13.7 — a) 

-hi A cos (147— a) ^A cos (77 — a) sin 77 cot ^* 

The numerical value of sin 77 cot ^ comes out 0.206; therefore, since 
cos (77 — a) is never numerically greater than unity, the mean of the 
fifteen ordinates (end ordinates half weight) cannot exceed ' .. . — 

= 0.015 -4, as stated on page 5. 

The formulas on pages 14 and 17 for clearing the mean latitude 
over a giren period ^om the eflFect of one component are easUy 
derived from formula (4). Let us consider first the formula on 
page 14. 

Let us represent the 14-month component whose eflFect is to be 
removed from the mean over a year bv J? cos (k^— fO where f' 
is the f when the time is reckoned from the beginning of the year. 
If the time is reckoned in months k=25?33526, and i takes on the 
values 0, 1, 2, 3, . . . 11, 12. We have therefore in (4), 

a= -f, /3 = 25?33526, n=12, 

and the mean of the 13 terms (first and last half weight) is 

1^5 cos (6k -fO sin 61c cot \k, (5) 

The product ^ sin 61c cot ^it is independent of the harmonic con- 
stants of the 14-month component; its value is 0.174, of which the 
logarithm is 9.2405 — 10. To connect the f' with the f reduced to 
1900.00, the latter being indicated by the absence of a prime, we 
have for the change of phase in one year 12ic = 304.023. Therefore the 
phase falls behind in one year by 360°~304?023 = 55?977, which is 
the amount by which the f reckoned from the beginning of any year 
exceeds the f reckoned from the beginning of the previous year, or 

f'-f-h55?977n, (6) 
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whjere n is the number of years elapsed since 1900.00. Since 6 ic== 
152?0, expression (5) becomes 0.174 -B cos(152?0-'f-55?977n). 
This expression is contained in the mean as taken; to clear the mean 
it must be subtracted, or the angle may be increased by 180° and 
the resulting expression added. This process gives for the quantity 
to be added algebraically to clear the mean over a year from the 
effect of the 14-month component the expression: 

+ 0.174 fl cos (332?0-f-55?977 n), (7) 

A similar process is used for deriving the expression to clear the 
mean over a 14-month (432.5 dav) period from the effects of the 
annual component. The result, wnen 13 values enter into the mean, 
with half weight to the first and last, is 

+ 0.142 R cos (213?l-r + 360y). (8) 

In (8) , i? and f are the amplitude and epoch of the annual compo- 
nent ; / is the interval expressed as a fraction of a year from the begin- 
ning of the year to the beginning of the series* that is, to the tune 
of the first value of the set of wmch the mean is taken. 

SECTION 2. HARMONIC ANALYSIS — ELIMINATION OF THE EFFECT OF 

ONE COMPONENT ON ANOTHER 

It may be well at this point to call attention to the two-fold use of 
the word ''component.'' The context will show which use is meant. 
There is the ordinary use according to which we speak of the com- 
ponent of any vector quantity in a given direction and there is another 
use, very frequent in this report, borrowed from the theory of the 
harmonic analysis of the tides; in this sense a component is simply 
one of the periodic terms into which the mathematical expression 
for a periooic phenomenon may be separated. One component is 
distinguished from another by its period, and a component of given 
period may comprise all parts of the mathematical expression having 
that period or its submiutiples but not multiples of tne period.^ In 
practice, however, these periodic terms or components are limited 
to harmonic terms, and, therefore, are not separable into terms of 
shorter period. 

In this report, in conformity with the practice adopted in the har- 
monic analysis of tides, these components are assumed to be of the 
general form R cos (fli — f ) . The quantities 5, a, and f, are constants; 
t is the time. R is called the amplitude of the component, and is 
always taken positive; its physical dimensions are the same as those 
of the quantity represented. The quantity a is called the speed and 
is related to the period P of the component by the equation aP=27r 
or 360° according as radians or degrees are used in the computation. 
The speed of one component may be a simple multiple of another, or 
the two speeds may be practically incommensurable. The quantity f 
is called tne epoch and it evidently depends on the origin from which 
the time t is reckoned, f is the interval from the origin of time to the 
first maximum of the term R cos (a^— f), the interval being of course 
in degrees or radians of which there are respectively 360° or 27r to a 

^ For example, an annual ccanponent has a term with a period of 1 year and may have terms, ii the 
latter are not especially excluded, 'with periods of 6 montl^, 4 months, 3 montiis, etc. A term with a 
period of two or three years would not be a i^rt ol the annual compotient . 
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period. The quantity 360° — f represents the phase at the origin of 
time of the penodic variation represented. 

The quantities B and f are called the harmonic constants of the 
particular component in question. In this report all f's as finally 
given are reduced to the beginning of the year 1900 as origin of time, 
regardless of whether the observations from which they are deduced 
include the year 1900 or not; the application of this reduction to 
1900.00 assumes of course that the pnenomenon represented repeats 
itself indefinitely. The expression Bcosiat — ^) may be expanded 
iato RcobS^ cosat+Rsiafsvaat'=ccos(U+s sin at, where c=fe cosf 
and s^Rsinf. The quantities c and s are also sometimes called 
harmonic constants. 

The relations among these four quantities R, f , c, and «, may be 
represented graphically. Take rectangular axes OC^ and 08 (fig. 11). 
If € and s be taten as the rectangular coordinates of a point P, then 
R and t will be its polar coordinates, the pole being and the initial 




O ^ M 

FiQ. 11.— Belatioia among yarioas seta of hannonic constants. 
The r*a in the figure represent the R of the text. 

line OC. A change in the origin time changes the f by an amount 

f)roportional to the interval l^tween the oW and new origins, but 
eaves the R unchanged. If the quantities reckoned from a new 
origin of time be given primed letters to distinguish the corresponding 
quantities reckoned from the old origin and written without primes, 
we have for the difference between the «'s, 



We have also, 



^ = f-f' = a(^-r)=aT, 



(1) 



c' = 5 cos f' = ff cos (f-^) 



Similarly, 



^R cos f cos $-\-R sin f sin ^ 
= c cos S + s sin 6. 

«'«— csin^ + s cos 6, 



(2) 
(3) 



Other writers sometimes assume as standard forms for the har- 
monic terms, R cos {at + Tj), R sin {ot—rDf R sin (at+ri). The S's 
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are of course the same in all four fonns. Wanach uses the last- 
named form in the Resultate des Intemationalen Breitendienstes, 
volume V. The relation between the f used in this report and 
Wanach's ty, is 

f+t7 = 90°or450°. (4) 

Let us now suppose for simplicity that the effect to be analyzed 
has two components only; that is, we are given a number of values 
yo> Vv Vif Vsj etc., of the expression 



y^A cos (Tn^—fi) +B cos (r^— fj) 1 
= Cia cos mt+Sjo. sin mt+Cn cos 7U-{'8n sin nt. J 



(5) 



The speeds m and n are known. The object is to determine the 
values of the amplitudes A and B and the epoch's fi, and ^2; or else 
the c's and s% from which the former set of quantities may be 
deduced. In theory four values of y would suffice to determine the 
four quantities required, but in practical applications of the harmonic 
analysis the y's are so affected by accidental errors that a determina- 
tion from four values of y would be of little use, particularly if the 
corresponding times were near together. In practice either the values 
of y are given at imif orm intervius of time or else a curve is available 
from which may be read the value of y at any time. The latter case 
is the one treated here. 

To find Cm and «m (or A and f^) it is convenient to take a imiform 
interval between the values of y equal to some submultiple of the 
period of the first term. For denniteness we shall carry the reasoning 
through with the number most frequently used for the purpose in 
this report, namely, 12; that is, we take curve readings at intervals 
of one-twelfth of the period of the first term. For this interval 

m^= 360^12 = 30°. 

We denote these readings in order by y^j j/j, y^, etc., and arrange 
them in 12 groups as follows: 



Group number. 


Value of mi. 


Readings. 


Mean 
Talue. 







90»-fj(360«) 


Vo, Vis* yHrVM, . . . . 
yi, yn, y», yw, . . . . 

Iff, yi4, y«,yn, 

yt, yti, ya,y», 


Y. 

1 


(6) 


1 


2 




3 








11 


330«+i(3W) 


Vih yn, Vi5» ir«7» . . . . 


Yn 

* 









j is an integer equal to zero for the y^s in the first subcolumn in the 
column of readings, to unity for those in the next, to two for the 
next, and so on. It is easy to see that the values of y in each group, 
in so far as they depend only on the first term A cos (tm^— f), are all 
alike because the values of (mt — f ) differ by 360° and the term repeats 
itself. Although within a group the first term repeats, the second 
term does not do so; but if we suppose the group to cover a long 
period of time, the second term will be now positive, now negative, 
and its mean value wiU be nearly zero. How nearlv zero will be 
considered presently; for the present we simply neglect the mean 
value of the second term in the several groups^ and suppose the mean 
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value of the second term in the several groups to be what any one of 
its members would be if the first term were the only one present. 
Using the second form of (5) and considering for simplicity the mean 
value of a group to be that of its first member we get, 



YQ = Cm cos 0°-hSm siu 0®, 

Fi = Cm COS 30°-l-«m siu 30°, 
Fj = ^m COS 60° + 8m siu 60°, 

^11 = Cm COS 33d°+Sxa siu 330' 



(7) 



These are 12 equations for only two imknowns Cm and 5m; we treat 
them as observation equations in the method of least squares. The 
normal equation for Cm is evidently 

[cos^ (k 30°)] Cm + [cos (Tc 30°) sin (k 30°)] Sjn^lylc (cos Jc 30°)]. (8) 

Square braclcets indicate siumnation with respect to 1c, which takes 
on integral values from to 11, inclusive. Bjr writing cos^ (k 30°) = 
i -f i cos (k 60°) and sin (k 30°) cos (Jc 30°) = i sin (A: 60^) and applying 
formulu (2) on page 71 we see that the coefficient of Cm is simply 6 ana 
the coefficient of ^m is zero. We have then 

6cm = [rikcos(fc30°)]. (9) 

In a similar fashion the normal equation for Sm simjiifies into 

6sm = [i^*sin(A:30°)]. (10) 

where k as before is successively 0, 1, 2, .... 11. 
It may be convenient to suppose a constant term h added to the 
expression (5) for y. It will evidently affect all values of Y^ aUke in 
(6) and have no effect on the correctness of formulas (9) and (10), 
since [h cos (k 30°)] and [h sin (k 30°)] are clearly zero. Formulas (9) 
and (10) are the essence of the harmonic analysis and the forms for 
the analvsis are designed simply to furnish convenient forms for com- 
puting the indicated summations in (9) and (10) and in similar equa- 
tions for terms whose periods are submultiples of those of the nrst 
term and for deducing the R^8 and f *s from the ^'s and c's. 

In taking up the effect of the omitted second term it is convenient 
to make a slignt change in the form of equations (9) and (10). Con- 
sider the different values of the v^s in Yq. According to formula (9) 
they are all multiplied by cos ; it is more convenient, however, to 
consider t/o ^ being multiphed by cos 0°, 2/12 by cos 360°, ^24 by cos 
2 X 360°, etc. Sinmarly for the multipUers of the y's in Y^ we may 
use successively cos 30°, cos (30° + 360°), cos (30° + 2x360°) etc. 
Let i represent the subscript of any y in the array (6) , and suppose the 
curvo-readings used cover p periods, so that the greatest value of i is 
12 p— 1. Then it is easy to see that (9) and (10) may be written 

6cm==- IViCos (ix30°)], 1 

? (11)' 

6 5m = -[yiSin (ix80°)]. 

The values of i in the summations range from to 12 p— 1. Equa- 
tions (11) do not give the true value of Cm and «m but the values of 
certain approximations to them; let us denote these approximations 
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as computed from (11) by Cm! and s^ and reserve the symbols without 
tie pnme for the true values. From the second fonn of (5), by 

writing X for — 30° 

Vi^c^ cos (i 30°) 4-«m sin {i 30°) +Cn cos (ix) cos {i 30°) 

+ 8n sin (ix) , 
yi cos (i 30°) =c„j cos^ {i 30°) +«m sin (i 30°)cos (i 30°) 

+ Cn COS (ir) cos (i 30°) + «n sin (ix) cos (i 30°) (12) 

= iCm + iCmCOS (t 60°) + i «m siu (t 60°) + i Cn COS t (30° + aj) 

+ icncos'i(30°-aj)+i5nsini(30°+a;)-i«nsini(30°~a;). (13) 

By performing the summations indicated in (11) using formulas (1) 
and (2), page 71, we see that the terms in sin {i 60°) and cos (i 60°) 
give a zero result in the sum and that we get 

6 c„-==6 0.+ ^[,„ cos [(6 p-1) (30°+.)] tin?(30°+tf ^ 

Vta 1\ fono si mi [6 p (30° - X)] 

+ CnCOsL(6p-l) (30 -X)j 3i^/(30O_^)^ 

+.sin[(e,-I)(30o..)] tjyff°/^f J 
-.sin[(6p-l)(30O-.)]«4^X^^}. a4) 

For brevity let us call the andes (6p- l)(30°+a;) and (6p- 1) (30°-a;), 
6 and <t>, respectively, and tne factors 

sm[6p(30°+x )] sin[6p(30°~x) ]^ respectivelv 

sin i (30° +x) ^^^ sin i (30° ^x) ^ ^^^-'^ respectively. 

Then (14) may be written 

6cm' = 6cm + 2^ Cn {F COS 6 -^f COS <p) -fSn (FamS-fsin (p) . (15) 
A similar calculation gives 

Cn (Fsin ^+/sin <p) +«n (/cos ip— F cos 6) . 



QStr/'^QSm-^ 



2p 



(16) 



The quantities within the brackets in (15) and (16) are always finite 
and as p increases we have more and more nearly 



Cm. O Cm, 

and 

Sja ^^ O Sail 



(17) 



as was assumed without strict proof at the beginning of this section. 
Even if jp is small, equations (17) may still give a good approximation. 
If the time covered by the I2p readings, which by hypothesis corre- 
spond just p p^ods of the first term of (5), also corresponds approx- 
hnately to an integral number of periods of the second term then the 
angles </> and B are nearly 0° or 180° and their sines are nearly zero. 
The quantity -P is in the cases that arise in practice much smaller 
than^ the factor sin i (30° +x) in the denominator being larger than 
sin i (30°— x). In the tidal work the series are sufficiently long so 
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that terms in F are dropped entirely. With sin <^ = and terms in 
F dropped the approximation (17) still holds good. 

Even when these conditions do not hold, equations (15) and (16) 
are still useful. Let us evaluate known portions for a case needed in 
the variation of latitude, namely, when the first term of (5) repre- 
sents the annual component and the second term the 14-month com- 
ponent and when p = 7. The value of x is 26?33526. We find 

^=136?41, 
</>«166?41, 
F= 0.591, 
/== 6.740. 

If the 14-month component had a period of precisely 14 months 
(14/12 years) instead of 432.5 days, the 7 years would represent 
exactly 6 periods of 14 months, the value of <t> would have been 180° 
and terms in sin <^, instead of being merely small, would have been 
zero. With these values using the subscript a instead of m for the 
annual component and the subscript h instead of n for the 14-month 
component, we get from (15) and (16) 

6ca + 0.4374 Cb~ 0.0840 Sb = 6ca', \ n9\ 

6sa + 0.1422 Cb + 0.4985 Sb = 65a'. J ^ ^ 

The right-hand sides of (18) represent the results of the ordinary 
uncorrected means of the groups shown in array (6) . 

Let Cm and 5m now refer to the 14-month component and Cn and 
Sn to the annual, and let us suppose curve-readings taken at inter- 
vals of one-twelith of the 14-month period, the first reading, v©, 
coinciding with the first reading for the annual component. If we make 
2) = 6, we are covering an interval of 6x432.5 days, or a little over 7 
years. We find in this way, remembering that we must use h for m 
and a for n in the subscripts of (15) and (16) 



6Cb + 0.5852 Ca + 0.1429 5a = 6Cb', 
65b -0.1668 Ca + 0.4886 Sa = 65b'. 



(19) 



The right-hand sides of (19) are known from the uncorrected results 
of the usual process of forming groups and means and analyzing them 
for the 14-month component. Equations (18) and (19) give four lin- 
ear equations for deducing the true values of Ca, Cb, 5a, and Sb from their 
uncorrected values cj, Cb , 5a', and 5b'. The solution of them gives 



Ca=-h0.1679(6CaO+0.0001(65aO-0.0122(6CbO+0.0023(65bO, 
5a=+0.0000(6CaO + 0.1679(65aO-0.0040(6CbO-0.0139(65bO, 
Cb=-0.0164(6CaO-0.0040(65aO + 0.1680(6CbO+0.0001(65bO, 
5b- -f0.0047(6CaO~0.0137(65aO-0.0000(6CbO+0.1679(65bO. J 



ym 



These are equations used in clearing one component from the effect of 
another for a so-called 7-vear series, which means that between the 
first and the last reading for the annual component there are 6 years, 
11 months (6.917 years) and between the first and the last reading for 
the 14-month component 7.005 years. The results are written in 
terms of 6ca', 6cb',etc., rather than in terms of plain Ca', Cb', because 
the former quantities appear on the forms for harmonic analysis, 
and the form adopted saves a division by 6. 

The approximate eauality obtained m this way between the 
periods of time covered in the two sets of readings is desirable in 
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view of the fact that the so-called harmonic constants are, in fact, 
somewhat variable and that the values obtained for them should be 
considered merely as mean values for the periods in question. An 
equality holding to about the same degree of approximation as the 
one just considered is 6 years = five 14-month periods. The interval 
between the first and last curve readings is 5 years 11 months (5.917 
years) for the annual component and 4^X432.5 days = 5.821 years 
for the 14-month component. This is the so-called G-yefu* series. 
The equations for deducing the corrected values of Ca, Cb, etc., from 
their imcorrected values are found for the 6-year series to be 

Ca= +0.1676(6Ca') -0.0001(6«a0 +0.0106(6CbO +0.0032(6«b0, 1 
«a= +0.0000(6ca0 +0.1677(6«aO -0.0027(6CbO +0.0126(60, l ro^\ 
Cb= +0.0146(6ca0 -0.0051(60 +0.1677 (6CbO -0.0001 (6sb0, | ^ ^ 
«b= +0.0033(6ca0 +0.0128(6«a0 +0.0000(6CbO +0.1677(6^0. J 

Equations (20) and (21) are the formulas used in passing from 
columns 3 and 4 of Table 3, pa^e 10, to columns 5 ana 6. 

In ordinary tidal work witn snort-period components the process 
of eUminating effects of components, other than the one analyzed, 
corresponds approximately to setting F= and using on the lef t-nand 
side of (18) the imcorrected values of Cb and Sh instead of the true 
values Cb and Sb- This process is not sufficiently accurate in the 
discussion of a year's observations on the long-period components 
and a process is used due to Darwin^ that resembles the one 
here described down to the derivation of equations (18) and (19). 
The equations, however, contain 10 unknowns, the s^s and c's of the 5 
long-period tides. It is suggested that these equations be solved 
by successive approximations, since one coefficient in each is much 
larger than the others. It woidd be entirely feasible to give a gen- 
eral solution of these equations in terms oi the known quantities, 
which are analogous to 6Ca', 6«a', 6cb', 6«b', as was done for equations 
(20). The amoimt of work in calculating such a general solution 
would be considerable, which is perhaps the reason wty Darwin does 
not give such a solution. However^ tne amoimt of work in solving 
the equations by successive approximation is considerable also, and 
has to be repeated for each year of observation discussed, while with 
a general solution available the computation is relatively brief. 

In the analysis to obtain possible 6-year, 3-year, and 2-year terms 
the 6-year period was divided into 24 parts instead of 12. The 
corresponding changes in formulas (9) ana (10) are easily made. 

In treating the 6-year component the mean of three readings a 
month apart was taken and the result used as a single reading. The 
necessary correction may be deduced as follows : For any component 
let at^ be the value of at for the middle reading. Let readmgs be 
taken at intervals t before and after the time t^. Then the three 
values of the term R cos (at — f ) are : 

B cos [a(^i-T)~f], B cos [a^i-f], B cos [a(^+T)-fl 

By combining the first and third the mean may be written 



» r -. o ri+2 COSiZr'l 

B COS [at^ - f] g h 



* The Harmonic Analysis of Tidal Observations, sec. 10, Sdentifio Papers, Vol. I, of British Associa- 
tion for the Advancement of Science, report for 1883. 
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which differs from the true quantity desired, namely, the middle 

reading, by the factor ^ +^ ^ °^ which for smaU values of r is 

nearly unity. This factor reduces all the values of the y's that 
enter^into tie harmonic analysis in the same ratio; the epoch r is 
therefore unaffected, but the amplitude Bis reduced from its true value 
in the ratio (1+2 cos ar):3. The amplitude deduced from compu- 
tations where a mean of three readmgs is treated as a single readmg 

3 
must be multiplied by yto '^^^ quantity may be called an 

X "T~ ^ COS QTo 

augmenting factor f which is the name applied in the harmonic analysis 
of tides where factors of this same general nature are used. With a 
6-year component and a month's interval aT = 5®. For a 3-year 
component ar = 10°, and for a 2-year component ar = 15°. The 
common logarithms of the augmentmg factors for the 6-year,3-year, 
and 2-year components are, respectively, 0.0011, .0.0044, and 0.0100. 

SECTION 3. THE CORRECTION FOR '' SLOPE'' IN THE HARMONIC 

ANALYSIS. . 

The expression for a function by a Fourier series takes the form 

/ W ==iCoH-Ci cos x + Cj cos 2a: • • • • l^x 

-\-8i sinx + Sj sin 2x .... j^ ^ 

1 r^' 

where Cn^"- I /(x) cos nxdx, 

T^J o 
1 P' 

and «n=- I /(x)sin nxdx. 

No provision is made in (1) for a term Jcx proportional to x. Such a 
term, if present, would mean that the sine curves of the various 
harmonics, instead of being laid off upward or downward from the 
level base line y=i Co, according to the usual graphic representation, 
would be laid off from an inclined base hne havmg the slope 1c. In 
discussing the variation of latitude there was found, besides the 
expected periodic terms, an apparent progressive change of latitude, 
roughly proportional to the time ana represented graphically by a 
line inclined to the axis of abscissas (the time-axisj; that is, there 
was foimd something analagous to the term tx, which we have sup- 
posed to be introduced into (1). 

A great variety of functions in no way periodic can be represented 
by a Fourier series within the range from x = to x = 27r, but if we 
suppose the function represented to be in fact essentially periodic 
except for a term Jcx and if we apply the ordinary formulas (2) to 
determine the coefficients in the expansion of such a fimction, neglect- 
ing the presence of the term ix, then the coefficients of the periodic 
terms will be falsified by this neglect. These coefficients, in addition 
to their proper values, will contain quantities which may be found 
for any harmonic term of order n by putting lex for/ (x) in (2) . Denot- 
ing by 7n and <rn the additional quantities introduced into Cn and 
«n we have 



(2) 



7ii=- I X cos nxcZx = 0, 



JeH' . . 2Je f 

(rn=- I xsm nxdx= • 



(3) 
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These are the quantities abeady in the coefficients as determined 
by the ordinary process; therefore, by reversing the signs we have as 

corrections to Cn and «n the quantities +0 and H — i respectively. 

Hf(x) is given only by its numerical values at uniform intervals 
from to 27r, the coefficients are then deduced by the methods of 
the harmonic analysis, the integrations in (2) being replaced by 
finite summations. It can be shown without much difficulty by the 
methods of the calculus of Finite DiflFerences, or simply by using the 
numerical values of Jcx in their appropriate places in any of the usual 
forms for harmonic analysis and by working the form through, that 
when the period is divided into 12 parts we have in place of (3) 

tJc 

Tfecot( 15^Xn) f(^) 

an g- 

For a period divided into 24 parts 

Tfccot(7i«'Xn) H5) 



I' 



7n- ^2' 



<^ii = 



12 



It is easily seen that for small values of n the value of a^ from (4) 
or (5) is not very diflFerent from that given by (3). The difference 
in the corresponding values of y is somewhat larger. It is explained 
by the well-known fact that at points of discontinuity of /(x), say 
for a; = Xi, the Fourier series obtained from (2) gives the mean of the 
two values of /(x,) or more precisely gives 

lim 



%i[^f(x,-^)+Ax, + €)2 



Such discontinuities really occur at a; = and x = 2w. Evidently, as 
far as the series (1) goes, /(0)=/(2ir), but the value of lex is zero 
iufone case and 2ir m the other. Hie series (1) therefore gives 
i (0 + 2ir)i=irt in both cases. The values of 7n in (4) and (6) were 
obtained with the first ordinate, or /"(O), taken equal to zero; if 
instead we take the first ordinate as tk and work through the form 
for harmonic analysis, we get 7n=0, whatever be the nunmer of parts 
into which the period is divided. Thi3 is merely to explain the differ- 
ence between the values of 7n as obtained from integration in (2) 
and from summations in (4) or (5). In the actual practice of cor- 
recting the results of harmonic analysis for a constant slope of the 
base Ime the Quantity 7n would have a value as given by (4) or (5), 
or whatever tne value might be for the number of ormnates used, 
and would not be zero. 

It is convenient to use instead of the slope Jc the difference between 
the ordinates of the inclmed base at the beginmnff and the end of a 

{>eriod; denote this quantitv by m. Then we nave m = 27rfe and 
rom (4) and (5) we get the following table of corrections. For com- 
parison, the values from equation (3) are also given; the integration 
m (3) corresponds to an innnite number of ordmates. 
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Correction to Fourier coefficients for inclination of hose line, 
[The quantity m is to be taken positive if the base line rises as x increases, and negative if the base line falls.] 



Ckirrection 
to- 



*i. 

«3. 
9i. 



12 
ordinates. 


24 
ordinates. 


Infinite 
number of 
ordmates. 


4-0. 311m 
4- .144m 
4- .083m 
+ .018m 
4- .022m 


4-0. 316m 
4- .155m 
4- .101m 
4- •072m 
4- .054m 
+ .042m 


+0.318m 
4- .159m 
4- .106m 
+ .080m 
4- .061m 
+ .053m 





Correction 
to— 



Cl 

C: 
Cn 
Ci 
C5 



12 
ordinates. 


24 
ordinates. 


+0.083m 
+ .083m 
+ .083m 
4- .083m 
+ .083m 
+ .083m 


4- 0.012m 
+ .012m 
+ .012m 
+ .012m 
+ .012m 
+ .012m 



Infinite 
number of 
ordinates. 












SECTION 4. FORMULAS FOR THE ELEMENTS DEFINING THE MOTION 

OF THE POLE. 

Throughout this section, except in the three concluding para^aphs 
on page 93, it will be assumed that the variation of latitude is Sue 
solely to the motion of the pole of rotation and furthermore that 
this motion is expressible in terms of harmonic constants of the ordi- 
nary sort. How far these assumptions are justified has been briefly 
considered in Chapter IV. 

The assumed harmonic motion of the North Pole is referred to 
coordinate axes passing through a mean position of the pole of rota- 
tion. In conformity with the practice oi the International Latitude 
Service, the positive direction of the x axis is taken southward along 
the meridian of Greenwich and the positive direction of the y axis 
southward along the meridian of 90 west of Greenwich. The har- 
monic motion for any given period is specified by x and y, the com- 
ponents along the corresponding axes; these components are taken 
m the form 

x = acos (Kt — a), \ (i^ 

y^h cos {d-fi). J ^^^ 

The time, t, is reckoned from some convenient origin; the quantity 
K is called the speed and is related to the period, P, of the harmonic 
oscillation by the relation 

360° 

or 

according as angles are reckoned in degrees or radians. The angles 
a and j8 are the mtervals from the assumed origin of time to the first 
maxima of the respective component oscillations. These intervals 
are reckoned in parts of the period, which is then considered as 360^ 
or as 27r according as a and j3 are in degrees or radians. 

The equations (1) considered as parametric equations in terms of 
t define tne curve in which the pole moves only as far as oscillations 
of the period in question are concerned. As is well known, the curve 
is an ellipse obUque, in general, to the coordinate axes. The equa- 
tion of this eUipse and formulas for the directions of its principal axes 
and the times when the pole is at the end of such axes, are derived 
without difficulty by the ordinary processes of analytic geometry. 
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The formulas that present themselves in this way are not well adapted 
to logarithmic calculation and other formulas will be derived that 
seem rather more convenient in this respect.' 

The distance of the pole from the origin at anytime is A<^ = V^^ +3/' 
and the maximum numerical value of A<^ corresponds to a position 
at the ends of the major axis of the ellipse, the minimum numerical 
value to a position.at the ends of the minor axis. From (1) 

A<t>^ = aj* +t/2 = (a^ cos* a + 6 cos^ jS) cos^ kI -f (a* sin^ a + &* sin* fi) sin' kI 

4- 2 (a* cos a sin a + 6* cos j8 sin p) cos kI sin kL 

By expressing functions erf ic^ in terms of 2Kt with a corresponding 
transformation for a and P we get 

A(^ = }(a» + &*) +i(«^ cos 2a + ¥ cos 2p) cos 2irf4-J(a* sin 2a 

-f- 6* sin 2jS) sin 2x1. 

This may be written in the form 

A0* = J (a* + 52) -h i c* cos 2 (kI -m) (2) 

by assuming 

a* cos 2ac + 6* cos 2/3 =c* cos 2/x 1 

(3) 
a* sin 2a -h 6* sin 2iS =c* sin 2ju. J 

When Kt=n OT Kt= n + lSO^ the value of A0* is evidently a maximum; 
when ic^=M±90° the value of A<^* is a minimum. Since n is deter- 
mined from 2/x, both m and n± 180® are included in the values of /x 
less than 360° that are given by (3). 

To adapt (3) to logarithmic calculation let us form values of 
c'cos [2m — (a -\- P)] and c *sin[2At — (a + jS)]. By expanding cos 
[2/i— (a+jS)] and using (3) there results, 

c* cos [2m -(a+jS)] = a* cos [2a -(a+j8)] + 6*cos[2/S -(a+j8)] 

= (a* + 62) cos (a-/3). (4) 

In a similar manner 

c^ sin [2m - (a +i3)] = (a' -6^) sin (a -j8) . (5) 

By dividing (5) by (4) 

tan [2m -(a +/3)] ^^ripja tan (a -/3) . (6) 

Assume the following equations for determining the auxiUary 

quantities d and 7 : j 

^ acos7=a, 

d sin 7 = bj 

where y may be taken in the first quadrant since a and h are positive. 
Then (6) may be written 

tan [2m - (a +/3)] = cos 2 7 tan (a -/3) . (8) 

Equation (8) taken alone determines four possible values of m less 
than 360*\ not only the m and the m± 180° foimd from (3) but also 
M±90°. To exclude the latter we note that by (4) 2M--(a+i5) 
must be so taken that its cosine has the same sign as cos {a—fi). 

s These formulas were proposed by the author as a problem in the American Mathematical Monthly for 
January, 1921, p. 36. 



(7) 
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By this criterion we exclude the values of m corresponding to /x ± 90® 
in (3) and get only the values of m corresponding to maximum values 
of Ail>\ 

If we denote by 6 the angle which the resultant displacement 
makes with the x axis we have 

A<^ cos ^*=x, 

A<^ sin d=y, 
or 

A<<>*sm 2^=2x1/, 1 

(9) 

A^*cos2d=a:*-^. J 

If in the expressions (1) for x and y we put ici=/x and use the result- 
ing values of x and y in (9)^ we get values of A<f>^ and 2$ correspond- 
ing to maximum values of A<^; that is, to the ends of the major 
axes. In this way we find 

2a^»2a& cos (kI —a) cos (ict — j8) »a&{cos [2ict — (a — j8)]+C08 (a — j8)}, 
or from (4) on putting ic^=m, 

Ail>* sin 2^= --F (a' + 6' + ^) cos (a -j8) . (10) 

c 

Again we have for the second equation of (9) 

7?-f=a^ cos^ {Kt-a)-¥ cosM«c^-j8)=i («'-&*)+ h (<3^* cos 2a 
-62 cos 2j8)cos 2/c«+i(ct' sin 2a -6* sin 2/3) sin 2icf, 

or by putting k^^/x and using (4), 

o/, 1 / , x^xj 1 a* cos* 2a -&* cos^ 2j8 
A</»* cos 2^ = 2 (^ -* ) + 2 c^ 

l a^sin* 2a -¥ sin* 2/ 3 
■^2 c* 



From (10) and (11) 



tan 2^=^33j3 cos va -/5). (12) 

The use of the auxiliary y again simpUfies the calculation and (12) 
may be written ^ ^fl^tan 2y cos (a -/3). (13) 

It will be foimd that (13) by itself gives four values of B less than 
360°, two corresponding to the ends of the major axis and two to 
the minor axis. The criterion for distinguishing those for the 
major axis is found in (10). Evidently, since oft, ^, a' + ft' + c* and 
A<^^ are all positive, for a maximiun, 20 must be taken so that 
sin 2d has the same sign as cos (a — /3). 
The value of A0 from (2) when id^ii gives the semi-major axis of 

the ellipse, or A* ~^^=W' +¥+<?). (14) 
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From (3) by squaring and adding 

c*=(i^ + J* + 2a2j3 cos 2{a-B) 
= (a^ -^h^y- 4a^l^ sin' (a - p) 

= (a'4-D')»[l-(^ sin («-«)']. (16) 

By using again the auxiliary angle 7, (16) may be written 

c* = (a^ + h^y [1 - sin» 27 sin» (a- jS)]. (16) 

We may introduce the further auxiliary S defined by 

sin 8= I sin 27 sin (a— jS) | . 

Taking S in the first ouadrant (16) becomes 

c^=(Pcos8, (17) 

smce by (7) d^=a' + 5». 

Then (14) becomes a a^ j t * /^ o\ 

A^ A<^x =» d cos ^5. (18) 

For the semi-minor axis B^A4>min it is easy to prove 

B = d sin ii. (19) 

The eccentricity of the ellipse (E) is found from 

A'-B' _ cosg^ . 

The values of A, B, and 6 may also be found by eliminating t from 
equations (1), which gives an equation of the second degree, which 
is easily shown to represent an eUipse and to which the standard 
methods may be applied for obtaining its dimensions and position; 
that is, the quantities A, B, and 6. It is interesting to note that the 
relations between 27, 2d, (a— j8), 5, and 2/i— (a+j8) are the same as if 
these quantities formed the parts of a right spherical triangle, the 
hypothenuse being 27 and the sides inclosing the right angle, 20 and 
d, the angles opposite these sides being, respectively, 90°--2M + «+i3 
and a— j8. Tms fact may be used to establish additional relations 
among the above quantities. 

The direction in which the pole rotates about the origin may also 
be deduced from the values 01 a and j8. - 

If 180°<a— j8<360°, the rotation is positive for the axes used, i. e., 
from east to west. Ii 0<a--j8<180°, the rotation is negative. 

Recapitulating the formulas and criteria we have: Motion of pole 

° ^ x^a cos (id —a) for motion along x axis, 1 /jx 

y=h cos (id — p) for motion along y axis, J 

a and b always positive. 

Compute the auxiliaries d and 7 by 

dcosy^a^ 1 /jjx 

dsmy^b, J 

d being always positive. 
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The time of maximum displacement (position at end of major axis) 
occurs when irf — Mj M being given by 



tan [2/i— (a+j8)]==cos 27 tan (a— j8). 



(Ill) 



Take 2/i— (a +j8) in such a quadrant that cos[2m— (a+j8)] has the 
same sign as cos (a— j8). 

The an^le which the direction of maximum displacement makes 
with X axis, or ^, is given by 



tan 2^= tan 27 cos (a— jS). 



(IV) 



Take 2B in such a quadrant that sin 2B has the same sign as 

COs(a — j8). 



Compute the auxiliary angle h by 

sin 5= I sin 27 sin (a— j8) | . 



(V) 



Take 5 in the first quadrant regardless of the algebraic sign of the 
ri^t-hand side of (V) . 

The semi-major and semi-minor axes A and B of the ellipse and 
its eccentricity E are given by 



A^d cos ^5, 
B=^d sin \h^ 

^ cos' \h 



(VI) 



The point whose motion is given by (I) describes the ellipse in a posi- 
tive direction (from the positive direction of the x axis towards the 
positive direction of the y axis) if 

or in a negative direction if 

0<a-j8<180^ 

There will be two values of m and two of d from (III) and (IV) 
each differing from the other member of the pair by 180^, even after 
the criteria nave been applied. Let us take for definiteness that 
value of B which lies on tne positive side of the x axis and affect all 

Juantities corresponding to it by a subscript zero, and the times and 
irections correspondii^ to the other maxima and the minima 
by the subscripts 1, 2, 3, according to the order of time in which 
they occur. Then the value of /up corresponding to ^o is chosen by 
the test that cos (mq— i^) is positive. We then have the following 
correspondences between the m's and the ^'s as they occur in order 
of time, according as the direction of rotation is positive or negative. 



Positive rotation 
180* <a-/5< 360* 


Negative rotation a 
0<a-/5<lW 


iti^ B^ Max. 
M1-P0+ 90* tfi-M- 90* Min. 
Mr-po+180* «i-0b±18O* Max. 
Mr-Po+270* «i-M-270* Min. 


/!« 9% Max. 
M1-/I0+ 90* tfi-fl^ 90* Min. 
Mf-^+180* «i-db±180* Max. 
M-/I0+270* ^-M- 90* Min. 



aThe direction in which the pole of rotation actually moves about the origin is negative for the axes here 
used. The formulas were originally devised for computations with the harmonic constants of tidal cor* 
rents and for these the rotation may have either direction. 
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The most convenient method of obtaining the elements of the 
ellipse from the harmonic constants of the latitude variation at the 
individual stations seems to be the indirect one of finding first the 
adjusted harmonic constants of the variation of latitude along the 
coordinate axes, that is, by finding the quantities a, 6, a and fi of 
equations (1). 

The observations at a station give data for obtaining that part of 

the variation of latitude with the period P= in the fonn 

A<I>==R cos (k^— f) 

= R cos f cos Kt+R sin f sin d. (21) 

But we also have by resolving the values of x and y along the meridan 
of the station Gongitude X, west longitudes positive) 

A<l)=x cos \+y sin X. (22) 

The X and y are those of equations (1) and on substituting their 
values in (18) we get 

A<^ = a cos (Kt—a) cos X + & cos-^k^— j8) sin X, 

= (a cos a cos X + J cos P sin X) cos Kt+ (a sin a cos X 

+ & sin j3 sin X) sin d, (23) 

Let us use the abbreviations 



a cos a=m, 
a sin a =71, 
b cos P—p, 
b sin P = q. 



J 



(24) 



By substituting these expressions in (23) and comparing with the 
second form of (21) we get for any one station, 

m cos \+p sin X = fi cos f =c, 1 /gc) 

n cos X + g sin X = -B sin f = 5. i 

Two stations at which R and f have been found give four equations, 
which are just sufficient to determine m, n, p, and a. From these 
the Quantities a, 6, a, and p are found by (24) and the elements of 
the eUipse by equations (I) to (VI). 

To treat the case arising whwi the harmonic constants of more 
than two stations are available some criterion must be foimd for 
testing the accuracy with which the adopted elements of the eUiptic 
motion represent the harmonic constants at the several stations. 
If we use the letters m, n, p, and g, to signify the adopted values of 
these quantities and denote by ct the error at any tmie in the re- 
presentation of the latitude oi a station due to the difference be- 
tween its harmonic constants as f oimd directly and as deduced from 
m, Uy Pj and q, then ct is given by the formula, 

€t= (m cos \-\-p sin X — c) cos /rf+ {n cos X + jsin X—s) sin d. (26) 

The error ct assumes in the course of a period (as d varies from 
to 2ir) all values between— €m and + €m, where €m is the maximum 
niunerical value of ct and is given by the equation 

€ni*= (m cos X+p sin X— c)*+ (n cos X+j sin X~«)^ (27) 
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The frequency of the errors of €t of a given size at a particular 
station does not vary at aU according to the Gaussian law, lar^e 
errors near the limits ±Cm being relatively more frequent than the 
smaller ones, but we may assume with some degree of plausibility 
that when the values of the €m's at a number of different stations 
are compared, the relative frequency of Cm's of different sizes will 
vary according to the Gaussian law, except of course as regards 
sign, the e^s being essentially positive. The mean square of €t 
or Ca' may be foimd by taking the mean value of €t' over a period, 
or 



t being in radians and the value of k being taken accordingly. It 
is easy to show that 

€a'==i€m^ (29) 

In obtaining the adjusted values of the elements of the polar motion 
thafHt most nearly to the harmonic constants of the several stations, 
the criterion used will be 

[6a'] = minimum. 

The square brackets indicate summation, as is usual in the theory 
of least squares, and the summation extends over all stations used 
in the adjustment. From (27) and (29) the quantity S to be mini- 
mized may be written 

S' = M = i{[(w^cosX+2>sinX-c)2+(7lcosX + gsinX--s)2]} (30) 

In order to have a minimum value m, n, p, and q must be so taken 
as to satisfy the equations 

dm 



dp 

dn 
dS 



= 0, 



or, 



dq 

m *" 



= 0. 



cos' X] + p [sin X cos X] = [c cos X], T 

m [sin X cos X] + p [sin' X] = [c sin X], Lr?i ^ 

n [cos' X] + g [sin X cos X] = [« cos X], rvoi; 

n [sin X cos X] + 2 t^in' X] = [« sin X]. J 

Square brackets indicate summation extending over all the stations 
used. 

Equations (31) would be the normal equations that would arise 
if we set out from observation equations oi the form 

m cos X+p sin X— c = 0, 1 ,09% 

n cos X + g sm X-« = 0, J ^^^^ 
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but these equations are not observation equations in the strict sense 
of the term. If we call the residuals of equations (32), v^ and V2 
respectively, that is, if v^ and V2 are the results of substituting in the 
left-hand sides of (32) the values of m, n, p, and g, deduced from 
(32), then v^ and Va have no particular significance separately, as 
their values depend on the epoch from which t is reckoned, but the 
sum of their squares is significant, for evidently 

Vi^4-V3' = €m' = 2€a^ (33) 

In estimating the probable errors of m, n, p, and q, as derived 
from (32), or flie probable error of a function of these quantities, 
the following conventions wiU be adopted, which are a natural 
extension of the convention by which the minimum value of M was 
taken as the criterion for obtaining values of m, n, p, and q. The 
probable errors of m, n, p, and q are computed just as if equations 
(32) were proper observation equations, but the mean square error 
of a single observation is not made to depend on v^ alone or V2 alone, 
but the following equation is used : 

(mean-square error)'= ^[^= ^^£3^' (^^) 

n being the number of ''observation equations,'' which equals twice 
the number of stations; the number 01 unknowns is of course four. 
The probable error, the Gaussian law being assumed, is the mean- 
square error multiplied by 0.6745. 

The difference between the harmonic constants at a station obtained 
from the least-squares adjustment, and those coming directly from 
observations, may be considered as representing simply errors of 
observation or as representing a variation of the apparent latitude 
(having the period m question) peculiar to the station, and not 
explainable by the motion of the pole. The existence of such indi- 
vidual variations of latitude with yearly period is well established. 
In so far as the residuals at all stations are alike their common value 
may be considered as the representation of the periodic annual por- 
tion of the Eamura term by harmonic constants. (See also p. 61.) 

The normal equations (31) fall into two groups of two equations 
each, the groups being identical in form out involving aifferent 
unknowns. 

With the longitudes of the 6 stations of the International Latitude 
Service the normal equations become numerically 



(35) 



+ 2.1433 m-0.1995p = ^, 
-0.1995 m + 3.8567 p = 5, 

+ 2.1433 71-0.1995 3= C, 1 ,«^s 

-0.1995 n+3.8567g = Z>, P"^^^ 

where for compactness the following abbreviations are used: 



A = [R cos f cos X] = [c cos X], 



c sin X 

8 COS X 



H37) 



B^[R cos f sin X] = 

C=^[R8m l cos X] 

D=^[R sin f sin X] = [5 sin X], 

The Aj of course, has no connectioit with the same letter used to 
denote the semi-major axis in equations (14), (18), (20), or elsewhere. 
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The solution of (35) and (36) gives 

m= +0.4688^ + 0.0243 B, \ ,oo\ 

2> = + 0.0243 A + 0.2605 B, J ^"^^^ 

n= +0.4688 (7+0.0243 D, 1 ,oqn 

g= +0.0243 (7+0.2605 Z>. J ^^^^ 

When only the three stations Miznsawa, Carloforte, and Ukiah, are 
used with equal weights equation (31) becomes numerically 



+ 1.8854 m-0.1128p = -A, 
-0.1128 m+1. 1146 p = 5, 

+ 1.8854 71-0.1128?= (7, 
-0.1128 n+ 1.1146 ? = Z). 



The solution of (40) and (41) gives 



m= +0.5336 ^ + 0.0540 5, 
p= +0.0540 ^ + 0.9026 J?, 

n= +0.5336 (7+0.0540 Z), 
2= +0.0540 (7+0.9026 D. 



}(40) 
}(41) 

}(42) 
}(43) 



By substituting the values of the individual X's in equations (38) 
to (43) the values of m, p, etc., may readily be expressed in tenns of 
the s's and c's of the individual stations. 

If it be assumed in advance that the eccentricity of the eUipse is 
zero — that is, that the pole moves imiformly in a circle — the adjust- 
ment is extremely simple. It has been noted (p. 10) that for motion 
of this sort the quantity f + X should be constant for all stations,* and 
the ampUtude, n, is clearly constant also. The process, then, is to 
take the arithmetic mean of the values of R at the several stations as 
the adjusted value of R and the mean of the several values of f +X as 
the adjusted value of f + X, and from it deduce the adjusted values 
of the rs. If there is considerable range in the values of R and f + X, 
and if in spite of this the assumption of circular motion is made, it 
is rather better to form the values of R cos (f +X) and R sin (f +X) 
at the several stations, to use the arithmetic means of these quantities 
as their adjusted values and from these adjusted values to deduce 
the adjusted values of R and f + X. It is not difficult to see that this 
latter process is the one that results from specializing for the circle 
tiie discussion just given for the more general case of the ellipse. 

It is not very difficult to obtain, on the assumption already made, 
the formulas for the probable errors of the elements of the ellipse 
deduced by formulas (l) to (VI) (pp. 85 and 86) from the adjusted 
values of m, n, p, and q. Since, nowever, the element of especial 
interest is the angle 6 that gives the direction of the major axis of 
the ellipse (see p. 86) we shall confine ourselves to deaucing for- 
mulas for the probable error of this quantity. 

For this purpose we need the quantities g— » g-» g-» and g- • From 

(IV) by expanding tan 27 and cos (a — 0) and expressing the result 
in terms of m, n, p, and q, by means of (II) and (20) we readily find 

tan2«=-4fn^±#^=^, (44) 

m^ + n^ — p^ — g^ D 

* This applies to the negative rotation of the annual or 14-month components; for positive rotation f— X 
is constant. 
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where N and D stand for the numerator and denominator of the pre- 
cedmg fraction. This D has no coimection with the D previously used 
for [s sin X]. By differentip,ting both sides with respect to m 

am D^ 

or substituting in the left-hand side for tan 20 and simplifying, 

bS ^ dm bm 

with similar expressions for ^» etc. 

The quantity N^ H- D^ may be expressed in terms of the auxiliary 
quantities already used and turns out to be equal to d* cos^ 5. We 

have g^ = 2p and ^= 2m and restoring for convenience the values of 
the auxiharies previously used, so that 

N=-2 ah cos (a-p) and D = a}-l^y 

we get 

d^ _ (fl^-6^) p~2 db m cos (a~i8) 

dm d* cos' 5 

_ (a'~6') r 2q6cos (a-/3) m "] 
d* COS'S L^ a' -6' J 

or by (7) and (12) 

d^ cos 27 r + ozn 

cos 2'v 
= d? cos' a COS 2 fP «°" 20 - m sin 2fl]. (45) 



where 



= Z (p cos 2^ - m sin 26) , (46) 

i7_ cos 27 _ 1 /47X 

^~d!' cos' 5 cos 2^ d' cos 5' ^^^ 

The second form of K may be derived by considering the spherical 
triangle mentioned on page 85. 
In a similar way we find 

g= Z (g cos 2e-n sm 26), (48) 

^-^K(m COS 26-\-p sm 26) (49) 

|^= K (n cos 2(? + g sin 26). (50) 

The independent quantities of which ^ is a fimction are not ulti- 
mately m, n, Pj and jcr, but the values of R cos f = c and R sin f =s 
in equations (32). The quantities c and s are treated as the inde- 
pendent observed quantities in the observation equations. 
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The formulas applying to such a case are given in almost anj com- 

f>rehensive work on the method of least squares.* In applying the 
ormulas, however, the probable error of one of the so-called observa- 
tion equations (32) must be taken as 0.6745 times the mean-square 
error from equation (34). If we call this probable error 6 and denote 
by €e the probable error of the adjusted value of ^, we find in this 

''' «--i(0«..-(^XI>"HI)"«' 

If ^ is to be in degrees instead of in radians, the factor (180°/x)* 
must be introduced into the right-hand side of (51). The quantities 
Qi\y Qnj fti^d Q22 are the coefficients of A and B in the solution of the 
normal equations. Thus, when aU six stations are used with equal 
weight, the Q's from equations (38) or (39) are as follows: 



+ 



(11 = + 0.4688, 
!i2 = + 0.0243, 
^22= -f 0.2605. 



(62) 



When only the 3 stations, Mizusawa, Carloforte, and Ukiah, are used 
with equal weight the Q^s from (42) and (43) are 



Gil =+0.5336, 

"!i2= +0.0540, 

(22= +0.9026. 



(53) 



The probable error of should not be affected by the axes chosen; 
that IS, if the x aijps were chosen in another meridian than that of 
Greenwich and if the longitudes of m, n, p, and 2> ai^d of the Q^s 
were computed for the new axes, these new values used in equations 
(46) to (51) should give the same value of te. This fact may be 
verified by direct transformation of no great intrinsic difficulty but 
of considerable length* in fact e^ is imchanged by transformation, and 
so is each line of (51) within the braces. Tms fact enables us to 
express ^ in terms of quantities depending on the dimensions of the 
adjusted ellipse itself, and not on the axes to which it is referred. 
The result, after introducing the factor (180/7r)^ so as to give ^ in 
degrees, is 

The square brackets here have no special significance as a sign of 
summation but are merely the usual symbol of aggregation. The 
quantities d and b are independent of the particular axes used, as is 
seen from equations (18) to (20), for it is easv to show that (P== 
A^'\-B^ andcP cos5 = ^^ —B^, The quantity witnin the braces is nu- 
merically invariant for a rotation 01 the axes through a given angle- 
in fact, the portions [QU + Q22] ^^^ [(di "622) cos 2^ + 2Qi3 sin 2^J 
will be foimd to be separately invariant. 

» For example, Wright and Hayford's Adjustment of Observations (2nd edition) p. 137, or Helmert's DU 
Ausgleidimigsrechnung nach der Methode der Kleinsten Quadrate (2d ed.) p. ISd 
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If we accept the hypothesis of a Kimura term that is constant at 
all stations, the natural procedure is then to combine the expression 
for the annual portion of the Kimura term with the expression for 
the variation oi latitude due to the motion of the pole. Let the 
annual portion of the Kimura term be representea by h cos kI 
-f i sin Kt and added to the right-hand side oi equation (23), it being 
tmderstood that only the annual part of the apparent variation oi 
latitude is under consideration. By a process similar to that on 
pages 87 to 89 we find, instead of the four '^normal" equations (31), 
six equations 



m [cos'X]H-p [sinX cos \] + h [cos X] = [c cos X]=-4, 
m [sin X cos X] + p [sin' X] + A [sin X] = [c sin X] = B, 
m [cos X] + p [sin X] + n h^[c]^E, 
n [cos' X] + a [sin X cos \] + l [cos X] = [s cos X] = C, 



n 
n 



M55) 



'sin X cos X] + q [sin' X] + Z [sin X] = [s sin X] = D, 
cos X] + g [sin X] + n Z = [s] = F. 



In (55) n is the number of stations and the square brackets indicate 
summation, as in the theory of least squares. 

These equations may be solved in terms of ^, B, E, etc. For all 
6 stations the solution, which corresponds to (38) and (39) is 






+ 0.4782 ^ + 0.0369 B -0.0412 E, 
+0.0369 ^ + 0.2776 5-0.0554 E, 
-0.0412 A -0.0554 B+ 0.1801 E, 



(56) 



(57) 



with similar equations for n, q, and I, namely, 

71= +0.4782 (7+0.0369 Z) -0.0412 F, 
2= +0.0369 (7+0.2776 Z> -0.0554 F, 
7= -0.0412 C -0.0554 Z? + 0.1801 F, 

When only 3 stations are used there are just sufficient conditions 
to determine the unknowns, but the least-squares type of solution 
may be retained for uniformity. For the stations ACzusawa, Carlo- 
forte, and Ukiah, we have corresponding to (40) and (41) 



m = 



r. 



+ 0.5442 ^ + 0.0516 B + 0.0600 E, 
+ 0.0516 .4 + 0.9032 B -0.0136 E, 
+ 0.0600 A -0.0136 5 + 0.3404 E, 



(58) 



71= +0.5442 C+ 0.0516 Z> + 0.0600 F, 
£= +0.0516 C+ 0.9032 D -0.0136 F 
7= +0.0600 (7-0.0136 Z> + 0.3404 F. 



(59) 



The values of m, n, p, g, Ji, and I in equations (56) to v59) may 
readily be expressed, if desired, in terms oi the individual s's and c's. 



SECTION 6. THE MATHEMATICAL THEORY OF 

LATITUDE. 



THE VARIATION OF 



An exhaustive treatment of the theory of the variation of latitude 
is not attempted here, but merely an outline, chiefly of the aspects 
most important in the present discussion. A number of references 
to the literature of the subject are given and these will supply the 
omitted proofs of the various statements here made, 
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The simplest mathematical problem that has sufficient resemblance 
to natural conditions to merit discussion is that of a rigid body set 
rotating but subject to the action of no forces. This problem was 
first adequately treated by Euler; • a more geometrical treatment 
was given by Poinsot.' The subject has aroused the interest of 
many mathematicians and will be found fully treated in standard, 
works on theoretic^ mechanics.® It was to this theory that astrono- 
mers looked for a possible periodic variation in latitude, and it was 
still appealed to alter the existence of such a variation had been 
estabhshed by observation but before the natiu*e of the variation was 
as well imderstood as it is now. Moreover, the simplified form of 
the equations of this theory can be made to apply to the actual case 
of a nonrigid body merely "by changing the numerical values of some 
of the constants. For present purposes the rotation may be suflB- 
ciently specified by referring it to rectangular axes fixed in the body 
but not fixed in space. The question of the position in space requires 
a separate treatment. Let the principal moments of mertia of the 
body be -4, B, and C in ascending order of magnitude, and let the axes 
of Xj y, and z coincide respectively with the axes of A, J?, and C, and 
let p, 2, and r denote respectivelv components of the rotation about 
these axes. For the case of no forces Euler 's equations are: 






(1) 



The general solution of these equations can be expressed in terms of 
Jacobi's eUiptic fimctions. sinam u, cosam u, etc., but for any condi- 
tions resemoling those oi nature the modulus is so very small that 
the eUiptic fimctions may be replaced by the ordinary trigonometric 
functions. It is more convenient to make the reauisite simplifying 
assumptions at the start and this we shall do, altnough we thereby 
lose the advantages of greater generality and of ease, in estimating 
the effect of omitted terms. 

In applying equations (1) to the earth the northerly direction of 
the z axis wilTbe treated as positive. The z axis or axis about which 
the moment of inertia is (7, the greatest of the principal moments, is 
the so-called axis of figure, determined by the form of the earth, or 
rather by its moments of inertia^ it is to be distinguished from the 
continually shifting axis of rotation, the direction of which is given 
by the instantaneous value of the resultant of the component rota^- 
tions p, J, and r. Since the latitudes vary but little, the axis of 
rotation must never depart greatly from the axis of figure; that is, 
the component rotation r about tne axis of figure must represent in 
magnitude nearly the entire rotation of the earth and p and gr must 
be small® in comparison with r. The axes of x and y are m the 

• Theoria Motas Corporain Bolidonmi sea RMdomm. Greifswald. 1766; second ed. Oreifswald, 1190, 
Chaps . xn and XIH. A possible application of the mathematics to the variability of terrestrial latitudes, 
at that time of course unknown to observation, is mentioned in Chap, xn, Prob. 74, Bcholiimi 3 (p. 29d 
of the edition of 1790). 

1 Th^orie nouvelle de la rotation des corps. Paris, 1851. partlcularlypart 2, pp. 65-122. 
> For example, Routh's Advanced Higla Dynamics, flf tn ed., p. 88. See also Helmert; Hdhefe Oeod&ie, 
Vol. II, p. 390, or Appel, Trait6 de M^canique Rationelle, 3rd Ed. Vol. II, p. 162. 

* The ratios p : r or 9 : r would be quantiaes of the order of magnitude of A in radians, A^ being tbe 
variation in latitude from the mean, that is, quantities of the order of one part in a million. 
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equator of figure, the x axis beii^ the axis about which the moment 
CI inertia is the less or along which the radius of the geoid is the 
, 'eater. (See fig. 12.) It is convenient to make the rotation about 
he z axis, whidi is the principal portion of the earth's rotation, 
positive; tnat is, from the positive du*ection of x toward the positive 
direction of y. If we take the axes of x and y accordingly, and com- 
pare these axes with the axes resulting from adding a north-pointing 
z axis to the x and y axes of the International Latitude Service 
(p. 33), we see that the two sets of axes can not be made to coincide, 




Fig. 12.— Coordinate axes for geometrical treatment of the motion of the North Pole and for the dynamical 

problem of the variation of latitude. 

Axes are of same nature as those of International Latitude Service except that the x axis is in the me- 
ridian plane of least equatorial moment of inertia instead of in meridian plane of Greenwich. In the dy- 
namic problem the axes of x and y are parallel to those in figure with ori^ at the center of the earth and 
with the positive direction of the y axis reversed as coi^ipared with the figure. 

F^pole of inertia. 

FMin angular units»x<-p/a>. 

MP in angular urdts^^^q/ta. 

but that if the positive axes x and z coincide in the two sets, the axes 
of y are oppositely directed. This fact must be remembered in pass- 
ing from me axes here used for the problem in mechanics to the axes 
of the International Latitude Service. 

Since © and q^ are so small we neglect their product in the fitrst of equa- 
tions (1). This neglect is further justified by the fact that B—Aia 
small for the earth as compared with O—B or C—A in the other two 
equations of the set. If A =5, the process of dropping the term in 
pq is of course rigorous. From the nrst of the equations (1) we find 



or 



r= const. =« = angular velocity of earth's rotation. 

89016*'— 22 7 



(2) 
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Substitute from (2) in the second of (1) and differentiate with respect 
to t, the time; the result may be written 

Substitute for ^ in (3) from the third of (1) and we get 

d'p . iC-A) (C-B) ,^ ^ 
-J+ -^ «»p=0. 

or 

^+Pp=0 (4) 

where P is written for -j^ co^, a quantity necessarily 

positive. The general solution of (4) is 

p = h cos (kt-\-fi)y (5) 

where h and P are constants of integration, h being restricted to small 
values on accoimt of the assumptions made. The corresponding 
value of a comes from the second of (1) and is 

_ ~A dp 

*~((7-5)co dV 
or 

The direction cosines of the resultant axis of rotation are in general 

VF#+?' W+F+^' Vp' + g' + r" °' '^'' '^' '^'^^'^ °^ ^ ^'^^ * 

are to be neglected the direction cosines are-? -, 1. To state the 

CO O) 

matter otherwise, the quantities ^; • are the sines of angular distances 

from the pole of rotation to the planes zy and zx, and being small 
may be taKen as the angular distances themselves in radians. When 
expressed in seconds of arc and projected orthogonally, on a plane 

tangent to the earth at the pole of figure, 2, - are the rectangular 

components of the variation of latitude of the same sort as are dealt 
with in section 4, though referred to a different set of rectangular 
axes, in that the x axis in section 4 is in the meridian plane of Green- 
wich and not, as here, in the meridian plane of the axis of the smaller 
principal equatorial moment of inertia. (See fig. 12.) 

From equation (5) or (6) the period of the variation of latitude is 
seen to be 27r/it. If we write out the value of t in full and take the 
sidereal day as unit of time, so that w = 27r, we get for the general case 

Period in sidereal days = ^ j-j-^ — .. ,^ — ^tt- (7) 

V \0—A) kC—B) 

When ^ = 5, we get • 

Period in sidereal days = ^_ i • (8) 
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According to De Sitter's*® values the period from (8) is 304.1 
sidereal days or 303.3 mean solar days; that is, about 10 months. A 
table of periods deduced from (7) is given below. This corresponds 
to the free period of the variation of latitude, a period dependent 
only on the properties of the body itself and not on any external 
forces. For a rigid earth the more nearly spherical it is the longer 
would be the free period. The difference between this free period 
and the observed free period, 14 months (432.5 mean solar days), 
is dealt with on page 98. 

If we consider the motion of the pole as given in rectangular 
coordinates, x=^p/o)y y^^^Q/oo, p and q being taken from (5) and (6), 
we see that the motion of the pole is of the harmonic elhptic variety 
treated in section 4. This means that in its free motion, tJie pole of 
rotation describes in the body an ellipse around the pole oj figure. The 
semiaxes of the el lipse, whic h are along the x and y axes, are, respec- 
tively, b and b \ //T^x ^ ' I^ is easy to see that the second quan- 
tity is larger than the first. The longer axis of the ellipse of polar mo- 
tion is in the meridian of the y axis, that is, in the meridian of the axis 
of the larger of the two principal equatorial moments of inertia or the 
meridian of the shorter equatorial axis of the geoid. This is the same 
condition as prevails in the more general case referred to later in this 
discussion and also on pages 62-64. 

For definiteness let us assume the quantity C— i{A + B) to be 
fixed and compute the ratio of the axes of the latitude variation 
ellipse and the period of the variation for different values of the 

B—A 
ratio p_ x(A-\-B) (^^^^^ / ^^^ simplicity). 

The results are shown in the table that follows: 



- 


Ratio of Axes 

liC-A) A 
M{C-B) B' 


Period. 


C-^(A+B)°^^' 






Mean Solar 






Days. 





1.000 


303.3 


1/100 


1.005 


303.3 


1/50 


1.010 


303.3 


1/46 


1.011 


303.3 


1/25 


1.020 


303.3 


1/10 


1.051 


303.7 


1/6 


1.105 


304.8 


1/4. 


1.133 


305.7 


1/3 


1.183 


307.6 


2/5 


1.224 


309.5 


1/2 


1.290 


313.2 


3/5 


1.361 


317.9 


2/3 


1.413 


321.7 


3/4 


1.481 


327.2 


4/6 


1.526 


330.9 


1 


1.729 


350.2 


6/5 


1.996 


379.1 


7/5 


2 375 


424.7 


8/5 


2.992 


505.5 


17/10 
9/5 


3.502 


575.7 


4.346 


695.7 



10 



On Isostasy, the Moments of Inertia and the Compression of the Earth. Koninklikje Akademie van 
^enschappen te Amsterdam. Proceedings of the meeting of Apr. 23, 1915 (Vol. XVII). 
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The fourth line in the table, argument 1/46, corresDonds to the 

B — A 
value of fi^i(AjL.'Q\ deduced from constants of Helmert's adopted 

gravity formula already referred to (p. 63). The corresponding 
ratio of the axes, 1.011, gives an ellipticity of 0.010, or an eccen- 
tricity of 0.15, which is less than the ellipticity deduced by Schweydar 
for a nonrigid earth, but the effect of an inequality in the equatorial 
moments is of the same general nature for a rigid earth as for a 
nonrigid earth. 

The observed free period of latitude variation of about 14 months 
was very puzzling to the early students of the variation of latitude, 
who were confidently expectmg a period of about 10 months, and 
could see no plausible reason for any other period. To be sure, an 
inequality in the two principal equatorial moments, A and 5, would 
cause a lengthening of the period but this inequality must be verv 
pronoimced indeed to produce any appreciable effect on the period, 
as may be seen from an inspection of the last column of the table. 
A period of 432 days falls between the arguments/=7/5 and /= 8/5- 
any such inequality in the moments of inertia would make itseli 
felt very emphatically in geodetic operations, both in the values of 
gravity and m the deflections of the plumb Une.^^ The corresponding 
ellipse of polar motion comes out much elongated instead of the 
almost circular eUipse actually observed. The difference between 
10 months and 14 is then not due in any appreciable degree to the 
inequality between the equatorial moments of inertia. The now 
generally accepted explanation of the lengthening of the period was 
given by Newcomb,*^ and the matter may be presented very clearly 
m his own words: 

Mr. Chandler's discovery gives rise to the question whether there can be any defect 
in the theory which assigns 306 days as the time of rotation. The object of this paper 
is to point out that there is such a defect, namely, the failure to take accoimt of the 
elasticity of the earth itself and of the mobility oi the ocean. 

The mathematical theory of the rotation of a solid body, on which conclusions 
hitherto received have been based, presupposes that the body is absolutely rigid. 
As the earth and ocean are not absolutely rigid, we have to inquire whetiier their 
flexibility appreciably affects the conclusions. That it does can be shown very 
simply from the following consideration: 

Imagine the earth to be a homogeneous spheroid, entirely covered by an ocean of 
the same density with itself. It is then evident that, if the whole mass be set in uni- 
form rotation aroimd any axis whatever, the ocean will assume the form of an oblate 
elKpsoid of revolution, whose smaller axis coincides with that of rotation. Hence, the 
axes of rotation and of figure will be in perfect coincidence imder all circiunstances. 

To apply a similar reasoning to the case of the earth, imagine that the axis of rota- 
tion is oisplaced by 0^20 from that of greatest moment of inertia, which I shall call 
the axis of figure. Then, with an ocean of the same densitv as the earth, its equator 
would be displaced by the same amount. The ocean level would change in middle 
latitudes by about one inch at the maximum. But this change would have for its 
effect a corresponding chang[e in the axis of figure. As ihe ocean covers only three- 
fourths of the earth the axis would be displaced by three-fourths of the distance 
between the two axes were ocean jand earth of equal densitv. But as the density of 
the earth is some five times as great the actual change woula be only one-fifth of this. 
It would even be less than one-fifth, because the displacement of the ocean equator 
would be resisted by the attraction of the earth itself. The exact amount of this 
resistance can not be accurately given, but I think the displacement would thereby 

11 For /= 7/5 the inequality in gravity at the Equator would be over 1.1 centimeters above ot below its 
mean vuue, a variation from minimum to maximum about as great as the mean change in gravity between 
Equator and pole apart from the effect of the centrifugal force. The maximum deflection of the vertioEil 
would be about 8 nunutes of arc. 

i> '' On the djmamics of the earth's rotation with respect to the variation of latitude," Monthly Notices 
of the Royal Astronomical Society, vol. 52 (1892), p. 336. 
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be reduced to one-half. I therefore think that one-fourteenth would be an approxi- 
mate estimate of the displacement of the axis of figure *^ in consequence of the move- 
ment of the ocean. As Mr. Chandler's period requires a displacement of two-sevenths 
the ocean displacement only accounts for one-fourth of the difference. 

The remainder is to be attributed to the elasticity of the earth itself. It is evident 
that the flexure caused by the noncoinddence of the two axes tends to distort the earth 
into a spheroid of the same form as that which the ocean assumes, and thus to bring 
the two axes together. 

We have now to show how this deformation of the earth changes the time of revolu- 
tion. Let us imagine ourselves to be looking down upon the North Pole and let P 
be the actual mean pole of the earth when the two axes are in coincidence and R the 
end of the axis of rotation. Then, in consequence of the rotation around iJ, the 
actual pole will be displaced to a certain point, P^. Now the law of rotation of i2 is 
such that it constantly moves around the instantaneous position of P^ itself. In 
other words, the angular motion of R at each moment is that which it would have if 
P^ had remained at rest. Hence the angular motion as seen from P is less than that 
from P' in the i-atio of P'R .- PR. 

But as R rotates, P^ continually changes its position and rotates also, remaining 
on the straight line PR. Thus the time of revolution of R around P is increased in 
the same ratio. 

Newcomb's argument may be illustrated by figure 13. Let B and 
7?i be ''consecutive'' positions of tbe pole oi rotation, separated 
by an infinitesimal interval of time dt The corresponding positions 
of the actual pole of figure, as affected by the yielding of the earth, 
are P' and P/ lying on the lines joining u and R. with P^, the point 
where the pole of figure would«be if the pole oi rotation were to 




FiQ. 13.— Relation between the undisturbed pole of figure, the actual pole of figure, and the pole of rotation. 

coincide with P. Let e be the angular velocity of the point R about 
P' and K its angular velocity about P. Then for the infinitesimal 
distance RR^ we have 

RR, = {PR)Kdt = {P'R)€dt 

^^ k:€=P'R:PR 

as stated by Newcomb. There are thus three different poles that 
appear in the theory of the variation of latitude: (1) The pole of the 
axis of instantaneous rotation; this is the pole that is located by the 
astronomical observations of the International Latitude Service; in 
figure 13 it is the point R or R^, (2) The point that would be the 
pole of figure if the pole of rotation should nappen to coincide with 
it. This pole may be called the undisturbed pole of jiaure or undis- 
turbed pole of inertia. The requirement that the pole of rotation 
should coincide in imagination with the pole of figure is necessary in 
order to be rid of the centrifugal force arising from the noncoinci- 
dence of the two poles. See under pole (3). The changes in the position 
of this pole are due solely to changes in the distribution of matter 
within and on the earth apart from any stress arising from the rota- 
tion. It is represented by P in figure 13. (3) Between P and R is 

" See p . 100. Without the effect of the yielding of the ocean the free period of the latitude variation would, 
•ocording to the estimate, be some 40 days less than the observed value. 
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the actual pole of figure, represented by P' or P/. The displacement 
of the actual pole away from the undisturbed pole is due to the 
slight centrifugal force resulting from the noncoincidence of the pole 
of rotation wiui the imdisturb^ pole of figure and to the yielding of 
the earth imder this centrifugal force, a vidding that tends to deform 
the earth so as to displace the actual pole of figure nearer to the pole 
of rotation. 

The equations of motion of a slowly changing body have been 
treated by various authors." Convenient expositions and develop- 
ments of their work will be f oimd in various textbooks and treatises.** 
There have been several treatments with special reference to the 
variation of latitude; one of the most recent and thorough is due to 
Schweydar.^® For our purposes use will be made of an article by 
Larmor ^^ rather too long to reproduce here, in which he shows by 
verv general reasoning how the problem of the motion of a yielding 
body may be reduced to that of a rigid one. Larmor proves that 
for the actual elastic earth we may suDstitute a rigid earth of about 
the same volume but with so much of its eauatorial bulge removed 
as would be due to the yielding of the eartn under the entire cen- 
trifugal force of rotation.^^ This would leave the earth less flattened, 
and since the differences of the moments C-A and C-B serves to 
measure the flattening, we must replace these differences by some- 
thing smaller by about 3/10 of itseli. While the relative changes in 
the differences of the moments are considerable, the change in any 
one moment is relatively small when compared with the moment 
itself. The smaller flattening of the equivalent rigid body corre- 
sponds to a longer free period. The free period is thus about 
303/(1 - 0.3) = 433 mean solar days. 

A similar calculation on the ratio of the axes of the eUipse of polar 
motion shows that the yielding of the earth makes the ellipse flatter. 
The effect of removing the centrifugal force is represented by a 
quantity proportional to PjCsin <^), where Pg is the iWendre's zonal 
harmonic of the second degree and <^ is the latitude. The inequahty 
in the two moments A and S may be represented by another surface 
spherical harmonic and the two harmomcs may be superposed when 
neither is large, which is the supposition here made. The quantity 
to be deducted from each difference of moments 0—A and C—B is 
then the same for small differences between A and B. The rigid body 
equivalent to the yielding body will then have for a ratio of axes of 
the polar eUipse 



l {0-A-x)A 
\{C-B-x)B 



14 Liouville: Additions k la Connaissance des Temps pour 1859, or Journal des Math^atiques puies et 
applique, vol. 3 (1858). Gylden: Recherches sur la rotation de la Terre; Actes de la Soci^t^ royale des 




»» Routh: Advanced Rigid Dynamics; 6th ed., p. 17. Lamb: Hi^er Mechanics (Cambridge, Ens., 
192UXP.171. Tisserand: M^canique C^este; VoL il, Chap. XXX. fielmert: H6here Gtood&sie; voiTsi, 

u Die Bewegung der Drehachse der Elastischen Erde in Erdkdrper und im Ranme: A s trouomiaohe 
Nachrichten, voL 203 (1916). p. 101. 

» The relation of the eartn^s free processional nutation due to its resistance against tidal deftnmatioo, 
Ftoceedings of the Royal Society of Ixmdon, sec. A, vol. 82 (1909), p. 89. 

18 If the entire centrifugal force were removed, it is to be presumed that the earth would in time ykld 
plastically and take on a nearly spherical form. This, however, is not what is meant. The yielding is 
to be taicen as elastic, that is, proportional to the force and not increasing with the time, except forthe 
yielding of the oceanic waters already mentioned. The effect of the rearrangement of the water on the equa- 
torial bulge, or more precisely on C—Ay must be allowed for in estimating the purely elastic yielcung. 
The material of the earth being what it is, if the cenlarifngal force were removed, the equatorial bulge woufi 
be reduced by about 3/10 of its present value. 



" • 
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instead of 



V 



{C-A)A 
(C-B)B' 



The changes in A and B themselves are neglected in comparison with 
the change in C—A and C—B. For the earth the quantity x is 
about 3/10 of the quantity i {A + B). If the ratio f=-(B-A)/ 
[0'-i(A-\'B)] is equal to 1/46, as deduced by Helmert f rom gravity 
observations, then the ratio of the axis for a rigid earth is 1.011 as 
already given in the table on page 97 and for a yielding earth 1.016; 
tie corresponding eccentricities are, respectively, 0.15 and 0.18. This 
result of this rough and ready calculation agrees with Schweydar's 
result as based on a more elaborate method of calculation. 

The reduction of 0.3 for the yielding of the earth includes the effect 
of the yielding of the ocean waters. On account of the irregular dis- 
tribution of land and water on the globe this effect is not quite sym- 
metrical about the earth's axis of rotation; this dissymmetry affects 
the ratio, the shape of the ellipse of polar motion, and also tne direc- 
tion of its major axis, so that the meridian in which the latter lies 
would no longer coincide with the meridian of the larger principal 
equatorial moment of inertia. The problem of correcting for the 
irregular distribution of land and water on a rigid earth is somewhat 
similar to the problem of correcting the equiubrium theory of the 
tides for the distributio^ of land and water, and has been practically 
solyed by Darwin and Turner" in accordance with the ideas of Lord 
Kelvin. The work of Darwin and Turner shows that the equilibrium 
tides at any point in the ocean are, almost the same as if the earth 
were entirely covered with water and a preliminary computation for 
the analogous problem of the effect of tne arrangement of land and 
water on the motion of the pole shows that this effect is compara- 
tively small also. The problem is further complicated by the elastic 
yielding of the earth under the irregularly distributed pressure of the 
ocean water. The results on page 60 are, of course, uncorrected. 

The general problem has been treated by BriU^® but his results do 
not appear to be immediately appUcable to the problem in hand. 
Brill does, however, reach the conclusion that Newcomb's estimate 
(p. 99) of the effect of the yielding of the ocean water on the free 
period of latitude variation, thou^ confessedly a rough estimate, 
is in close agreement with Brill's own more elaborate calculations. 
An attempt will be made to evaluate the correction for a distribution 
of land and water more closely conforming to the actual one than the 
one used for simplicity in the preliminary estimate, and also to estimate 
the effect due to the elastic yielding ox the earth under the irregular 
distribution of water pressure. The results, when applied as correc- 
tions to the quantities specifying the shape and position of the 
ellipse of free polar motion, should render tne results as to the mo- 
ments of inertia of the earth derived from the corrected ellipse 
directly comparable with those from gravity observations. 

When the moments of inertia are imequal the actual pole of figure 
(P' in figure 13) is no longer exactly on the line RP joining the 

>* "On the correctiom of the equilibrium theorv of tides for the oontinents, " Proceedings of the Royal 
Society of London, vol. 40 (1886), p. 303, or Darwin's Scientific Papers, Vol. I (Cambridge, England, 1907, 
p. 328). The theory does not treat the self-attraction of the water, or rathw it supposes the seu-attraction 
to bethe same as if the earth were completely covered with watw. 

M Uber die Elastizitat der Erde (doctor's thesis), Gdttingen, 1908. 
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undisturbed pole of figure with the pole of rotation. Schwe^dar 
computes that when the pole of rotation describes an ellipse having a 
ratio of axes equal to 1.016, the actual pole of figure describes an 
ellipse in which the ratio is 1.038, a somewhat flatter ellipse than 
that described by the pole of rotation. 

So far the discussion has been of the free vibration only. Owing 
to changes taking place in the distribution of matter in the earth 
the undisturbed poles of figure shift slightly. The most conspicuous 
changes in the distribution of masses are approximately periooic with 
an annual period, and are due to changes oi Uke period in barometric 
pressures, rainfall, ocean currents, etc. 

Corresponding to terms of any period in the motion of the imdis- 
turbed pole of inertia there will be terms of like period in the motion 
of the pole of rotation. The portion of the polar motion expressed 
by terms of the same period as those that express the motion of the 
undisturbed pole of figure are said to represent the forced motion 
of the pole of rotation. The form of path represented by the motion 
in the free period is always the same, that is, it is always an ellipse 
of fixed eccentricity lying in a fixed direction or a circle as the case 
may be, regardless of the form of the forced motion, but the size of 
the ellipse or circle and the initial phase of the free motion do depend 
on the forced motion. 

When a term is present in the forced motion with a period approxi- 
mately equal to the free period, the familar phenomenon oi reso- 
nance^' occurs and the amplitude of the free motion is greatly increased. 
Hie annual term, which represents the principal part of the forced 
motion in the case of the earth, is near enough in period to the 432.5 
days of the free period to cause a considerable increase in the ampli- 
tude of the free motion. The precise amount of the increase depends 
on initial conditions, but the ampUtude of the free motion is prooably 
five or six times as great as it would be if the forced motion has a 
much longer period. It is an interesting speculation to consider 
whether, at some time in the past, the free period owing to a diflFerent 
flattening of the earth and different internal conditions might not 
have approximated closely to a year. The free motion of the pole 
would tnen become large, and various peculiar conditions might 
conceivably arise which, for a yieldiug body, can not be followed oy 
the present theory, since the latter is limited to small displacements 
of the pole." 

The equations connecting the motion of the pole of rotation with 
that of the undisturbed pole of figure are stated below without proof. 
Let X and y be the rectangular coordinates of the pole of rotation 
projected on a plane tangent to the earth at some point near the mean 
pole, and let | and t\ be the rectangular coordinates of the imdis- 
turbed pole of figure referred to the same axes. The equations in 
question are then , 



I=-*-e- 



(9) 



B For a li^ body the solution of eqoatioiis (1) in tenns of the elliptic functions (p. 94) is perfectly general. 
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Here K^-p, P being the free period expressed in ihB same imit as 

the time t It should be said that the axes used are similar to those 
of the International Latitude Service used in Chapter IV in that the 
y axis is so taken that the rotation of the earth as seen jErom the North 
tole and the free motion of the pole both involve negative rotations. 
The direction of the x axis, however, is not necessarily toward Green- 
wich as it is for the x axis of the International Latitude Service. It 
should be said also that equations (9) apply to the case of equal 
equatorial moments of inertia, that is, to a curcular free path of the 
pole of rotation when the forced motion is periodic. The more 

general case of imequal equatorial moments does not appear to have 
een fully ii^orked out for the forced motion. The expressions for f 
and 17 are not, however, limited to periodic terms but may be of any 
nature whatever. The solution of the two equations (9) simultane- 
ously, supposing J and 17 given, evidently leads to a linear differential 
equation of the second order with constant coefficients and a known 
second member. Equations (9) are fully discussed by Wanach^* 
who gives three forms of the general solution and considers a number 
of curious special cases that arise from different assumed forms for 
f and 17. The two examples that follow are adapted from Wanach. 
Example 1. — Suppose tnat the undisturbed pole of inertia has a 
progressive motion given by 

f = fo + g,^ 1(^0) 

Evidently Jo ^^^ ^0 ^^ ^^^ coordinates at the origin of time and a 
and b the velocities along the coordinate axes. The motion of the 
pole of rotation is given by 



x = ^o — Q+<^^ + ^ cos(irf— f), 



a 



where 

R cos f = a;o-?o + r' 



(11) 



fisinf = yo-^o-r- 



(12) 



The terms at and ht in (11) indicate that the progressive motion of 
the pole of rotation takes place in the same direction and at the same 
rate as the progressive motion of the undisturbed pole of figure. 
The free motion given by the periodic terms in ic^ is evidently uni- 
form and circular, the ampUtude and phase depending on the initial 
conditions. The effect of the progressive motion on ampUtude and 
phase is shown by the terms h/K and a//c in (12). 

Let us suppose a progressive motion equal to that derived in Chap- 
ter II, namely, 0f005 a year. The year being the unit of time, 
«=27r/1.184. For definiteness we take the motion as ^tirely alon^ 
the X axis, so that a = Of 005 and 6 = 0. The terra— a/K is then equ^ 

to 0f00094. This quantity is small in comparison with an amphtude 

. . — ■ — ■ — ■ ' 'I ' ■ 

s> Die Chandlersche uad die Newcombsche Periode der Polbewegmig; Zentralbureau der Intenuitioiialfin 
Erdmessung, publication no. 34, Berlin, 1919. 
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of 052000, which is about that of the free motion. The maximum 

effect on the amplitude is ± Of 0009 and on the phase about one-third 

of a degree. 

Example 2. — Suppose a forced motion of the undisturbed pole of 

inertia given by . , ^ v . 

^ '^ f=a cos (m<— a), 1 /iQX 

i7«6cos {mt-fi). I ^^^' 

Equations (13) represent a general elliptic motion of the type dis- 
cussed in section 4 of this appendix, tne period being 2irlm. The 
values of x and y are 



— K 



aj=— 2— 3 [Kaco&(m^— a) + m6sin (m^— i8)] + i?cos (<rf--f), 

ifv n 

y°» 2^ji [«6co8(m<— i8)--masin(m^— a)]— fisin (irf~f), 
where 

R cos f =Xo+ 2_,^ ('^^ cosa— mJsin/S), 

R sin f =yoH — 2~;5 ('^^ coBP+masiaa). 



(14) 



(15) 



Equations (15) show the effect of a possible ''resonance" between 
the periods of the free and the forced motion. If the periods approach 
equality, irF—t? decreases toward zero and the amplitude frona (15) 
increases indefinitely. 

The pole whose motion is followed by means of the observations 
of the International Latitude Service is the pole of rotation. When 
its coordinates x and y are known the coordinates of the undisturbed 
pole of figure may be deduced from (9). We get 



1^2/ 
, 1 dx 



(16) 



The values of f and r\ thus depend on the derivatives of x and y, 
and when x and y are deduced directly from observation and are 
therefore affected oy errors of various sorts, the derivatives deduced 
from the erroneous x and y are usually still more \mcertain and 
erroneous. Equations (16) give unsatisfactory results when they 
are used to foUow the motions of the imdisturoed pole of figure by 
using individual observed values of x and y. About all that can be 
done is to get the mean periodic portions oi ? and rj by using instead 
of the actual values of x and y, their smoothed out values obtained 
from harmonic constants.^"* If desired, instead of the constants 
themselves, the 12 or 24 mean results could be used that were put 
through the process of analysis in order to obtain the harmonic 
constants. 

A numerical illustration of the relation between the pole of rota- 
tion and the undisturbed pole of figure is given by figures 14 and 15. 
The harmonic constants of the pole of rotation are taken from 
Table 12 below the horizontal line and represent the results when the 

n The higher hannonio terms generally have little significance in this connection. (See p. 108.) 
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Kimura term is included in the adjustment. The semiannual terms 
are taken from page 58. When only the annual terms are included 
both poles describe ellipses, which are shown in full lines. The 
smaller and flatter ellipses belong to the undisturbed pole of figure. 
The major axes of the two elUpses are at right angles; it is easy to 
show that this is true in the general case. The roman numerals 

I, II, III, etc., indicate the position of the poles at the beginning of 
each twelfth of a year, or approximately on the first day of January, 
February, etc. 

If the harmonic constants of the pole of rotation used in the two 
figures be compared, it will be seen that the chief difference is in the 
value of /3. This changes the direction of the major axis by some 32®, 
but the ellipses are so nearly circular that the general cnaracter oi 
the motion of the pole of rotation is little affected. The effect on the 
path of the pole of figure is much more noticeable, since the ellipses 
are flatter; moreover the direction m which the pole of figure de- 
scribes its elUptic orbit in figure 14 is opposite to what it is in %ure 15. 

The results of including the semiannual terms is shown by the 
interrupted lines in the figures. The effect on the path of the pole of 
rotation is comparatively small, as the paths stUl resemble ellipses; 
on the path of the pole of figure, however, the effects are relatively 
great and there is little resemblance to the original ellipses left, nor 
do the two dotted curves for the path of the pole of figure at all re- 
semble each other. 

For convenience the formula by which the curves were plotted are 
given here. For figure 14 (series 1900-11, 6 stations): 
«= -0r0390 cos ^- Of 0788 sm ^ + 0.^0011 cos 2^ -Of 0012 sin 26. 
y= -0.^0738 cos ^ + 0f0194 sm 0- Of 0038 cos 2d + Of 0004 sm 26. 

The resulting values of the coordinates of the undisturbed pole of 
figure are from (16) : 

{ = - Of 01 60 cos 6 + Of 0086 sin 6 -f Of 0020 cos 26 -f Of 0078 sm 26. 
rj= +0f0195 cos d-0f0268 sin d-OfOOlO cos 2^ + 0.^0030 sin 26. 

For figure 15 (series 1900-17, three stations) : 
«« -Of 0385 cos d- Of 0774 sm d-|-0f0045 cos 2d -Of 0007 sin 26. 
y= -Of 0675 cos d-f pf0354 sm d + 0.^0007 cos 2d + 0f0007 sin 26. 

{ = + Of 0034 cos d + Of 0025 sin d + Of 0062 cos 26 - Of 0024 sin 26. 
n= + Of 0242 cos d-0f0102 sin d + 0f0024 cos 2d+0f0114 sin 26. 

The quantity 6 represents ^^, ',t being the number of days from the 

beginning of the year. The values of 6 corresponding to the points I, 

II, III, etc., are respectively 0°, 30'', 60^, etc. 

The numerical coeflS^cients of cos 6 and sin d in the expressions for 
X and y are of course, the values of a cos a, a sin a, etc., corresponding 
to the values of a, a, 6, and jS, in Table 12, or to put tne matter more 

Srecisely, the coefficients in question are values of m, n, p, and g, as 
educed by an adjustment oi the harmonic constants of the annual 
component at the several stations. The figures in Table 12 were 
deduced from the values of m, n, p, and g, given above. The rectan- 
gular coordinates, x and y, of the mean annual path of the pole of 
rotation over a given period may be accurately represented by taking 
a sufficient numoer of terms of two Fourier series in d, but tms repre- 
sentation may be valid in a mathematical sense only and then only 
for X and y themselves. The derivatives of x aad y with regard to 



